!P.CHARSIZE = 1.5
  !P.CHARTHICK = 2.5
  !P.THICK = 5.5

fidTDiffs_CAMRT_Eclipse = ncdf_open('TDiffs_LatTrace_CAMRT_001.nc')
fidTDiffs_CAMRT_After = ncdf_open('TDiffs_LatTrace_After_CAMRT_001.nc')
fidTDiffs_PHOT_Eclipse = ncdf_open('TDiffs_LatTrace_PHOT_002.nc')
fidTDiffs_PHOT_After = ncdf_open('TDiffs_LatTrace_After_PHOT_002.nc')
fidTDiffs_CO2NIR_Eclipse = ncdf_open('TDiffs_LatTrace_CO2NIR_003.nc')
fidTDiffs_CO2NIR_After = ncdf_open('TDiffs_LatTrace_After_CO2NIR_003.nc')
fidTDiffs_EUVRATES_Eclipse = ncdf_open('TDiffs_LatTrace_EUVRATES_004.nc')
fidTDiffs_EUVRATES_After = ncdf_open('TDiffs_LatTrace_After_EUVRATES_004.nc')
fidTDiffs_QTHERMAL_Eclipse = ncdf_open('TDiffs_LatTrace_QTHERMAL_005.nc')
fidTDiffs_QTHERMAL_After = ncdf_open('TDiffs_LatTrace_After_QTHERMAL_005.nc')

ncdf_varget, fidTDiffs_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidTDiffs_CAMRT_Eclipse, 'data', TDiffs_CAMRT_Eclipse
ncdf_varget, fidTDiffs_CAMRT_After, 'data', TDiffs_CAMRT_After
ncdf_varget, fidTDiffs_PHOT_Eclipse, 'data', TDiffs_PHOT_Eclipse
ncdf_varget, fidTDiffs_PHOT_After, 'data', TDiffs_PHOT_After
ncdf_varget, fidTDiffs_CO2NIR_Eclipse, 'data', TDiffs_CO2NIR_Eclipse
ncdf_varget, fidTDiffs_CO2NIR_After, 'data', TDiffs_CO2NIR_After
ncdf_varget, fidTDiffs_EUVRATES_Eclipse, 'data', TDiffs_EUVRATES_Eclipse
ncdf_varget, fidTDiffs_EUVRATES_After, 'data', TDiffs_EUVRATES_After
ncdf_varget, fidTDiffs_QTHERMAL_Eclipse, 'data', TDiffs_QTHERMAL_Eclipse
ncdf_varget, fidTDiffs_QTHERMAL_After, 'data', TDiffs_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_TDiffs_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X Temperature Differences Eclipse-Baseline Longitude=267.50'
xTitle = 'Latitude (Degrees)'
yTitle = 'Temperature Difference (K)'
  
plot, lats,TDiffs_CAMRT_After,linestyle=2,yrange=[-40,10],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, TDiffs_CAMRT_Eclipse, thick=6
oplot, lats, TDiffs_PHOT_Eclipse, color=3, thick=6
oplot, lats, TDiffs_PHOT_After, linestyle=2, color=3, thick=6
oplot, lats, TDiffs_CO2NIR_Eclipse, color=6, thick=6
oplot, lats, TDiffs_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, TDiffs_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, TDiffs_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, TDiffs_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, TDiffs_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-85,-70],[-22,-22], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (18:15Z)', charthick=4
oplot, [-85,-70],[-24,-24], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', color=1, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, '+ PHOT ', color=3, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, '+ CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, '+ EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.20, '+ QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close

fidNmF2Diffs_CAMRT_Eclipse = ncdf_open('NmF2Diffs_LatTrace_CAMRT_001.nc')
fidNmF2Diffs_CAMRT_After = ncdf_open('NmF2Diffs_LatTrace_After_CAMRT_001.nc')
fidNmF2Diffs_PHOT_Eclipse = ncdf_open('NmF2Diffs_LatTrace_PHOT_002.nc')
fidNmF2Diffs_PHOT_After = ncdf_open('NmF2Diffs_LatTrace_After_PHOT_002.nc')
fidNmF2Diffs_CO2NIR_Eclipse = ncdf_open('NmF2Diffs_LatTrace_CO2NIR_003.nc')
fidNmF2Diffs_CO2NIR_After = ncdf_open('NmF2Diffs_LatTrace_After_CO2NIR_003.nc')
fidNmF2Diffs_EUVRATES_Eclipse = ncdf_open('NmF2Diffs_LatTrace_EUVRATES_004.nc')
fidNmF2Diffs_EUVRATES_After = ncdf_open('NmF2Diffs_LatTrace_After_EUVRATES_004.nc')
fidNmF2Diffs_QTHERMAL_Eclipse = ncdf_open('NmF2Diffs_LatTrace_QTHERMAL_005.nc')
fidNmF2Diffs_QTHERMAL_After = ncdf_open('NmF2Diffs_LatTrace_After_QTHERMAL_005.nc')

ncdf_varget, fidNmF2Diffs_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidNmF2Diffs_CAMRT_Eclipse, 'data', NmF2Diffs_CAMRT_Eclipse
ncdf_varget, fidNmF2Diffs_CAMRT_After, 'data', NmF2Diffs_CAMRT_After
ncdf_varget, fidNmF2Diffs_PHOT_Eclipse, 'data', NmF2Diffs_PHOT_Eclipse
ncdf_varget, fidNmF2Diffs_PHOT_After, 'data', NmF2Diffs_PHOT_After
ncdf_varget, fidNmF2Diffs_CO2NIR_Eclipse, 'data', NmF2Diffs_CO2NIR_Eclipse
ncdf_varget, fidNmF2Diffs_CO2NIR_After, 'data', NmF2Diffs_CO2NIR_After
ncdf_varget, fidNmF2Diffs_EUVRATES_Eclipse, 'data', NmF2Diffs_EUVRATES_Eclipse
ncdf_varget, fidNmF2Diffs_EUVRATES_After, 'data', NmF2Diffs_EUVRATES_After
ncdf_varget, fidNmF2Diffs_QTHERMAL_Eclipse, 'data', NmF2Diffs_QTHERMAL_Eclipse
ncdf_varget, fidNmF2Diffs_QTHERMAL_After, 'data', NmF2Diffs_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_NmF2Diffs_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X NmF2 Differences Eclipse-Baseline Longitude=267.50'
xTitle = 'Latitude (Degrees)'
yTitle = 'NmF2 Difference (cm-3)'

plot, lats,NmF2Diffs_CAMRT_After,linestyle=2,yrange=[-300000,100000],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, NmF2Diffs_CAMRT_Eclipse, thick=6
oplot, lats, NmF2Diffs_PHOT_Eclipse, color=3, thick=6
oplot, lats, NmF2Diffs_PHOT_After, linestyle=2, color=3, thick=6
oplot, lats, NmF2Diffs_CO2NIR_Eclipse, color=6, thick=6
oplot, lats, NmF2Diffs_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, NmF2Diffs_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, NmF2Diffs_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, NmF2Diffs_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, NmF2Diffs_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-85,-70],[-160000,-160000], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (18:15Z)', charthick=4
oplot, [-85,-70],[-172000,-172000], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', color=1, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, '+ PHOT ', color=3, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, '+ CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, '+ EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.20, '+ QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close

fidTDiffs_001_CAMRT_Eclipse = ncdf_open('TDiffs_LatTrace_CAMRT_001.nc')
fidTDiffs_001_CAMRT_After = ncdf_open('TDiffs_LatTrace_After_CAMRT_001.nc')
fidTDiffs_001_PHOT_Eclipse = ncdf_open('TDiffs_LatTrace_PHOT_001.nc')
fidTDiffs_001_PHOT_After = ncdf_open('TDiffs_LatTrace_After_PHOT_001.nc')
fidTDiffs_001_CO2NIR_Eclipse = ncdf_open('TDiffs_LatTrace_CO2NIR_001.nc')
fidTDiffs_001_CO2NIR_After = ncdf_open('TDiffs_LatTrace_After_CO2NIR_001.nc')
fidTDiffs_001_EUVRATES_Eclipse = ncdf_open('TDiffs_LatTrace_EUVRATES_001.nc')
fidTDiffs_001_EUVRATES_After = ncdf_open('TDiffs_LatTrace_After_EUVRATES_001.nc')
fidTDiffs_001_QTHERMAL_Eclipse = ncdf_open('TDiffs_LatTrace_QTHERMAL_001.nc')
fidTDiffs_001_QTHERMAL_After = ncdf_open('TDiffs_LatTrace_After_QTHERMAL_001.nc')

ncdf_varget, fidTDiffs_001_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidTDiffs_001_CAMRT_Eclipse, 'data', TDiffs_001_CAMRT_Eclipse
ncdf_varget, fidTDiffs_001_CAMRT_After, 'data', TDiffs_001_CAMRT_After
ncdf_varget, fidTDiffs_001_PHOT_Eclipse, 'data', TDiffs_001_PHOT_Eclipse
ncdf_varget, fidTDiffs_001_PHOT_After, 'data', TDiffs_001_PHOT_After
ncdf_varget, fidTDiffs_001_CO2NIR_Eclipse, 'data', TDiffs_001_CO2NIR_Eclipse
ncdf_varget, fidTDiffs_001_CO2NIR_After, 'data', TDiffs_001_CO2NIR_After
ncdf_varget, fidTDiffs_001_EUVRATES_Eclipse, 'data', TDiffs_001_EUVRATES_Eclipse
ncdf_varget, fidTDiffs_001_EUVRATES_After, 'data', TDiffs_001_EUVRATES_After
ncdf_varget, fidTDiffs_001_QTHERMAL_Eclipse, 'data', TDiffs_001_QTHERMAL_Eclipse
ncdf_varget, fidTDiffs_001_QTHERMAL_After, 'data', TDiffs_001_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_TDiffs_001_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X Temperature Differences Eclipse-Baseline Longitude=267.50'
xTitle = 'Latitude (Degrees)'
yTitle = 'Temperature Difference (K)'
  
plot, lats,TDiffs_001_CAMRT_After,linestyle=2,yrange=[-30,5],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, TDiffs_001_CAMRT_Eclipse, thick=6
oplot, lats, TDiffs_001_PHOT_Eclipse, color=3, thick=6
oplot, lats, TDiffs_001_PHOT_After, linestyle=2, color=3, thick=6
oplot, lats, TDiffs_001_CO2NIR_Eclipse, color=6, thick=6
oplot, lats, TDiffs_001_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, TDiffs_001_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, TDiffs_001_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, TDiffs_001_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, TDiffs_001_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-82,-69],[-17.5,-17.5], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (18:15Z)', charthick=4
oplot, [-82,-69],[-19,-19], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', color=1, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, 'PHOT ', color=3, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, 'CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, 'EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.20, 'QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close

fidNmF2Diffs_001_CAMRT_Eclipse = ncdf_open('NmF2Diffs_LatTrace_CAMRT_001.nc')
fidNmF2Diffs_001_CAMRT_After = ncdf_open('NmF2Diffs_LatTrace_After_CAMRT_001.nc')
fidNmF2Diffs_001_PHOT_Eclipse = ncdf_open('NmF2Diffs_LatTrace_PHOT_001.nc')
fidNmF2Diffs_001_PHOT_After = ncdf_open('NmF2Diffs_LatTrace_After_PHOT_001.nc')
fidNmF2Diffs_001_CO2NIR_Eclipse = ncdf_open('NmF2Diffs_LatTrace_CO2NIR_001.nc')
fidNmF2Diffs_001_CO2NIR_After = ncdf_open('NmF2Diffs_LatTrace_After_CO2NIR_001.nc')
fidNmF2Diffs_001_EUVRATES_Eclipse = ncdf_open('NmF2Diffs_LatTrace_EUVRATES_001.nc')
fidNmF2Diffs_001_EUVRATES_After = ncdf_open('NmF2Diffs_LatTrace_After_EUVRATES_001.nc')
fidNmF2Diffs_001_QTHERMAL_Eclipse = ncdf_open('NmF2Diffs_LatTrace_QTHERMAL_001.nc')
fidNmF2Diffs_001_QTHERMAL_After = ncdf_open('NmF2Diffs_LatTrace_After_QTHERMAL_001.nc')

ncdf_varget, fidNmF2Diffs_001_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidNmF2Diffs_001_CAMRT_Eclipse, 'data', NmF2Diffs_001_CAMRT_Eclipse
ncdf_varget, fidNmF2Diffs_001_CAMRT_After, 'data', NmF2Diffs_001_CAMRT_After
ncdf_varget, fidNmF2Diffs_001_PHOT_Eclipse, 'data', NmF2Diffs_001_PHOT_Eclipse
ncdf_varget, fidNmF2Diffs_001_PHOT_After, 'data', NmF2Diffs_001_PHOT_After
ncdf_varget, fidNmF2Diffs_001_CO2NIR_Eclipse, 'data', NmF2Diffs_001_CO2NIR_Eclipse
ncdf_varget, fidNmF2Diffs_001_CO2NIR_After, 'data', NmF2Diffs_001_CO2NIR_After
ncdf_varget, fidNmF2Diffs_001_EUVRATES_Eclipse, 'data', NmF2Diffs_001_EUVRATES_Eclipse
ncdf_varget, fidNmF2Diffs_001_EUVRATES_After, 'data', NmF2Diffs_001_EUVRATES_After
ncdf_varget, fidNmF2Diffs_001_QTHERMAL_Eclipse, 'data', NmF2Diffs_001_QTHERMAL_Eclipse
ncdf_varget, fidNmF2Diffs_001_QTHERMAL_After, 'data', NmF2Diffs_001_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_NmF2Diffs_001_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X NmF2 Differences Eclipse-Baseline Longitude=267.50'
Title = 'Latitude (Degrees)'
yTitle = 'NmF2 Difference (cm-3)'
  
plot, lats,NmF2Diffs_001_CAMRT_After,linestyle=2,yrange=[-300000,200000],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, NmF2Diffs_001_CAMRT_Eclipse, thick=6
oplot, lats, NmF2Diffs_001_PHOT_Eclipse, color=3, thick=6
oplot, lats, NmF2Diffs_001_PHOT_After, linestyle=2, color=3, thick=6
oplot, lats, NmF2Diffs_001_CO2NIR_Eclipse, color=6, thick=6
oplot, lats, NmF2Diffs_001_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, NmF2Diffs_001_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, NmF2Diffs_001_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, NmF2Diffs_001_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, NmF2Diffs_001_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-82,-69],[-120000,-120000], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (18:15Z)', charthick=4
oplot, [-82,-69],[-140000,-140000], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', color=1, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, 'PHOT ', color=3, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, 'CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, 'EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.20, 'QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close


fidTDiffsLon285_CAMRT_Eclipse = ncdf_open('TDiffs_LatTrace_CAMRT_001_Lon285.nc')
fidTDiffsLon285_CAMRT_After = ncdf_open('TDiffs_LatTrace_After_CAMRT_001_Lon285.nc')
fidTDiffsLon285_PHOT_Eclipse = ncdf_open('TDiffs_LatTrace_PHOT_002_Lon285.nc')
fidTDiffsLon285_PHOT_After = ncdf_open('TDiffs_LatTrace_After_PHOT_002_Lon285.nc')
fidTDiffsLon285_CO2NIR_Eclipse = ncdf_open('TDiffs_LatTrace_CO2NIR_003_Lon285.nc')
fidTDiffsLon285_CO2NIR_After = ncdf_open('TDiffs_LatTrace_After_CO2NIR_003_Lon285.nc')
fidTDiffsLon285_EUVRATES_Eclipse = ncdf_open('TDiffs_LatTrace_EUVRATES_004_Lon285.nc')
fidTDiffsLon285_EUVRATES_After = ncdf_open('TDiffs_LatTrace_After_EUVRATES_004_Lon285.nc')
fidTDiffsLon285_QTHERMAL_Eclipse = ncdf_open('TDiffs_LatTrace_QTHERMAL_005_Lon285.nc')
fidTDiffsLon285_QTHERMAL_After = ncdf_open('TDiffs_LatTrace_After_QTHERMAL_005_Lon285.nc')

ncdf_varget, fidTDiffsLon285_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidTDiffsLon285_CAMRT_Eclipse, 'data', TDiffsLon285_CAMRT_Eclipse
ncdf_varget, fidTDiffsLon285_CAMRT_After, 'data', TDiffsLon285_CAMRT_After
ncdf_varget, fidTDiffsLon285_PHOT_Eclipse, 'data', TDiffsLon285_PHOT_Eclipse
ncdf_varget, fidTDiffsLon285_PHOT_After, 'data', TDiffsLon285_PHOT_After
ncdf_varget, fidTDiffsLon285_CO2NIR_Eclipse, 'data', TDiffsLon285_CO2NIR_Eclipse
ncdf_varget, fidTDiffsLon285_CO2NIR_After, 'data', TDiffsLon285_CO2NIR_After
ncdf_varget, fidTDiffsLon285_EUVRATES_Eclipse, 'data', TDiffsLon285_EUVRATES_Eclipse
ncdf_varget, fidTDiffsLon285_EUVRATES_After, 'data', TDiffsLon285_EUVRATES_After
ncdf_varget, fidTDiffsLon285_QTHERMAL_Eclipse, 'data', TDiffsLon285_QTHERMAL_Eclipse
ncdf_varget, fidTDiffsLon285_QTHERMAL_After, 'data', TDiffsLon285_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_TDiffsLon285_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X Temperature Differences 250km Eclipse-Baseline Lon=285.0'
xTitle = 'Latitude (Degrees)'
yTitle = 'Temperature Difference (K)'
  
plot, lats,TDiffsLon285_CAMRT_After,linestyle=2,yrange=[-40,10],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, TDiffsLon285_CAMRT_Eclipse, thick=6
oplot, lats, TDiffsLon285_PHOT_Eclipse, color=3, thick=6
oplot, lats, TDiffsLon285_PHOT_After, linestyle=2, color=3, thick=6
;oplot, lats, TDiffsLon285_CO2NIR_Eclipse, color=6, thick=6
;oplot, lats, TDiffsLon285_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, TDiffsLon285_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, TDiffsLon285_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, TDiffsLon285_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, TDiffsLon285_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-82,-69],[-22,-22], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (19:05Z)', charthick=4
oplot, [-82,-69],[-24,-24], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, '+ PHOT + CO2NIR ', color=3, charsize=1.5, charthick=4
;xyouts, /NORMAL, 0.25, 0.26, '+ CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, '+ EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, '+ QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close

fidNmF2DiffsLon285_CAMRT_Eclipse = ncdf_open('NmF2Diffs_LatTrace_CAMRT_001_Lon285.nc')
fidNmF2DiffsLon285_CAMRT_After = ncdf_open('NmF2Diffs_LatTrace_After_CAMRT_001_Lon285.nc')
fidNmF2DiffsLon285_PHOT_Eclipse = ncdf_open('NmF2Diffs_LatTrace_PHOT_002_Lon285.nc')
fidNmF2DiffsLon285_PHOT_After = ncdf_open('NmF2Diffs_LatTrace_After_PHOT_002_Lon285.nc')
fidNmF2DiffsLon285_CO2NIR_Eclipse = ncdf_open('NmF2Diffs_LatTrace_CO2NIR_003_Lon285.nc')
fidNmF2DiffsLon285_CO2NIR_After = ncdf_open('NmF2Diffs_LatTrace_After_CO2NIR_003_Lon285.nc')
fidNmF2DiffsLon285_EUVRATES_Eclipse = ncdf_open('NmF2Diffs_LatTrace_EUVRATES_004_Lon285.nc')
fidNmF2DiffsLon285_EUVRATES_After = ncdf_open('NmF2Diffs_LatTrace_After_EUVRATES_004_Lon285.nc')
fidNmF2DiffsLon285_QTHERMAL_Eclipse = ncdf_open('NmF2Diffs_LatTrace_QTHERMAL_005_Lon285.nc')
fidNmF2DiffsLon285_QTHERMAL_After = ncdf_open('NmF2Diffs_LatTrace_After_QTHERMAL_005_Lon285.nc')

ncdf_varget, fidNmF2DiffsLon285_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidNmF2DiffsLon285_CAMRT_Eclipse, 'data', NmF2DiffsLon285_CAMRT_Eclipse
ncdf_varget, fidNmF2DiffsLon285_CAMRT_After, 'data', NmF2DiffsLon285_CAMRT_After
ncdf_varget, fidNmF2DiffsLon285_PHOT_Eclipse, 'data', NmF2DiffsLon285_PHOT_Eclipse
ncdf_varget, fidNmF2DiffsLon285_PHOT_After, 'data', NmF2DiffsLon285_PHOT_After
ncdf_varget, fidNmF2DiffsLon285_CO2NIR_Eclipse, 'data', NmF2DiffsLon285_CO2NIR_Eclipse
ncdf_varget, fidNmF2DiffsLon285_CO2NIR_After, 'data', NmF2DiffsLon285_CO2NIR_After
ncdf_varget, fidNmF2DiffsLon285_EUVRATES_Eclipse, 'data', NmF2DiffsLon285_EUVRATES_Eclipse
ncdf_varget, fidNmF2DiffsLon285_EUVRATES_After, 'data', NmF2DiffsLon285_EUVRATES_After
ncdf_varget, fidNmF2DiffsLon285_QTHERMAL_Eclipse, 'data', NmF2DiffsLon285_QTHERMAL_Eclipse
ncdf_varget, fidNmF2DiffsLon285_QTHERMAL_After, 'data', NmF2DiffsLon285_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_NmF2DiffsLon285_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X NmF2 Differences 250km Eclipse-Baseline Lon=285.0'
xTitle = 'Latitude (Degrees)'
yTitle = 'NmF2 Difference (cm-3)'
  
plot,lats,NmF2DiffsLon285_CAMRT_After,linestyle=2,yrange=[-300000,100000],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, NmF2DiffsLon285_CAMRT_Eclipse, thick=6
oplot, lats, NmF2DiffsLon285_PHOT_Eclipse, color=3, thick=6
oplot, lats, NmF2DiffsLon285_PHOT_After, linestyle=2, color=3, thick=6
;oplot, lats, NmF2DiffsLon285_CO2NIR_Eclipse, color=6, thick=6
;oplot, lats, NmF2DiffsLon285_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, NmF2DiffsLon285_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, NmF2DiffsLon285_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, NmF2DiffsLon285_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, NmF2DiffsLon285_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-82,-69],[-160000,-160000], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (19:05Z)', charthick=4
oplot, [-82,-69],[-175000,-175000], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, '+ PHOT + CO2NIR ', color=3, charsize=1.5, charthick=4
;xyouts, /NORMAL, 0.25, 0.26, '+ CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, '+ EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, '+ QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close

fidForcingsLon285_CAMRT_Eclipse = ncdf_open('Forcings_LatTrace_CAMRT_001_Lon285.nc')
fidForcingsLon285_CAMRT_After = ncdf_open('Forcings_LatTrace_After_CAMRT_001_Lon285.nc')
fidForcingsLon285_PHOT_Eclipse = ncdf_open('Forcings_LatTrace_PHOT_002_Lon285.nc')
fidForcingsLon285_PHOT_After = ncdf_open('Forcings_LatTrace_After_PHOT_002_Lon285.nc')
fidForcingsLon285_CO2NIR_Eclipse = ncdf_open('Forcings_LatTrace_CO2NIR_003_Lon285.nc')
fidForcingsLon285_CO2NIR_After = ncdf_open('Forcings_LatTrace_After_CO2NIR_003_Lon285.nc')
fidForcingsLon285_EUVRATES_Eclipse = ncdf_open('Forcings_LatTrace_EUVRATES_004_Lon285.nc')
fidForcingsLon285_EUVRATES_After = ncdf_open('Forcings_LatTrace_After_EUVRATES_004_Lon285.nc')
fidForcingsLon285_QTHERMAL_Eclipse = ncdf_open('Forcings_LatTrace_QTHERMAL_005_Lon285.nc')
fidForcingsLon285_QTHERMAL_After = ncdf_open('Forcings_LatTrace_After_QTHERMAL_005_Lon285.nc')

ncdf_varget, fidForcingsLon285_CAMRT_Eclipse, 'x', lats

ncdf_varget, fidForcingsLon285_CAMRT_Eclipse, 'data', ForcingsLon285_CAMRT_Eclipse
ncdf_varget, fidFrocingsLon285_CAMRT_After, 'data', ForcingsLon285_CAMRT_After
ncdf_varget, fidForcingsLon285_PHOT_Eclipse, 'data', ForcingsLon285_PHOT_Eclipse
ncdf_varget, fidForcingsLon285_PHOT_After, 'data', ForcingsLon285_PHOT_After
ncdf_varget, fidForcingsLon285_CO2NIR_Eclipse, 'data', ForcingsLon285_CO2NIR_Eclipse
ncdf_varget, fidForcingsLon285_CO2NIR_After, 'data', ForcingsLon285_CO2NIR_After
ncdf_varget, fidForcingsLon285_EUVRATES_Eclipse, 'data', ForcingsLon285_EUVRATES_Eclipse
ncdf_varget, fidForcingsLon285_EUVRATES_After, 'data', ForcingsLon285_EUVRATES_After
ncdf_varget, fidForcingsLon285_QTHERMAL_Eclipse, 'data', ForcingsLon285_QTHERMAL_Eclipse
ncdf_varget, fidForcingsLon285_QTHERMAL_After, 'data', ForcingsLon285_QTHERMAL_After

set_plot, 'ps'
psfile = 'eclipse_ForcingsLon285_LatTrace.ps'
device, /color, /land, bits_per_pixel=8,filename=psfile
LOADCT, 39, ncolors = 16

Title0 = ' WACCM-X Forcing Differences 250km Eclipse-Baseline Lon=285.0'
xTitle = 'Latitude (Degrees)'
yTitle = 'Forcing Difference (K)'
  
plot, lats,ForcingsLon285_CAMRT_After,linestyle=2,yrange=[-40,10],thick=6,xrange=[-90.0,90.0],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle
oplot, lats, ForcingsLon285_CAMRT_Eclipse, thick=6
oplot, lats, ForcingsLon285_PHOT_Eclipse, color=3, thick=6
oplot, lats, ForcingsLon285_PHOT_After, linestyle=2, color=3, thick=6
;oplot, lats, ForcingsLon285_CO2NIR_Eclipse, color=6, thick=6
;oplot, lats, ForcingsLon285_CO2NIR_After, linestyle=2, color=6, thick=6
oplot, lats, ForcingsLon285_EUVRATES_Eclipse, color=12, thick=6
oplot, lats, ForcingsLon285_EUVRATES_After, linestyle=2, color=12, thick=6
oplot, lats, ForcingsLon285_QTHERMAL_Eclipse, color=14, thick=6
oplot, lats, ForcingsLon285_QTHERMAL_After, linestyle=2, color=14, thick=6

oplot, [-82,-69],[-22,-22], thick=6
xyouts, /NORMAL, 0.25, 0.40, 'ECLIPSE (19:05Z)', charthick=4
oplot, [-82,-69],[-24,-24], linestyle=2, thick=6
xyouts, /NORMAL, 0.25, 0.37, 'AFTER ECLIPSE (24:00Z)', charthick=4
xyouts, /NORMAL, 0.25, 0.32, 'CAMRT ', charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.29, '+ PHOT + CO2NIR ', color=3, charsize=1.5, charthick=4
;xyouts, /NORMAL, 0.25, 0.26, '+ CO2NIR ', color=6, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.26, '+ EUVRATES ', color=12, charsize=1.5, charthick=4
xyouts, /NORMAL, 0.25, 0.23, '+ QTHERMAL ', color=14, charsize=1.5, charthick=4
device, /close