;+
; NAME: ircal_getj2000 
; PURPOSE:
;	Convert the current-epoch IRCAL coords into J2000 coords
;
; INPUTS:
; 	hdr			a FITS header
; KEYWORDS:
; OUTPUTS:
; 	ra, dec		J2000 coords.
;
; HISTORY:
; 	Began 2006-06-21 01:33:34 by Marshall Perrin 
;-

PRO ircal_getj2000,hdr,ra,dec

	ra = ten_string(sxpar(hdr,"RA"))*15
	dec = ten_string(sxpar(hdr,"DEC"))

 	; IRCAL's native RA and DEC keywords are in the current epoch at time of
    ; observation.
    date = sxpar(hdr,"DATE-OBS")
    datev = date_conv(date,"V")
    epoch = double(datev[0]) + (datev[1]+(datev[2]+(datev[3]+double(datev[4]))/60.)/24.)/365.24
	print,"input coords are "+adstring(ra,dec)
	print, "input epoch is "+strc(epoch)
	precess,ra,dec,epoch,2000.0,/print

end