There are quite a few examples of code around that convert between different timescales. Probably the best examples are the IDL source code here. Note that BJD is barycentric JD which, for most VS purposes, is the same as HJD. I think the difference only matters if you are trying to detect planets going around pulsars…
British Astronomical Association
Supporting amateur astronomers since 1890