pro main

common alldata,emap4,emap304,aug21_171,aug21_195,aug21_284,a_he1,a_si12,a_mg10,a_o5,a_fe19,sxt_qmap,mdi_map,bbso_map,map_pics1,map_pics2,map_chip

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

spawn,'gunzip alldata_21.sav.gz'
restore,'alldata_21.sav'
spawn,'gzip alldata_21.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
;

set_line_color
!p.multi=0
!p.charsize=.8
!p.ticklen=-.02

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_21.ps',/color

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

parts1 = str_sep(aug21_195.time,' ')
parts2 = str_sep(emap4.time,' ')
usetitle='Aug 21 TRACE:195 '+parts1(1)+' EIT:195 '+parts2(1)
plot_map,aug21_195,/square,xrange=[100, 500], yrange=[100, 500],title=usetitle
plot_map,emap4,/over,/square, levels=[1000,2000, 3000, 4000],lcolor=3,trans=[-1,4.],/drotate
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
;  **TRACE shift = [1,-4] w.r.t. EIT**
;  EIT shift w.r.t TRACE = [-1,4]
;

;
; 2) TRACE 171 vs CDS Mg10
;

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

levmg10 = [0.,.8]

parts1 = str_sep(aug21_171.time,' ')
parts2 = str_sep(a_mg10.time,' ')
usetitle='Aug 21 TRACE:171 '+parts1(1)+' CDS:Mg10 '+parts2(1)
plot_map,aug21_171,/square, xrange=[100, 500], yrange=[100, 500],title=usetitle
plot_map,a_mg10,/over,/square,/drotate,levels=levmg10,lcolor=4,trans=[-3,14]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; CDS w.r.t. TRACE = [-3,14]
;

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

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

parts1 = str_sep(emap4.time,' ')
parts2 = str_sep(a_mg10.time,' ')
usetitle='Aug 21 EIT:195 '+parts1(1)+' CDS:Mg10 '+parts2(1)
plot_map,emap4,/square, xrange=[100, 500], yrange=[100, 500],title=usetitle
plot_map,a_mg10,/over,/square,/drotate,levels=levmg10,lcolor=4,trans=[-2,10]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 4) Now SXT vs CDS Fe 19
;

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

levfe19 = [0.,.79]

parts1 = str_sep(sxt_qmap(3).time,' ')
parts2 = str_sep(a_fe19.time,' ')
usetitle='Aug 21 SXT '+parts1(1)+' CDS:Fe19 '+parts2(1)
plot_map, sxt_qmap(3),/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,a_fe19,/over,/square,/drotate,levels=levfe19,trans=[15,25],lcolor=3
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; Note: this one is difficult to do, because the resolution is so
; bad for both - but at best guess CDS w.r.t. SXT is [15,25]
;

;
; 5) as another test, plot SXT w.r.t. TRACE 284 - should be [-18,-11]
; - show also with 295 and CDS Mg10 and Fe19 - because
; this illustrates how TRACE 284 and SXT are really rather
; different
;

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

parts1 = str_sep(aug21_284.time,' ')
parts2 = str_sep(sxt_qmap(3).time,' ')
parts3 = str_sep(a_mg10.time,' ')
usetitle='Aug 21 TRACE:284 '+parts1(1)+' SXT '+parts2(1)+' CDS:Mg10,Fe19 '+parts3(1)
usextitle='TRACE:284 base, SXT yellow contours, Mg10 green, Fe19 red"
plot_map,aug21_284,/square, xrange=[100, 500], yrange=[100, 500],title=usetitle,xtitle=usextitle
plot_map, sxt_qmap(3), /over, /square,levels=[200, 300, 500, 1000, 1900],/drotate,trans=[-18,-11],lcolor=2
plot_map,a_mg10,/over,/square,/drotate,levels=levmg10,lcolor=4,trans=[-3,14]
plot_map,a_fe19,/over,/square, /drotate,levels=levfe19,trans=[-3,14],lcolor=3
if ntype eq 0 then stop
if ntype eq 1 then device,/close

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

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

parts1 = str_sep(emap4.time,' ')
parts2 = str_sep(sxt_qmap(3).time,' ')
parts3 = str_sep(a_mg10.time,' ')
usetitle='Aug 21 EIT:195 '+parts1(1)+' SXT '+parts2(1)+' CDS:Mg10,Fe19 '+parts3(1)
usextitle='EIT:195 base, SXT yellow contours, Mg10 green, Fe19 red"
plot_map,emap4,/square,  xrange=[100, 500], yrange=[100, 500],title=usetitle,xtitle=usextitle
plot_map, sxt_qmap(3), /over,/square,  levels=[200, 300, 500, 1000, 1900],/drotate,trans=[-17,-15],lcolor=2
plot_map,a_mg10,/over,/square, /drotate,levels=levmg10,lcolor=4,trans=[-2,10]
plot_map,a_fe19,/over,/square, /drotate,levels=levfe19,lcolor=3,trans=[-2,10]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 7) now get MLSO CHIP 10830 lined up with CDS He1
;

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

levhe1 = [0.,.73]

parts1 = str_sep(map_chip.time,' ')
parts2 = str_sep(a_he1.time,' ')
usetitle='Aug 21 CHIP:He10830 '+parts1(1)+' CDS:He1 '+parts2(1)
plot_map,map_chip,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,a_he1,/over,/square,/drotate,levels=levhe1,lcolor=3,trans=[4,17]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; so CHIP w.r.t. CDS is [-4,-17]
;

;
; ***this should mean that CHIP w.r.t. EIT is [-6,-7] ***
; 8) plot this:
;

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

parts1 = str_sep(emap4.time,' ')
parts2 = str_sep(map_chip.time,' ')
parts3 = str_sep(a_he1.time,' ')
usetitle='Aug 21 EIT:195 '+parts1(1)+' CHIP:He10830 '+parts2(1)+' CDS:He1 '+parts3(1)
usextitle='EIT:195 base, CHIP:10830 yellow contour, CDS:He1 red'
plot_map,emap4,/square, xrange=[100, 500], yrange=[100, 500],title=usetitle,xtitle=usextitle
plot_map,map_chip,/over,/square,/drotate,trans=[-6,-7],levels=[-500,-1000,-1500],lcolor=2
plot_map,a_he1,/over,/square,/drotate,levels=levhe1,lcolor=3,trans=[-2,10]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 9) Line up PICS H-alpha 19:18 image with EIT 304
;

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

parts1 = str_sep(emap304.time,' ')
parts2 = str_sep(map_pics2.time,' ')
usetitle='Aug 21 EIT:304 '+parts1(1)+' PICS:H-alpha '+parts2(1)
plot_map,emap304,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,map_pics2,/over,/square,/drotate,lcolor=3,trans=[9,-3]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; ***so PICS w.r.t. EIT should be [9,-3] ***
;

; 
; 10) now line up BBSO with PICS 18:22 image
;

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

parts1 = str_sep(bbso_map.time,' ')
parts2 = str_sep(map_pics1.time,' ')
usetitle='Aug 21 BBSO:H-alpha '+parts1(1)+' PICS:H-alpha '+parts2(1)
plot_map,bbso_map,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,map_pics1,/over,/square,/drotate,lcolor=3,trans=[2,2]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; so PICS w.r.t. BBSO is [2,2]
; and BBSO w.r.t. EIT should be [7,-5]
;

;
; 11) ASSUME MDI and EIT are aligned
; And line up BBSO and MDI using [7,-5]
;

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

parts1 = str_sep(mdi_map.time,' ')
parts2 = str_sep(bbso_map.time,' ')
usetitle='Aug 21 MDI '+parts1(1)+' BBSO:H-alpha '+parts2(1)
plot_map,mdi_map,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,bbso_map,/square,/over,/drotate,lcolor=3,trans=[7,-5],levels=[0,100,200]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 12) Also show MDI on BBSO
;

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

parts1 = str_sep(bbso_map.time,' ')
parts2 = str_sep(mdi_map.time,' ')
usetitle='Aug 21 BBSO:H-alpha '+parts1(1)+' MDI '+parts2(1)
plot_map,bbso_map,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,mdi_map,/over,/square,/drotate,lcolor=4,trans=[-7,5],nlevels=3
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 13) Also show PICS on MDI using [9,-3]
;

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

parts1 = str_sep(mdi_map.time,' ')
parts2 = str_sep(map_pics1.time,' ')
usetitle='Aug 21 MDI '+parts1(1)+' PICS:H-alpha '+parts2(1)
plot_map,mdi_map,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,map_pics1,/over,/square,/drotate,lcolor=3,trans=[9,-3]
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 14) Also show MDI on PICS
;

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

parts1 = str_sep(map_pics1.time,' ')
parts2 = str_sep(mdi_map.time,' ')
usetitle='Aug 21 PICS:H-alpha '+parts1(1)+' MDI '+parts2(1)
plot_map,map_pics1,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,mdi_map,/over,/square,/drotate,lcolor=4,trans=[-9,3],nlevels=3
if ntype eq 0 then stop
if ntype eq 1 then device,/close

;
; 15) finally check how PICS lines up with CHIP
; should be [15,4]
;

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

parts1 = str_sep(map_pics1.time,' ')
parts2 = str_sep(map_chip.time,' ')
usetitle='Aug 21 PICS:H-alpha '+parts1(1)+' CHIP:He 10830 '+parts2(1)
plot_map,map_pics1,/square,xrange=[100,500],yrange=[100,500],title=usetitle
plot_map,map_chip,/over,/square,/drotate,trans=[-15,-4],levels=[-500,-1000,-1500],lcolor=2
if ntype eq 1 then device,/close
if ntype eq 0 then stop

end