FUNCTION read_head_from_drms, series, sunum

; EXAMPLE
; IDL> head = read_head_from_drms('mps_schunker.MT_256x256x61[2002.01.20_00:15:00_TAI]')
; IDL> head = read_head_from_drms('mps_schunker.MT_256x256x61')
; IDL> head = read_head_from_drms('mps_schunker.V2_mag1_512x512x1919[2010.04.03_22:53:15_TAI],'281474977109500')

;;;;;;;;;;;   SPAWNING USING DRMS   ;;;;;;;;;;;;;;;;;;;;;
if keyword_set(sunum) then begin
    spawn,'show_info ' + series + ' -P sunum=' + sunum + '| grep SUM', dir
    spawn,'show_info ' + series + ' -A sunum=' + sunum + '| grep fits', filename


; now header stuff
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=OBS_L0',OBS_L0
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=REF_L0',REF_L0
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=OBS_B0',OBS_B0
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=REF_B0',REF_B0
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=DAXIS1',DAXIS1
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=DAXIS2',DAXIS2
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=T_OBS_FI',T_OBS_FI
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=T_OBS_LA',T_OBS_LA
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=T_OBS_FI_step',T_OBS_FI_step
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=T_OBS_FI_epoch',T_OBS_FI_epoch
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=T_OBS_LA_step',T_OBS_LA_step
    spawn,'show_info ' + series + ' sunum=' + sunum + ' key=T_OBS_LA_epoch',T_OBS_LA_epoch
    




endif else begin
    spawn,'show_info ' + series + ' -P n=1 | grep SUM', dir
    spawn,'show_info ' + series + ' -A n=1 | grep fits', filename


; now header stuff
    spawn,'show_info ' + series + ' n=1 key=OBS_L0',OBS_L0
    spawn,'show_info ' + series + ' n=1 key=REF_L0',REF_L0
    spawn,'show_info ' + series + ' n=1 key=OBS_B0',OBS_B0
    spawn,'show_info ' + series + ' n=1 key=REF_B0',REF_B0
    spawn,'show_info ' + series + ' n=1 key=DAXIS1',DAXIS1
    spawn,'show_info ' + series + ' n=1 key=DAXIS2',DAXIS2
    spawn,'show_info ' + series + ' n=1 key=T_OBS_FI',T_OBS_FI
    spawn,'show_info ' + series + ' n=1 key=T_OBS_LA',T_OBS_LA
    spawn,'show_info ' + series + ' n=1 key=T_OBS_FI_step',T_OBS_FI_step
    spawn,'show_info ' + series + ' n=1 key=T_OBS_FI_epoch',T_OBS_FI_epoch
    spawn,'show_info ' + series + ' n=1 key=T_OBS_LA_step',T_OBS_LA_step
    spawn,'show_info ' + series + ' n=1 key=T_OBS_LA_epoch',T_OBS_LA_epoch
    
endelse


; read in fits file
head=headfits(dir+'/'+filename)




sxaddpar,head,'OBS_L0',float(OBS_L0(1)),' deg'
sxaddpar,head,'REF_L0',float(REF_L0(1)),' deg'
sxaddpar,head,'OBS_B0',float(OBS_B0(1)),' deg'
sxaddpar,head,'REF_B0',float(REF_B0(1)),' deg'
sxaddpar,head,'DAXIS1',float(DAXIS1(1)),' rad'
sxaddpar,head,'DAXIS2',float(DAXIS2(1)),' rad'
sxaddpar,head,'T_OBS_FI',T_OBS_FI(1)
sxaddpar,head,'T_OBS_LA',T_OBS_LA(1)
;sxaddpar,head,'T_OBS_FI_step',T_OBS_FI_step(1)
;sxaddpar,head,'T_OBS_FI_epoch',T_OBS_FI_epoch(1)
;sxaddpar,head,'T_OBS_LA_step',T_OBS_LA_step(1)
;sxaddpar,head,'T_OBS_LA_epoch',T_OBS_LA_epoch(1)
sxaddpar,head,'SUDIR',' ','  '+dir


return,  head


END