;;eis_prep;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;
;filename='eis_l0_20061206_191520.fits.gz'
;eis_prep,filename,/default,/save

;;getwindata;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;
filename2='eis_l1_20061206_191520.fits'
wd=eis_getwindata(filename2,195.12)
help,wd,/str
loadct,3
plot_image,total(wd.int[*,*,*],1)

;;auto_fit;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;
;eis_auto_fit,wd,fit,refwvl=195.12
restore,'fit.sav'
help,fit,/str
eis_fit_viewer, wd, fit
outfit=eis_tilt_correction(wd,fit)
outfit2=eis_orbit_spline(outfit,/setzero,/plotvel,spline=spline)
eis_fit_viewer, wd, outfit2
window,0,xs=512,ys=512*1.5  
!P.multi=[0,1,3]
loadct,3
plot_image,outfit2.int
load_vel
plot_image,outfit2.vel,max=-50,min=50
loadct,33
plot_image,outfit2.wid,max=0.05,min=0.0223
loadct,0
!P.multi=0

;;make map;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;
obj=obj_new('eis_data',filename2)
index_long=fitshead2struct((obj->gethdr())->getbte_hdr())
index_short=index_long
index_short.xcen=index_short.xcen-2.
index_short.ycen=index_short.ycen-16.
index2map,index_short,outfit2.int,map_intensity
index2map,index_short,outfit2.vel,map_velocity
index2map,index_short,outfit2.wid,map_width
map_velocity.data=-1.*map_velocity.data
!P.multi=[0,1,3]
loadct,3
plot_map,map_intensity,charsize=2
load_vel
plot_map,map_velocity,dmax=50.,dmin=-50.,charsize=2
loadct,33
plot_map,map_width,dmax=0.05,dmin=0.0223,charsize=2
loadct,0
!P.multi=0

;;spectral line plot;;;;;;;;;;;;;;;;;;;;;
;;;;
i=200
j=200
plot,wd.wvl[*]-outfit2.offset[i,j],wd.int[*,i,j],psym=10,xs=1
oplot,wd.wvl[*]-outfit2.offset[i,j],outfit2.int[i,j]/sqrt(2.*3.14)/$
outfit2.wid[i,j]*exp(-(wd.wvl[*]-outfit2.offset[i,j]-outfit2.cen[i,j])^2.$
/outfit2.wid[i,j]^2/2.)+outfit2.bg1[i,j]

;;;;practice;;;;;;;;;;;;;;;
filename='eis_l0_20090205_181427.fits.gz'
eis_prep,filename,/default,/save
filename2='eis_l1_20090205_181427.fits' 
wd=eis_getwindata(filename2,195.12)
loadct,3
plot_image,total(wd.int[*,*,*],1)
eis_auto_fit,wd,fit,refwvl=195.12
eis_fit_viewer, wd, fit
outfit=eis_tilt_correction(wd,fit)
outfit2=eis_orbit_spline(outfit,/setzero,/plotvel,spline=spline)
eis_fit_viewer, wd, outfit2
window,0,xs=512,ys=512*1.5
!P.multi=[0,3,1]
loadct,3
plot_image,outfit2.int
load_vel
plot_image,outfit2.vel-average(outfit2.vel),max=-20,min=20
loadct,33
plot_image,outfit2.wid,max=0.0446,min=0.0223
loadct,0
!P.multi=0

obj=obj_new('eis_data',filename2)
index_long=fitshead2struct((obj->gethdr())->getbte_hdr())
index_short=index_long
index_short.xcen=index_short.xcen-2.
index_short.ycen=index_short.ycen-16.
index2map,index_short,outfit2.int,map_intensity
index2map,index_short,outfit2.vel,map_velocity
index2map,index_short,outfit2.wid,map_width
map_velocity.data=-1.*map_velocity.data+average(map_velocity.data)
!P.multi=[0,3,1]
loadct,3
plot_map,map_intensity,charsize=2
load_vel
plot_map,map_velocity,dmax=20.,dmin=-20.,charsize=2
loadct,33
plot_map,map_width,dmax=0.0446,dmin=0.0223,charsize=2
loadct,0
!P.multi=0

;;;;;alignment;;;;;;
;;;;;;
;filename='eis_l0_20090205_211049.fits.gz'
;eis_prep,filename,/default,/save
filename2='eis_l1_20090205_211049.fits'
wd=eis_getwindata(filename2,195.12)
obj=obj_new('eis_data',filename2)
index_long=fitshead2struct((obj->gethdr())->getbte_hdr())
index_short=index_long
index_short.xcen=index_short.xcen-2.
index_short.ycen=index_short.ycen-16.
index2map,index_short,total(wd.int,1),map_intensity

read_eit,'/Users/imada/Work/HSC/EIT/20090205/eit_l1_20090205_212409.fits',data,index
index2map,data,index,map_eit

window,0,xs=512,ys=512
loadct,8
plot_map,map_eit,dmax=50,dmin=1
window,1,xs=256,ys=512
!P.multi=[0,2,1]
plot_map,map_eit,dmax=15,dmin=1,fov=map_intensity
loadct,3
plot_map,map_intensity
!P.multi=0

;;;;;;;
;filename='eis_l0_20090205_062423.fits.gz'
;eis_prep,filename,/default,/save
filename2='eis_l1_20090205_062423.fits'
wd2=eis_getwindata(filename2,195.12)
obj=obj_new('eis_data',filename2)
index_long=fitshead2struct((obj->gethdr())->getbte_hdr())
index_short=index_long
index_short.xcen=index_short.xcen-2.
index_short.ycen=index_short.ycen-16.
index2map,index_short,total(wd2.int,1),map_intensity2

read_eit,'/Users/imada/Work/HSC/EIT/20090205/eit_l1_20090205_091209.fits',data,index
index2map,data,index,map_eit2

window,0,xs=512,ys=512
loadct,8
plot_map,map_eit2,dmax=50,dmin=1
window,1,xs=1024,ys=512
!P.multi=[0,2,1]
plot_map,map_eit2,dmax=20,dmin=1,fov=map_intensity2
loadct,3
plot_map,map_intensity2
!P.multi=0

map_intensity2.xc=map_intensity2.xc+41.
map_intensity2.yc=map_intensity2.yc-49.

window,0,xs=512,ys=512
loadct,8
plot_map,map_eit2,dmax=50,dmin=1
window,1,xs=1024,ys=512
!P.multi=[0,2,1]
plot_map,map_eit2,dmax=20,dmin=5,fov=map_intensity2
plot_map,map_intensity2,/cnt,/over,levels=[5000]
loadct,3
plot_map,map_intensity2
plot_map,map_intensity2,/cnt,/over,levels=[5000]
!P.multi=0

;;;;;;;;;;;


map_intensity.xc=map_intensity.xc+41.
map_intensity.yc=map_intensity.yc-49.

window,0,xs=512,ys=512
loadct,8
plot_map,map_eit,dmax=50,dmin=1
window,1,xs=256,ys=512
!P.multi=[0,2,1]
plot_map,map_eit,dmax=15,dmin=3,fov=map_intensity
plot_map,map_intensity,/cnt,/over,levels=[5000]
loadct,3
plot_map,map_intensity
plot_map,map_intensity,/cnt,/over,levels=[5000]
!P.multi=0