PRO plot_waccm_wam_recon_lon_eq,zgm,tn,tnWAM,lon,lonWAM,twv1ZAllT,twv1ZAllTWAM ; This procedure was created on Tuesday, June 9th, 2009. ; This procedure was created by SOARS Protege McArthur (Mack) Jones Jr. ; This procedure will produce a contour w/ lat. on the x-axis and altitude on the y-axis ;newcolor !P.CHARSIZE = 1.5 !P.CHARTHICK = 2.5 !P.THICK = 5.5 LTimes = ['18','20','22','0','2','4','6'] twvtotal1to3 = REFORM(twv1ZAllT[*,*,*,1]) twvtotal1to6 = REFORM(twv1ZAllT[*,*,*,1]) twvtotal1to9 = REFORM(twv1ZAllT[*,*,*,1]) FOR iwvn=2,3 DO twvtotal1to3[*,*,*] = twvtotal1to3[*,*,*] + twv1ZAllT[*,*,*,iwvn] FOR iwvn=2,6 DO twvtotal1to6[*,*,*] = twvtotal1to6[*,*,*] + twv1ZAllT[*,*,*,iwvn] FOR iwvn=2,9 DO twvtotal1to9[*,*,*] = twvtotal1to9[*,*,*] + twv1ZAllT[*,*,*,iwvn] twvtotal1to3_shift = SHIFT(twvtotal1to3,72,0,0) twvtotal1to6_shift = SHIFT(twvtotal1to6,72,0,0) twvtotal1to9_shift = SHIFT(twvtotal1to9,72,0,0) twvtotal1to3WAM = REFORM(twv1ZAllTWAM[*,*,*,1]) twvtotal1to6WAM = REFORM(twv1ZAllTWAM[*,*,*,1]) twvtotal1to9WAM = REFORM(twv1ZAllTWAM[*,*,*,1]) FOR iwvn=2,3 DO twvtotal1to3WAM[*,*,*] = twvtotal1to3WAM[*,*,*] + twv1ZAllTWAM[*,*,*,iwvn] FOR iwvn=2,6 DO twvtotal1to6WAM[*,*,*] = twvtotal1to6WAM[*,*,*] + twv1ZAllTWAM[*,*,*,iwvn] FOR iwvn=2,9 DO twvtotal1to9WAM[*,*,*] = twvtotal1to9WAM[*,*,*] + twv1ZAllTWAM[*,*,*,iwvn] twvtotal1to3WAM_shift = SHIFT(twvtotal1to3WAM,96,0,0) twvtotal1to6WAM_shift = SHIFT(twvtotal1to6WAM,96,0,0) twvtotal1to9WAM_shift = SHIFT(twvtotal1to9WAM,96,0,0) tn_shift = SHIFT(tn,72,0,0,0) tnWAM_shift = SHIFT(tnWAM,96,0,0,0) lon_shift = lon - 180.0 lonWAM_shift = lonWAM - 180.0 ; nTAllW = dimsize[3] iLt = 47 ;Equator dimszgm = size(zgm) nx = dimszgm(1) ny = dimszgm(2) nT = dimszgm(4) tNZSub = FLTARR(nX,nT) tNZSub_shift = FLTARR(nX,nT) twvtotal1to3Sub = FLTARR(nX,nT) twvtotal1to3Sub_shift = FLTARR(nX,nT) twvtotal1to6Sub = FLTARR(nX,nT) twvtotal1to6Sub_shift = FLTARR(nX,nT) twvtotal1to9Sub = FLTARR(nX,nT) twvtotal1to9Sub_shift = FLTARR(nX,nT) ; ; Find altitude index for this time ; FOR iLn=0,nX-1 DO BEGIN FOR iT = 0, nT-1 DO BEGIN zWLnLtT = zgm[iLn,iLt,*,iT] zInd = WHERE(zWLnLtT GE 275.0 AND zWLnLtT LE 295.0) IF zInd[0] EQ -1 THEN BEGIN MESSAGE, 'Cannot find altitude index in WACCM' ENDIF ELSE BEGIN iZ = zInd[0] ENDELSE tNZSub[iLn,iT] = tN[iLn,iLt,iZ,iT] tNZSub_shift[iLn,iT] = tN_shift[iLn,iLt,iZ,iT] twvtotal1to3Sub[iLn,iT] = twvtotal1to3[iLn,iLt,iT] twvtotal1to3Sub_shift[iLn,iT] = twvtotal1to3_shift[iLn,iLt,iT] twvtotal1to6Sub[iLn,iT] = twvtotal1to6[iLn,iLt,iT] twvtotal1to6Sub_shift[iLn,iT] = twvtotal1to6_shift[iLn,iLt,iT] twvtotal1to9Sub[iLn,iT] = twvtotal1to9[iLn,iLt,iT] twvtotal1to9Sub_shift[iLn,iT] = twvtotal1to9_shift[iLn,iLt,iT] ENDFOR ENDFOR ; ; ;March Equator Mean Temperature ;March 20N Gravity Wave parameterization High Resolution Temperature ; ; SET_PLOT, 'PS' ; iLt = 58 ;20 North ; device, /color, /land, bits_per_pixel=8,filename='mtm_UT00-UT23_MarchMean_Equator_Tcomp_lon.ps' ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' LOADCT, 39, ncolors = 16 ; dimsize = size(twzonalmeandev0) xTitle = 'Longitude (Degrees)' yTitle = 'Temperature (K)' ; ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHR_20N_zmd_tcomp.ps' Title0 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT00' ; LOADCT, 39, ncolors = 10 plot, lon_shift, tNZSub_shift[*,0], xrange=[-90.0,90.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; Title0 = ' WAM Temperature Deviation from Zonal Mean March GWHR 20N 285km UT00' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,0])+550., color=14, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,0])+550., color=9, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,0])+550., color=3, thick=6 ; nTAllW = dimsize[3] axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title0 = ' WAM Temperature Reconstruction March Mean Equator 285km UT00' ; plot, lonWAM_shift, tnWAM_shift[*,46,142,0],xrange=[-90.0,90.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,0])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,0])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,0])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title1 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT01' ; LOADCT, 39, ncolors = 10 plot, lon_shift, tNZSub_shift[*,1], xrange=[-105.0,75.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title1,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; Title0 = ' WAM Temperature Deviation from Zonal Mean March GWHR 20N 285km UT00' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,1])+550., color=14, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,1])+550., color=9, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,1])+550., color=3, thick=6 ; nTAllW = dimsize[3] axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title1 = ' WAM Temperature Reconstruction March Mean Equator 285km UT01' ; plot, lonWAM_shift,tnWAM_shift[*,46,142,1],xrange=[-105.0,75.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title1,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,1])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,1])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,1])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title2 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT02' ; LOADCT, 39, ncolors = 10 plot, lon_shift, tNZSub_shift[*,2], xrange=[-120.0,60.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title2,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; Title0 = ' WAM Temperature Deviation from Zonal Mean March GWHR 20N 285km UT00' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,2])+550., color=14, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,2])+550., color=9, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,2])+550., color=3, thick=6 ; nTAllW = dimsize[3] axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title2 = ' WAM Temperature Reconstruction March Mean Equator 285km UT02' ; plot, lonWAM_shift, tnWAM_shift[*,46,142,2],xrange=[-120.0,60.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title2,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,2])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,2])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,2])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title3 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT03' ; plot, lon_shift, tNZSub_shift[*,3], xrange=[-135.0,45.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,3])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,3])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,3])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title3 = ' WAM Temperature Reconstruction March Mean Equator 285km UT03' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,3], xrange=[-135.0,45.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,3])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,3])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,3])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title4 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT04' ; plot, lon_shift, tNZSub_shift[*,4], xrange=[-150.0,30.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title4,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,4])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,4])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,4])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title4 = ' WAM Temperature Reconstruction March Mean Equator 285km UT04' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,4], xrange=[-150.0,30.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title4,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,4])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,4])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,4])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title5 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT05' ; plot, lon_shift, tNZSub_shift[*,5], xrange=[-165.0,15.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title5,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,5])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,5])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,5])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title5 = ' WAM Temperature Reconstruction March Mean Equator 285km UT05' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,5], xrange=[-165.0,15.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title5,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,5])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,5])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,5])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title6 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT06' ; plot, lon_shift, tNZSub_shift[*,6],xrange=[-180.0,0.0],yrange=[460,630],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,6])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,6])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,6])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title6 = ' WAM Temperature Reconstruction March Mean Equator 285km UT06' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,6], xrange=[-180.0,0.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,6])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,6])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,6])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title7 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT07' ; plot, lon, tNZSub[*,7],xrange=[165.0,345.0],yrange=[450,630],xstyle=9,ystyle=1,subtitle=Title7,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,7])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,7])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,7])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title7 = ' WAM Temperature Reconstruction March Mean Equator 285km UT07' ; plot, lonWAM, tnWAM[*,iLt,142,7], xrange=[165.0,345.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title7,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,7])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,7])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,7])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title8 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT08' ; plot, lon, tNZSub[*,8],xrange=[150.0,330.0],yrange=[460,630],xstyle=9,ystyle=1,subtitle=Title8,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,8])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,8])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,8])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title8 = ' WAM Temperature Reconstruction March Mean Equator 285km UT08' ; plot, lonWAM, tnWAM[*,iLt,142,8], xrange=[150.0,330.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title8,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,8])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,8])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,8])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title9 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT09' ; plot, lon, tNZSub[*,9], xrange=[135.0,315.0],yrange=[460,620],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,9])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,9])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,9])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title9 = ' WAM Temperature Reconstruction March Mean Equator 285km UT09' ; plot, lonWAM, tnWAM[*,iLt,142,9], xrange=[135.0,315.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,9])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,9])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,9])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title10 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT10' ; plot, lon, tNZSub[*,10], xrange=[120.0,300.0],yrange=[450,620],xstyle=9,ystyle=1,subtitle=Title10,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,10])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,10])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,10])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title10 = ' WAM Temperature Reconstruction March Mean Equator 285km UT10' ; plot, lonWAM, tnWAM[*,iLt,142,10], xrange=[120.0,300.0],yrange=[730,890],xstyle=9,ystyle=1,subtitle=Title10,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,10])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,10])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,10])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title11 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT11' ; plot, lon, tNZSub[*,11], xrange=[105.0,285.0],yrange=[460,620],xstyle=9,ystyle=1,subtitle=Title11,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,11])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,11])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,11])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title11 = ' WAM Temperature Reconstruction March Mean Equator 285km UT11' ; plot, lonWAM, tnWAM[*,iLt,142,11], xrange=[105.0,285.0],yrange=[740,890],xstyle=9,ystyle=1,subtitle=Title11,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,11])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,11])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,11])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title12 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT12' ; plot, lon, tNZSub[*,12], xrange=[90.0,270.0],yrange=[480,620],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,12])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,12])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,12])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title12 = ' WAM Temperature Reconstruction March Mean Equator 285km UT12' ; plot, lonWAM, tnWAM[*,iLt,142,12], xrange=[90.0,270.0],yrange=[750,890],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,12])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,12])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,12])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title13 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT13' ; plot, lon, tNZSub[*,13], xrange=[75.0,255.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title13,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,13])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,13])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,13])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title13 = ' WAM Temperature Reconstruction March Mean Equator 285km UT13' ; plot, lonWAM, tnWAM[*,iLt,142,13], xrange=[75.0,255.0],yrange=[750,890],xstyle=9,ystyle=1,subtitle=Title13,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,13])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,13])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,13])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 Title14 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT14' ; plot, lon, tNZSub[*,14], xrange=[60.0,240.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title14,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,14])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,14])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,14])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title14 = ' WAM Temperature Reconstruction March Mean Equator 285km UT14' ; plot, lonWAM, tnWAM[*,iLt,142,14], xrange=[60.0,240.0],yrange=[750,890],xstyle=9,ystyle=1,subtitle=Title14,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,14])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,14])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,14])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title15 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT15' ; plot, lon, tNZSub[*,15], xrange=[45.0,225.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,15])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,15])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,15])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title15 = ' WAM Temperature Reconstruction March Mean Equator 285km UT15' ; plot, lonWAM, tnWAM[*,iLt,142,15], xrange=[45.0,225.0],yrange=[760,910],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,15])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,15])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,15])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title16 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT16' ; plot, lon, tNZSub[*,16], xrange=[30.0,210.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title16,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,16])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,16])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,16])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title16 = ' WAM Temperature Reconstruction March Mean Equator 285km UT16' ; plot, lonWAM, tnWAM[*,iLt,142,16], xrange=[30.0,210.0],yrange=[760,910],xstyle=9,ystyle=1,subtitle=Title16,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,16])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,16])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,16])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title17 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT17' ; plot, lon, tNZSub[*,17], xrange=[15.0,195.0],yrange=[460,610],xstyle=9,ystyle=1,subtitle=Title17,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon, reform(twvtotal1to3Sub[*,17])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon, reform(twvtotal1to6Sub[*,17])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon, reform(twvtotal1to9Sub[*,17])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title17 = ' WAM Temperature Reconstruction March Mean Equator 285km UT17' ; plot, lonWAM, tnWAM[*,iLt,142,17], xrange=[15.0,195.0],yrange=[760,910],xstyle=9,ystyle=1,subtitle=Title17,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM, reform(twvtotal1to3WAM[*,46,17])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM, reform(twvtotal1to6WAM[*,46,17])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM, reform(twvtotal1to9WAM[*,46,17])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title18 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT18' ; plot, lon_shift, tNZSub_shift[*,18], xrange=[0.0,180.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,18])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,18])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,18])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title18 = ' WAM Temperature Reconstruction March Mean Equator 285km UT18' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,18], xrange=[0.0,180.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,18])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,18])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,18])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title19 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT19' ; plot, lon_shift, tNZSub_shift[*,19], xrange=[-15.0,165.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title19,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,19])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,19])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,19])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title19 = ' WAM Temperature Reconstruction March Mean Equator 285km UT19' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,19], xrange=[-15.0,165.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title19,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,19])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,19])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,19])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title20 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT20' ; plot, lon_shift, tNZSub_shift[*,20], xrange=[-30.0,150.0],yrange=[470,620],xstyle=9,ystyle=1,subtitle=Title20,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,20])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,20])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,20])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title20 = ' WAM Temperature Reconstruction March Mean Equator 285km UT20' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,20], xrange=[-30.0,150.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title20,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,20])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,20])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,20])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title21 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT21' ; plot, lon_shift, tNZSub_shift[*,21], xrange=[-45.0,135.0],yrange=[480,610],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,21])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,21])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,21])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title21 = ' WAM Temperature Reconstruction March Mean Equator 285km UT21' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,21], xrange=[-45.0,135.0],yrange=[760,900],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,21])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,21])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,21])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title22 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT22' ; plot, lon_shift, tNZSub_shift[*,22], xrange=[-60.0,120.0],yrange=[470,610],xstyle=9,ystyle=1,subtitle=Title22,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,22])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,22])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,22])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title22 = ' WAM Temperature Reconstruction March Mean Equator 285km UT22' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,22], xrange=[-60.0,120.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title22,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,22])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,22])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,22])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title23 = ' WACCM-X Temperature Reconstruction March Mean Equator 285km UT23' ; plot, lon_shift, tNZSub_shift[*,23], xrange=[-75.0,105.0],yrange=[480,610],xstyle=9,ystyle=1,subtitle=Title23,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lon_shift, reform(twvtotal1to3Sub_shift[*,23])+550., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lon_shift, reform(twvtotal1to6Sub_shift[*,23])+550., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lon_shift, reform(twvtotal1to9Sub_shift[*,23])+550., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' Title23 = ' WAM Temperature Reconstruction March Mean Equator 285km UT23' ; plot, lonWAM_shift, tnWAM_shift[*,iLt,142,23], xrange=[-75.0,105.0],yrange=[750,900],xstyle=9,ystyle=1,subtitle=Title23,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' oplot, lonWAM_shift, reform(twvtotal1to3WAM_shift[*,46,23])+835., color=14, thick=6 ; nTAllW = dimsize[3] oplot, lonWAM_shift, reform(twvtotal1to6WAM_shift[*,46,23])+835., color=9, thick=6 ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 oplot, lonWAM_shift, reform(twvtotal1to9WAM_shift[*,46,23])+835., color=3, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' xyouts, /NORMAL, 0.65, 0.8, 'Migrating Tides s=2-3 ', color=14, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.76, 'Migrating Tides s=2-6 ', color=9, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.72, 'Migrating Tides s=2-9 ', color=3, charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 xyouts, /NORMAL, 0.65, 0.68, 'All', charsize=1.5 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 device,/close ; oplot, lonworig, reform(uwvtotal[*,58,7,5]), color=9, thick=6 set_plot, 'x' ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT18' ; plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,6]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav' ; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT21' ; plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,7]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ;; twzonalmeandevorig21 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav' ; ; device,/close ; set_plot, 'x' ; ;March Equator Gravity Wave parameterization High Resolution Temperature ; ; iLt = 47 ;20 North ; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHR_Equator_zmd_tcomp.ps' ; LOADCT, 39, ncolors = 16 ; Title0 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT00' ; plot, lonw, twzonalmeandev0[*,iLt,0], xrange=[-90.0,90.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev0[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,0]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,0]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' ; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_Equator_T.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT03' ; plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,1]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,1]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_Equator_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT06' ; plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev6[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,2]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,2]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev6 = twzonalmeandev ;; twzonalmeandevorig6 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_Equator_T.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT09' ; plot, lonworig, twzonalmeandevorig9[*,iLt,0],xrange=[135.0,315.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig9[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,3]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,iLt,7,3]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev9 = twzonalmeandev ;; twzonalmeandevorig9 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_Equator_T.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT12' ; plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig12[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,4]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,iLt,7,4]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev12 = twzonalmeandev ;; twzonalmeandevorig12 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_Equator_T.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March Equator 285km UT15' ; plot,lonworig,twzonalmeandevorig15[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig15[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,5]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,iLt,7,5]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_Equator_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT18' ; plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,6]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,6]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_Equator_T.sav' ; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_31N_FullLon.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_Equator_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March GWHR Equator 285km UT21' ; plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,7]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,7]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ;; twzonalmeandevorig21 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_Equator_T.sav' ; ; device,/close ; set_plot, 'x' ; ;March 20N 3.1.9 Mean Temperature ; ; iLt = 58 ;20 North ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_March_20N_zmd_tcomp.ps' ; LOADCT, 39, ncolors = 16 ; Title0 = ' Temperature Deviation from Zonal Mean March 20N 285km UT00' ; plot, lonw, twzonalmeandev0[*,iLt,0], xrange=[-90.0,90.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev0[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_31N_FullLon.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT03_March_20N_T.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March 20N 285km UT03' ; plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,1]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT06_March_20N_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March 20N 285km UT06' ; plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev6[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,2]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev6 = twzonalmeandev ;; twzonalmeandevorig6 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT09_March_20N_T.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March 20N 285km UT09' ; plot, lonworig, twzonalmeandevorig9[*,iLt,0],xrange=[135.0,315.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig9[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,3]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,3]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev9 = twzonalmeandev ;; twzonalmeandevorig9 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT12_March_20N_T.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March 20N 285km UT12' ; plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig12[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,4]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev12 = twzonalmeandev ;; twzonalmeandevorig12 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT15_March_20N_T.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15' ; plot,lonworig,twzonalmeandevorig15[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig15[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,5]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT18_March_20N_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March 20N 285km UT18' ; plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,6]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_31N_FullLon_20N_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT21_March_20N_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March 20N 285km UT21' ; plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,7]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ;; twzonalmeandevorig21 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N_T.sav' ; ; device,/close ; set_plot, 'x' ; ;March Equator 3.1.9 Mean Temperature ; ; iLt = 47 ;Equator ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; SET_PLOT, 'PS' ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_March_Equator_zmd_tcomp.ps' ; LOADCT, 39, ncolors = 16 ; Title0 = ' Temperature Deviation from Zonal Mean March Equator 285km UT00' ; plot, lonw, twzonalmeandev0[*,iLt,0], xrange=[-90.0,90.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev0[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,0]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,47,7,0]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_March_31N_FullLon.sav' ; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT03_March_20N_T.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March Equator 285km UT03' ; plot, lonw, twzonalmeandev3[*,iLt,0], xrange=[-135.0,45.0],yrange=[-110,60],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev3[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,1]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,47,7,1]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_March_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT06_March_20N_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March Equator 285km UT06' ; plot, lonw, twzonalmeandev6[*,iLt,0], xrange=[-180.0,0.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev6[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,2]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,47,7,2]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev6 = twzonalmeandev ;; twzonalmeandevorig6 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_March_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT09_March_20N_T.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March Equator 285km UT09' ; plot, lonworig, twzonalmeandevorig9[*,iLt,0],xrange=[135.0,315.0],yrange=[-110,70],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig9[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,47,7,3]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,47,7,3]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev9 = twzonalmeandev ;; twzonalmeandevorig9 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_March_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT12_March_20N_T.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March Equator 285km UT12' ; plot, lonworig,twzonalmeandevorig12[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig12[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,47,7,4]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,47,7,4]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev12 = twzonalmeandev ;; twzonalmeandevorig12 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_March_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT15_March_20N_T.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March Equator 285km UT15' ; plot,lonworig,twzonalmeandevorig15[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig15[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,47,7,5]), color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,47,7,5]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_March_31N_FullLon_Equator_T.sav' ; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT18_March_20N_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March Equator 285km UT18' ; plot, lonw, twzonalmeandev18[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,60],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev18[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,6]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,47,7,6]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_March_31N_FullLon_Equator_T.sav' ;; ; restore, '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N.sav' ; dimsize = size(twzonalmeandev0) ; nTAllW = dimsize[3] ; ;; device, /color, /land, filename='idl_tmt_UT21_March_20N_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March Equator 285km UT21' ; plot, lonw, twzonalmeandev21[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,50],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev21[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,47,7,7]), color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,47,7,7]), color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ;; twzonalmeandevorig21 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_March_20N_T.sav' ; ; device,/close ; set_plot, 'x' ; ; SET_PLOT, 'PS' ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHRMean_20N_T_tcomp.ps' ; LOADCT, 39, ncolors = 16 ; ; Title0 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT00' ; plot, lonw, tn_shift[*,iLt,7,0], xrange=[-90.0,90.0],yrange=[450,590],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0])+530., color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' ; ;; device,/close ;; set_plot, 'x' ;; ;; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT03' ; plot, lonw, tn_shift[*,iLt,7,1], xrange=[-135.0,45.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1])+530., color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,1])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT06' ; plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,2])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev6 = twzonalmeandev ;; twzonalmeandevorig6 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_20N_T.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT09' ; plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,3])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,3])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev9 = twzonalmeandev ;; twzonalmeandevorig9 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_20N_T.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT12' ; plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,4])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev12 = twzonalmeandev ;; twzonalmeandevorig12 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_20N_T.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT15' ; plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,5])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT18' ; plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,6])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March Mean GWHR 20N 285km UT21' ; plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,7])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ;; twzonalmeandevorig21 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav' ; ; device,/close ; set_plot, 'x' ; ;March Equator GWHR Mean Temperature ; ; iLt = 47 ;Equator ; ; SET_PLOT, 'PS' ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchGWHRMean_Equator_T_tcomp.ps' ; LOADCT, 39, ncolors = 16 ; ; Title0 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT00' ; plot, lonw, tn_shift[*,iLt,7,0], xrange=[-90.0,90.0],yrange=[450,590],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,0])+530., color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,0])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_Equator_T.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT03' ; plot, lonw, tn_shift[*,iLt,7,1], xrange=[-135.0,45.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,1])+530., color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,1])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_Equator_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT06' ; plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,2])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,2])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev6 = twzonalmeandev ;; twzonalmeandevorig6 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_Equator_T.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT09' ; plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,3])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,iLt,7,3])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev9 = twzonalmeandev ;; twzonalmeandevorig9 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_Equator_T.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT12' ; plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,4])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,iLt,7,4])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev12 = twzonalmeandev ;; twzonalmeandevorig12 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_Equator_T.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT15' ; plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,iLt,7,5])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,iLt,7,5])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_Equator_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT18' ; plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,6])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,6])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_Equator_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March Mean GWHR Equator 285km UT21' ; plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,iLt,7,7])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,iLt,7,7])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_Equator_T.sav' ; ; device,/close ; set_plot, 'x' ; ; SET_PLOT, 'PS' ; ; device, /color, /land, bits_per_pixel=8, filename='idl_mtm_UT00-UT21_MarchMean_20N_T_tcomp.ps' ; LOADCT, 39, ncolors = 16 ; ; Title0 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT00' ; plot, lonw, tn_shift[*,iLt,7,0], xrange=[-90.0,90.0],yrange=[450,590],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,0])+530., color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,0])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' ; ;; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT03' ; plot, lonw, tn_shift[*,iLt,7,1], xrange=[-135.0,45.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,1])+530., color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,1])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT06' ; plot, lonw, tn_shift[*,iLt,7,2], xrange=[-180.0,0.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,2])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,2])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev6 = twzonalmeandev ;; twzonalmeandevorig6 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_20N_T.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT09' ; plot, lonworig, tn[*,iLt,7,3],xrange=[135.0,315.0],yrange=[440,590],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,3])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,3])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev9 = twzonalmeandev ;; twzonalmeandevorig9 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_20N_T.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT12' ; plot, lonworig,tn[*,iLt,7,4],xrange=[90.0,270.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,4])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,4])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev12 = twzonalmeandev ;; twzonalmeandevorig12 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon_20N_T.sav' ;; ; xTitleOrig = ' Longitude (Degrees East)' ;; device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_20N_T.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT15' ; plot,lonworig,tn[*,iLt,7,5],xrange=[45.0,225.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonworig, tn[*,iLt,iPlt] ; oplot, lonworig, reform(compd_sd_td_qd[*,58,7,5])+530, color=14, thick=6 ; oplot, lonworig, reform(uwvtotal[*,58,7,5])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev15 = twzonalmeandev ;; twzonalmeandevorig15 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT18' ; plot, lonw, tn_shift[*,iLt,7,6], xrange=[0.0,180.0],yrange=[460,600],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,6])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,6])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev18 = twzonalmeandev ;; twzonalmeandevorig18 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon_20N_T.sav' ;; ;; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March Mean 20N 285km UT21' ; plot, lonw, tn_shift[*,iLt,7,7], xrange=[-45.0,135.0],yrange=[460,590],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ;; for iPlt = 1,nTAllW-1 DO oplot, lonw, tn[*,iLt,iPlt] ; oplot, lonw, reform(compd_sd_td_qd_shift[*,58,7,7])+530, color=14, thick=6 ; oplot, lonw, reform(uwvtotal_shift[*,58,7,7])+530., color=9, thick=6 ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ;; twzonalmeandev21 = twzonalmeandev ;; twzonalmeandevorig21 = twzonalmeandevorig ;; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = '/local/i/joemci/SW/IDL/MTM/MarchSav/waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav' ; ; device,/close ; set_plot, 'x' ; return end