;+
; NAME:
;	smei_mkstdstar
; PURPOSE:
;	Convert grd files for standard stars to Fits files
;	used by smei_skyd
; CATEGORY:
;	camera/idl/star
; CALLING SEQUENCE:
	PRO smei_mkstdstar
; INPUTS:
; OPTIONAL INPUT PARAMETERS:
; OUTPUTS:
; OPTIONAL OUTPUT PARAMETERS:
; INCLUDE:
	@compile_opt.pro		; On error, return to caller
; CALLS:
;	who_am_i, grd_read, IsType, gzip_file
; MODIFICATION HISTORY:
;	JUN-2006, Paul Hick (UCSD/CASS)
;	SEP-2006, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
;		In addition to old PSF there now is an alternative PSF
;		(presumably with the background removed more accurately)
;		The old standard_star.fts is now standard_star_fixd.fts;
;		the new PSF is in standard_star_fill.fts.
;-

root = who_am_i(/dir)

psf = ['fixd','fill']

FOR n=0,1 DO BEGIN

	sts = grd_read( filepath(root=root,'stdstar_c1_'+psf[n]+'.grd'), s1 )
	sts = grd_read( filepath(root=root,'stdstar_c2_'+psf[n]+'.grd'), s2 )
	sts = grd_read( filepath(root=root,'stdstar_c3_'+psf[n]+'.grd'), s3 )

	out = filepath(root=root,'standard_star_'+psf[n]+'.fts')

	mkhdr	, hdr, IsType(s1), size(s1,/dim), /extend
	fxaddpar, hdr, 'CAMERA'	,  1, ' Camera'
	fxaddpar, hdr, 'CX'	, ((size(s1,/dim))[0]+1)*0.5, ' Array index (base 1) for RA=0'
	fxaddpar, hdr, 'CY'	, ((size(s1,/dim))[1]+1)*0.5, ' Array index (base 1) for DEC=0'
	fxaddpar, hdr, 'BD'	, 40, ' Bins per degree'

	writefits, out, s1, hdr

	mkhdr	, hdr, IsType(s2), size(s2,/dim), /image
	fxaddpar, hdr, 'CAMERA'	,  2, ' Camera'
	fxaddpar, hdr, 'CX'	, ((size(s2,/dim))[0]+1)*0.5, ' Array index (base 1) for RA=0'
	fxaddpar, hdr, 'CY'	, ((size(s2,/dim))[1]+1)*0.5, ' Array index (base 1) for DEC=0'
	fxaddpar, hdr, 'BD'	, 40, ' Bins per degree'

	writefits, out, s2, hdr, /append

	mkhdr	, hdr, IsType(s3), size(s3,/dim), /image
	fxaddpar, hdr, 'CAMERA'	,  3, ' Camera'
	fxaddpar, hdr, 'CX'	, ((size(s3,/dim))[0]+1)*0.5, ' Array index (base 1) for RA=0'
	fxaddpar, hdr, 'CY'	, ((size(s3,/dim))[1]+1)*0.5, ' Array index (base 1) for DEC=0'
	fxaddpar, hdr, 'BD'	, 40, ' Bins per degree'

	writefits, out, s3, hdr, /append

	sts = gzip_file(out,/force)

ENDFOR

RETURN  &  END