function jd2epoch, jd
;+
; NAME:
;	JD2EPOCH
;
;
; PURPOSE:
;	Convert Julian Day to SMEI epoch.
;
;
; CATEGORY:
;	Utils
;
;
; CALLING SEQUENCE:
;	epoch = jd2epoch(jd)
;
;
; INPUTS:
;	jd	double	The julian day number(s) to convert
;
;
; OUTPUTS:
;	epoch	double	The SMEI epoch(s) 
;
;
; MODIFICATION HISTORY:
;	Original: 25/3/03; SJT
;-

jd0 = doy2jd([2000, 1, 12])     ; Epoch is noon on 1/1/00

return, (double(jd)-jd0)*8.64d4

end