pro main

common alldata,emap4,emap304,aug19_171,aug19_195,aug19_284,a_he1,a_si12,a_mg10,a_o5,sxt_pmap,mdi_map1,mdi_map2,bbso_map,map_pics,map_chip

;
; You will need the idl save set with all the maps of the data
;

spawn,'gunzip alldata_19.sav.gz'
restore,'alldata_19.sav'
spawn,'gzip alldata_19.sav'

read,'hard or soft copies? (1 hard, 0 soft)', ntype

;
; NOTE: if you say hard copy, it will make a bunch of
; postscript copies (same as on web page).  If you say
; soft, it will pull up xwindows to show plots, and stop
; between each one - you will have a choice of writing a
; gif file from the xwindow
;

!p.multi=0
!p.charsize=.9
!p.ticklen=-.01

if ntype eq 1 then set_plot,'ps'
if ntype eq 0 then set_plot,'x'

if ntype eq 0 then window,xs=800,ys=800,retain=3
if ntype eq 1 then device,filename='PLOTS/trace_eit_19.ps',/color

;
; DO ALIGNMENTS AGAINST EIT
;
;
; 1) TRACE 195 vs. EIT 195
;

loadct,10
gamma_ct,0.9

parts1 = str_sep(emap4.time,' ')
parts2 = str_sep(aug19_195.time,' ')
usetitle='Aug 19 TRACE:195 '+parts1(1)+' EIT:195 '+parts2(1)
plot_map,aug19_195,/square,xrange=[-350,50], yrange=[100, 500],title=usetitle
plot_map,emap4,/over,/square, levels=[0,1000],color=3,trans=[0,7.],/drotate
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/trace_eit_19.gif',tvrd()
if ntype eq 1 then device,/close

;
;  EIT shift w.r.t TRACE = [0,7]
;  **TRACE shift = [0,-7] w.r.t. EIT**
;

;
; 2) TRACE 171 vs CDS Mg10
;

loadct,2

if ntype eq 1 then device,filename='PLOTS/trace_cds_19.ps',/color

parts1 = str_sep(aug19_171.time,' ')
parts2 = str_sep(a_mg10.time,' ')
usetitle='Aug 19 TRACE:171 '+parts1(1)+' CDS:Mg10 '+parts2(1)
plot_map,aug19_171,/square, xrange=[-350,50], yrange=[100, 500],title=usetitle
plot_map,a_mg10,/over,/square,/drotate,color=3,trans=[-4,11]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/trace_cds_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; CDS w.r.t. TRACE = [-4,11]
;

;
; 3) so plot **CDS w.r.t. EIT is [-4,4]**
;

if ntype eq 1 then device,filename='PLOTS/eit_cds_19.ps',/color

loadct,0
gamma_ct,0.4

parts1 = str_sep(emap4.time,' ')
parts2 = str_sep(a_mg10.time,' ')
usetitle='Aug 19 EIT:195 '+parts1(1)+' CDS:Mg10 '+parts2(1)
plot_map,emap4,/square, xrange=[-350,50], yrange=[100, 500],title=usetitle
plot_map,a_mg10,/over,/square,/drotate,color=4,trans=[-4,4]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/eit_cds_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; 4) Now CDS Si 12 vs SXT
;

if ntype eq 1 then device,filename='PLOTS/cds_sxt_19.ps',/color

loadct,1
set_line_color

parts1 = str_sep(a_si12.time,' ')
parts2 = str_sep(sxt_pmap(3).time,' ')
usetitle='Aug 19 CDS:Si12 '+parts1(1)+' SXT:AlMg '+parts2(1)
plot_map, a_si12,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,sxt_pmap(3),/over,/square,/drotate,trans=[5,-8],lcolor=2
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/cds_sxt_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; SXT w.r.t. CDS is [5,-8]
;

;
; ***so the upshot is SXT w.r.t. EIT is [1,-4]***
; 5) plot this with EIT 195 and CDS lines
;

loadct,0
gamma_ct,0.4

if ntype eq 1 then device,filename='PLOTS/eit_sxt_19.ps',/color

parts1 = str_sep(emap4.time,' ')
parts2 = str_sep(sxt_pmap(3).time,' ')
usetitle='Aug 19 EIT:195 '+parts1(1)+' SXT '+parts2(1)
plot_map,emap4,/square,xrange=[-350,50], yrange=[100, 500],title=usetitle
plot_map, sxt_pmap(3), /over,/square,/drotate,trans=[1,-4],color=4
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/eit_sxt_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; 6) Line up BBSO H-alpha image with EIT 304
;

loadct,0

if ntype eq 1 then device,filename='PLOTS/bbso_eit_19.ps',/color

parts1 = str_sep(bbso_map.time,' ')
parts2 = str_sep(emap304.time,' ')
usetitle='Aug 19 BBSO:H-alpha '+parts1(1)+' EIT:304 '+parts2(1)
plot_map,bbso_map,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,emap4,/over,/square,/drotate,color=3,trans=[1,-10]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/bbso_eit_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; ***so BBSO w.r.t. EIT should be [-1,10] ***
;

; 
; 7) now line up BBSO with PICS image
;

if ntype eq 1 then device,filename='PLOTS/bbso_pics_19.ps',/color

parts1 = str_sep(bbso_map.time,' ')
parts2 = str_sep(map_pics.time,' ')
usetitle='Aug 19 BBSO:H-alpha '+parts1(1)+' PICS:H-alpha '+parts2(1)
plot_map,bbso_map,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,map_pics,/over,/square,/drotate,lcolor=3,trans=[5,-16]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/bbso_pics_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; so PICS w.r.t. BBSO is [5,-16]
; and PICS w.r.t. EIT should be [4,-6]
;


;
; 8) now get MLSO CHIP 10830 lined up with MLSO PICS H-alpha
; not exactly same temperature, but same time
;

gamma_ct,0.8 
set_line_color

if ntype eq 1 then device,filename='PLOTS/chip_pics_19.ps',/color

parts1 = str_sep(map_chip.time,' ')
parts2 = str_sep(map_pics.time,' ')
usetitle='Aug 19 CHIP:He10830 '+parts1(1)+' PICS:H-alpha '+parts2(1)
plot_map,map_chip,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,map_pics,/over,/square,/drotate,lcolor=3,trans=[5,-5]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/chip_pics_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; so CHIP w.r.t. PICS is [-5,5]
;

;
; ***this should mean that CHIP w.r.t. EIT is [-1,-1] ***
; 9) plot this
;

if ntype eq 1 then device,filename='PLOTS/chip_eit_19.ps',/color

parts1 = str_sep(map_chip.time,' ')
parts2 = str_sep(emap4.time,' ')
usetitle='Aug 19 CHIP:10830 '+parts1(1)+' EIT:195 '+parts2(1)
plot_map,map_chip,/square, xrange=[-350, 50], yrange=[100, 500],title=usetitle
plot_map,emap4,/over,/square,/drotate,trans=[1,1],lcolor=3
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/chip_eit_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; 10) also compare MLSO CHIP 10830 with CDS He1
; similar temperature, but different times
; 

if ntype eq 1 then device,filename='PLOTS/chip_cds_19.ps',/color

parts1 = str_sep(map_chip.time,' ')
parts2 = str_sep(a_he1.time,' ')
usetitle='Aug 19 CHIP:He10830 '+parts1(1)+' CDS:He1 '+parts2(1)
plot_map,map_chip,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,a_he1,/over,/square,/drotate,levels=levhe1,lcolor=3,trans=[-3,5]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/chip_cds_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; so CHIP w.r.t. CDS is [3,-5]
;

;
; 11) ASSUME MDI and EIT are aligned
; And line up MDI and BBSO 14:29 using [-1,10]
;

loadct,0
gamma_ct,0.8 

if ntype eq 1 then device,filename='PLOTS/mdi_bbso_19.ps',/color

parts1 = str_sep(mdi_map2.time,' ')
parts2 = str_sep(bbso_map.time,' ')
usetitle='Aug 19 MDI '+parts1(1)+' BBSO:H-alpha '+parts2(1)
plot_map,mdi_map2,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,bbso_map,/square,/over,/drotate,lcolor=3,trans=[-1,10]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/mdi_bbso_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; 12) Also show MDI on BBSO
;

if ntype eq 1 then device,filename='PLOTS/bbso_mdi_19.ps',/color

parts1 = str_sep(bbso_map.time,' ')
parts2 = str_sep(mdi_map2.time,' ')
usetitle='Aug 19 BBSO:H-alpha '+parts1(1)+' MDI '+parts2(1)
plot_map,bbso_map,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,mdi_map2,/over,/square,/drotate,lcolor=4,trans=[1,-10],nlevels=3
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/bbso_mdi_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; 13) Also show PICS on MDI using [4,-6]
;

if ntype eq 1 then device,filename='PLOTS/mdi_pics_19.ps',/color

parts1 = str_sep(mdi_map2.time,' ')
parts2 = str_sep(map_pics.time,' ')
usetitle='Aug 19 MDI '+parts1(1)+' PICS:H-alpha '+parts2(1)
plot_map,mdi_map2,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,map_pics,/over,/square,/drotate,lcolor=3,trans=[4,-6]
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/mdi_pics_19.gif',tvrd()
if ntype eq 1 then device,/close

;
; 14) Also show MDI on PICS
;

if ntype eq 1 then device,filename='PLOTS/pics_mdi_19.ps',/color

parts1 = str_sep(map_pics.time,' ')
parts2 = str_sep(mdi_map2.time,' ')
usetitle='Aug 19 PICS:H-alpha '+parts1(1)+' MDI '+parts2(1)
plot_map,map_pics,/square,xrange=[-350,50],yrange=[100,500],title=usetitle
plot_map,mdi_map2,/over,/square,/drotate,lcolor=4,trans=[-4,6],nlevels=3
if ntype eq 0 then read,'want a gif? (0 no, 1 yes)',gifres
if gifres eq 1 then write_gif,'PLOTS/pics_mdi_19.gif',tvrd()
if ntype eq 1 then device,/close

end