dimen=[512,512,40]
dir='/data/SOLMAG/DATA/IVM/NA-SURVEY/'
inim=dir+'/20031026/sodium/IMG0_20031026.1821_df_rf_ic_pl_sc_ds_sub_av_st_ca_c'
;inim=dir+'/20031026/sodium/IMG0_20031026.1832_df_rf_ic_pl_sc_ds_sub_av_st_ca_c' ; flare one
; uncomment below to read in the first time
OPENR, lun, inim[0], /GET_LUN
im = FLTARR(dimen[0],dimen[1],4,dimen[2])
READU,lun,im
CLOSE,lun
FREE_LUN,lun

head=ivm_rd_head(dir+'/20031026/IMGR_20031026.1821')
;nwave = head.n_image_sets - 1
;wvec = head.fp_set(0:nwave-1)     ; and their FP settings
nwave = head.n_image_sets
wvec = head.fp_set(0:nwave-1)     ; and their FP settings
; Compute the relative wavelengths in Angstroms
wvec = IVM_FP_STEP(5896.) * (wvec-wvec(0))
wvec=shift(wvec,1) ; make continuum where it should be
wvec=wvec-max(wvec)/2.
wvec=wvec/10. ;(put in nm)
xtt='589.6 + !M'+string(68B)+string(108B)+'!X (nm)'

;single-point example from same data as in PAPER_IVM_eg.pro
xtarget=246
;ytarget=268
ytarget=272
sp1=shift(reform(im[xtarget,ytarget,*,*]),0,1)
xtarget=246
ytarget=120
;ytarget=397
sp2=shift(reform(im[xtarget,ytarget,*,*]),0,1)

set_plot,'ps'
device,file='PAPER_IVM_spect_eg.eps',bits=8,/times,xsi=20,ysi=5
loadct,0
!p.multi=[0,4,1]
!p.charsize=1.5
!x.margin=[3.5,2]
!x.omargin=[4,0]
!y.margin=[4,1]
!y.omargin=[0,0]
plot,wvec,sp1[0,*],thick=4,xtit=xtt,yran=[0,3.2e4],/yst,xran=[-0.05,0.05],/xst
oplot,wvec,sp2[0,*],col=150,thick=6
plot,wvec,sp1[1,*],thick=4,xtit=xtt,yran=[-0.1,0.1],xran=[-0.05,0.05],/xst
oplot,wvec,sp2[1,*],col=150,thick=6
plot,wvec,sp1[2,*],thick=4,xtit=xtt,yran=[-0.1,0.1],xran=[-0.05,0.05],/xst
oplot,wvec,sp2[2,*],col=150,thick=6
plot,wvec,sp1[3,*],thick=4,xtit=xtt,yran=[-0.2,0.2],xran=[-0.05,0.05],/xst
oplot,wvec,sp2[3,*],col=150,thick=6

device,/close
set_plot,'x'
reset