print, 'Running waccm_mtm_savin_1time_T' ; ; Set the pressure levels to process ; nZ = 1 zNeeded = [285.0] ; ; Set Latitude to process ; ; latNeeded = 20.0 ; ; Read first WACCM file to get dimensions of data time read ; ; waccmFilename = '/aim/d/waccm/sminyear/12/wa319_2x_sminyear.cam2.h1.1951-12-01-00000.nc' waccmFilename = '/aim/d/waccm/wa319_TIPHYS/MTM/WAX3548GWSMnQUP50_2x.cam2.h1.1954-03-01-00000.nc' ; print, 'Reading first input file to get dimensions ' ncdf_waccm_read,waccmFilename,lonW,latW,lev3D,timeOutW,nXW,nYW,nZW,nTW,nameCnst,$ GEOPTH=zW,XWind=uNW,YWind=vNW,ZWind=wNW,TEMP=tNW print, 'Restoring input WACCM data' ; RESTORE, 'waccm_mtm_UT0.sav' ; RESTORE, 'waccm_mtm_UT0_March.sav' ; RESTORE, 'waccm_mtm_UT3_March.sav' ; RESTORE, 'waccm_mtm_UT6_March.sav' ; RESTORE, 'waccm_mtm_UT9_March.sav' ; RESTORE, 'waccm_mtm_UT12_March.sav' ; RESTORE, 'waccm_mtm_UT15_March.sav' ; RESTORE, 'waccm_mtm_UT18_March.sav' ; RESTORE, 'waccm_mtm_UT21_March.sav' ; ; RESTORE, 'waccm_mtm_UT0_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT3_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT6_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT9_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT12_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT15_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT18_Dec_30S.sav' ; RESTORE, 'waccm_mtm_UT21_Dec_30S.sav' ; ; RESTORE, 'waccm_mtm_UT0_MarchGWHR.sav' RESTORE, 'waccm_mtm_UT3_MarchGWHR.sav' ; RESTORE, 'waccm_mtm_UT6_MarchGWHR.sav' ; RESTORE, 'waccm_mtm_UT9_MarchGWHR.sav' ; RESTORE, 'waccm_mtm_UT12_MarchGWHR.sav' ; RESTORE, 'waccm_mtm_UT15_MarchGWHR.sav' ; RESTORE, 'waccm_mtm_UT18_MarchGWHR.sav' ; RESTORE, 'waccm_mtm_UT21_MarchGWHR.sav' ; ; RESTORE, 'waccm_mtm_UT0_Sept.sav' ; RESTORE, 'waccm_mtm_UT3_Sept.sav' ; RESTORE, 'waccm_mtm_UT6_Sept.sav' ; RESTORE, 'waccm_mtm_UT9_Sept.sav' ; RESTORE, 'waccm_mtm_UT12_Sept.sav' ; RESTORE, 'waccm_mtm_UT15_Sept.sav' ; RESTORE, 'waccm_mtm_UT18_Sept.sav' ; RESTORE, 'waccm_mtm_UT21_Sept.sav' ; ; RESTORE, 'waccm_mtm_UT0_June.sav' ; RESTORE, 'waccm_mtm_UT3_June.sav' ; RESTORE, 'waccm_mtm_UT6_June.sav' ; RESTORE, 'waccm_mtm_UT9_June.sav' ; RESTORE, 'waccm_mtm_UT12_June.sav' ; RESTORE, 'waccm_mtm_UT15_June.sav' ; RESTORE, 'waccm_mtm_UT18_June.sav' ; RESTORE, 'waccm_mtm_UT21_June.sav' ; ; tNAllW = tNAllW1[*,*,*,0:29] ; zNAllW = zNAllW1[*,*,*,0:29] ; tNAllW = tNAllW2[*,*,*,0:29] ; zNAllW = zNAllW2[*,*,*,0:29] ; tNAllW = tNAllW3[*,*,*,0:29] ; zNAllW = zNAllW3[*,*,*,0:29] ; tNAllW = tNAllW4[*,*,*,0:29] ; zNAllW = zNAllW4[*,*,*,0:29] ; tNAllW = tNAllW5[*,*,*,0:29] ; zNAllW = zNAllW5[*,*,*,0:29] ; tNAllW = tNAllW6[*,*,*,0:29] ; zNAllW = zNAllW6[*,*,*,0:29] ; tNAllW = tNAllW7[*,*,*,0:29] ; zNAllW = zNAllW7[*,*,*,0:29] ; tNAllW = tNAllW8[*,*,*,0:29] ; zNAllW = zNAllW8[*,*,*,0:29] ; ; tNAllW = tNAllW1[*,*,*,0:30] ; zNAllW = zNAllW1[*,*,*,0:30] tNAllW = tNAllW2[*,*,*,0:30] zNAllW = zNAllW2[*,*,*,0:30] ; tNAllW = tNAllW3[*,*,*,0:30] ; zNAllW = zNAllW3[*,*,*,0:30] ; tNAllW = tNAllW4[*,*,*,0:30] ; zNAllW = zNAllW4[*,*,*,0:30] ; tNAllW = tNAllW5[*,*,*,0:30] ; zNAllW = zNAllW5[*,*,*,0:30] ; tNAllW = tNAllW6[*,*,*,0:30] ; zNAllW = zNAllW6[*,*,*,0:30] ; tNAllW = tNAllW7[*,*,*,0:30] ; zNAllW = zNAllW7[*,*,*,0:30] ; tNAllW = tNAllW8[*,*,*,0:30] ; zNAllW = zNAllW8[*,*,*,0:30] ; ; Find altitude level to process in WACCM data ; tZSubW = FLTARR(nXW,nYW,nTAllW) tWMonthMean = FLTARR(nXW,nYW) tWMonthMeanOrig = FLTARR(nXW,nYW) tWZonalMean = FLTARR(nYW,nTAllW) tWZonalMeanOrig = FLTARR(nYW,nTAllW) tWZonalMeanDev = FLTARR(nXW,nYW,nTAllW) tWZonalMeanDevOrig = FLTARR(nXW,nYW,nTAllW) FOR iLt=0,nYW-1 DO BEGIN IF iLt GT 0 THEN lonW = lonWOrig FOR iLn=0,nXW-1 DO BEGIN FOR iT = 0, nTAllW-1 DO BEGIN ; ; Find altitude index for this time ; zWLnLtTIn = REFORM(zNAllW[iLn,iLt,*,iT]) zWLnLtT = zWLnLtTIn * (1. + zWLnLtTIn/6370.0) FOR iZW = 1, nZW-1 DO BEGIN ; print, 'index is ', iLt, iLn, iT, iZW, zWLnLtT[iZW] , zWLnLtTIn[iZW] IF ZWLnLtT[iZW] LT 285.0 THEN BEGIN tDiffAll = zWLnLtT(iZW-1) - zWLnLtT(iZW) tDiffLowerWt = (zWLnLtT(iZW-1) - 285.0) / tDiffAll tDiffUpperWt = (285.0 - zWLnLtT(iZW)) / tDiffAll tZSubW(iLn,iLt,iT) = tNAllW(iLn,iLt,iZW,iT)*tDiffLowerWt + tNAllW(iLn,iLt,iZW-1,iT)*tDiffUpperWt iZ = iZW BREAK ENDIF ELSE IF ZWLnLtT[iZW] EQ 285.0 THEN BEGIN tZSubW(iLn,iLt,iT) = REFORM(tNAllW(iLn,iLt,iZW,iT)) iZ = iZW BREAK ENDIF ENDFOR ; 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 ; ; Here iZ is a single index ; ; tZSubW(iLn,iLt,iT) = REFORM(tNAllW(iLn,iLt,iZ,iT)) ENDFOR ;Time ENDFOR ;Longitude ; ; Need to convert from 0-360 longitude grid to -180.0 to 180.0 longitude grid ; nLonTemp = nXW+nXW/2 tZSubWTemp = FLTARR(nLonTemp,nYW,nTAllW) tZSubWOrig = tZSubW tZSubWTemp[0:nXW-1,iLT,*] = tZSubW[0:nXW-1,iLt,*] tZSubWTemp[nXW:nLonTemp-1,iLt,*] = tZSubW[0:nXW/2-1,iLt,*] tZSubW[0:nXW-1,iLt,*] = tZSubWTemp[nXW/2:nLonTemp-1,iLt,*] lonWOrig = lonW lonW = lonW - 180.0 ; ; Find month mean ; FOR iLn=0,nXW-1 DO BEGIN tWMonthMeanOrig[iLn,iLt] = MEAN(tZSubWOrig[iLn,iLt,*]) ; tWZonalMean[iT] = (TOTAL(tZSubW[0:iLonMax,iT])*(FLOAT(iLonMax)/FLOAT(iLonMin-iLonMax)) + $ ; TOTAL(tZSubW[iLonMin:143,iT])*((FLOAT(nXW)-FLOAT(iLonMax-1))/FLOAT(iLonMin-iLonMax))) tWMonthMean[iLn,iLt] = MEAN(tZSubW[iLn,iLt,*]) ; print, 'iLt, iT, twzonalmean[iLt,iT] ', iLt, iT, twzonalmean[iLt,iT] ENDFOR ;lon,lat,level,time zAmp ; zAmp = FLTARR(nXW,nYW,nP,nTAllW) ; tAmp = FLTARR(nXW,nYW,nP,nTAllW) ; ; Set zonal mean for appropriate longitude range ; ; iLonMin = 36 ; iLonMax = 108 ; iLonMin = 18 ; iLonMax = 90 ; iLonMin = 0 ; iLonMax = 72 ; iLonMin = 54 ; iLonMax = 126 ; iLonMin = 36 ; iLonMax = 108 ; iLonMin = 18 ; iLonMax = 90 ; iLonMin = 72 ; iLonMax = 143 ; iLonMin = 54 ; iLonMax = 126 iLonMin = 0 iLonMax = 143 ; ; Find zonal mean for appropriate longitude range ; FOR iT=0,nTAllW-1 DO BEGIN tWZonalMeanOrig[iLt,iT] = MEAN(tZSubWOrig[iLonMin:iLonMax,iLt,iT]) ; tWZonalMean[iT] = (TOTAL(tZSubW[0:iLonMax,iT])*(FLOAT(iLonMax)/FLOAT(iLonMin-iLonMax)) + $ ; TOTAL(tZSubW[iLonMin:143,iT])*((FLOAT(nXW)-FLOAT(iLonMax-1))/FLOAT(iLonMin-iLonMax))) tWZonalMean[iLt,iT] = MEAN(tZSubW[iLonMin:iLonMax,iLt,iT]) ; print, 'iLt, iT, twzonalmean[iLt,iT] ', iLt, iT, twzonalmean[iLt,iT] ENDFOR ; ; Need to calculate deviation from zonal mean for each day and time ; ; FOR iLn=0,nXW-1 DO BEGIN ; ; FOR iT=0,nTAllW-1 DO BEGIN ; ; tWZonalMeanDev[iLn,iLt,iT] = tZSubW[iLn,iLt,iT] - tWZonalMean[iLt,iT] ; tWZonalMeanDevOrig[iLn,iLt,iT] = tZSubWOrig[iLn,iLt,iT] - tWZonalMeanOrig[iLt,iT] ; print, 'iLt, iT2, twzonalmean[iLt,iT] ', iLt, iT, twzonalmean[iLt,iT] ; ; ; ENDFOR ; ; ENDFOR ENDFOR ;Latitude iLt = 58 ;20 degrees north latitude ; iLt = 32 ;30 degrees south latitude set_plot, 'ps' ; device, /color, filename='idl_tmt_UT00-UT06.ps' xTitle = ' Longitude (Degrees)' yTitle = 'Temperature Deviation from Zonal Mean (K) LTimes = ['18','20','22','0','2','4','6'] !P.CHARSIZE = 1.5 ; ;March 20N High Resolution Temperature ; iLt = 58 ;20 North ; device, /color, /land, filename='idl_tmt_UT00_March_20N_T.ps' ; Title0 = ' Temperature March 20N 285km UT00' ; plot, lonw, tzsubw[*,iLt,0], xrange=[-90.0,90.0],yrange=[420,600],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, tzsubw[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; tzsubw0 = tzsubw ; tzsubworig0 = tzsubworig ; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = 'waccm_mtm_zmd_UT00_March_20N_T.sav' ; device, /color, /land, filename='idl_tmt_UT03_March_20N_T.ps' Title3 = ' Temperature March 20N 285km UT03' plot, lonw, tzsubw[*,iLt,0], xrange=[-135.0,45.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] for iPlt = 1,nTAllW-1 DO oplot, lonw, tzsubw[*,iLt,iPlt] axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' tzsubw3 = tzsubw tzsubworig3 = tzsubworig save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = 'waccm_mtm_zmd_UT03_March_20N_T.sav' ; ;March 20N Gravity Wave parameterization High Resolution Temperature ; iLt = 58 ;20 North ; device, /color, /land, filename='idl_tmt_UT00_MarchGWHR_20N_T.ps' ; Title0 = ' Temperature March GWHR 20N 285km UT00' ; plot, lonw, tzsubw[*,iLt,0], xrange=[-90.0,90.0],yrange=[420,600],xstyle=9,ystyle=1,subtitle=Title0,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, tzsubw[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; tzsubw0 = tzsubw ; tzsubworig0 = tzsubworig ; save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = 'waccm_mtm_zmd_UT00_MarchGWHR_20N_T.sav' ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N_T.ps' Title3 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT03' plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title3,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' tzsubw3 = tzsubw tzsubworig3 = tzsubworig save, lonw, lonworig, latw, tzsubw0, tzsubworig0, file = 'waccm_mtm_zmd_UT03_MarchGWHR_20N_T.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N_T.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title6,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_MarchGWHR_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 GWHR 20N 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title9,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_MarchGWHR_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 GWHR 20N 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title12,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_MarchGWHR_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 20N 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title15,xtitle=xTitleorig,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_MarchGWHR_20N_T.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N_T.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title18,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_MarchGWHR_20N_T.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N_T.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-80,100],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_MarchGWHR_20N_T.sav' device,/close set_plot, 'x' ; ;June 20N Full Longitude Mean ; ; iLt = 58 ;20 degrees north latitude ; ; device, /color, /land, filename='idl_tmt_UT00_June_20N_FullLon.ps' ; Title0 = ' Temperature Deviation from Zonal Mean June 20N 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev0 = twzonalmeandev ; twzonalmeandevorig0 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_June_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT03_June_20N_FullLon.ps' ; Title3 = ' Temperature Deviation from Zonal Mean June 20N 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev3 = twzonalmeandev ; twzonalmeandevorig3 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev3, twzonalmeandevorig3, file = 'waccm_mtm_zmd_UT03_June_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_June_20N_FullLon.ps' ; Title6 = ' Temperature Deviation from Zonal Mean June 20N 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_June_20N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_June_20N_FullLon.ps' ; Title9 = ' Temperature Deviation from Zonal Mean June 20N 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_June_20N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_June_20N_FullLon.ps' ; Title12 = ' Temperature Deviation from Zonal Mean June 20N 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_June_20N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_June_20N_FullLon.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_June_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_June_20N_FullLon.ps' ; Title18 = ' Temperature Deviation from Zonal Mean June 20N 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_June_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_June_20N_FullLon.ps' ; Title21 = ' Temperature Deviation from Zonal Mean June 20N 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_June_20N_FullLon.sav' ; ; device,/close ; set_plot, 'x' ; ;June Equator Full Longitude Mean ; ; set_plot, 'ps' ; ; iLt = 47 ; ; device, /color, /land, filename='idl_tmt_UT00_June_Eq_FullLon.ps' ; Title0 = ' Temperature Deviation from Zonal Mean June Equator 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT03_June_Eq_FullLon.ps' ; Title3 = ' Temperature Deviation from Zonal Mean June Equator 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT06_June_Eq_FullLon.ps' ; Title6 = ' Temperature Deviation from Zonal Mean June Equator 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_June_Eq_FullLon.ps' ; Title9 = ' Temperature Deviation from Zonal Mean June Equator 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_June_Eq_FullLon.ps' ; Title12 = ' Temperature Deviation from Zonal Mean June Equator 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_June_Eq_FullLon.ps' ; Title15 = ' Temperature Deviation from Zonal Mean June Equator 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT18_June_Eq_FullLon.ps' ; Title18 = ' Temperature Deviation from Zonal Mean June Equator 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT21_June_Eq_FullLon.ps' ; Title21 = ' Temperature Deviation from Zonal Mean June Equator 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,80],xstyle=9,ystyle=1,subtitle=Title21,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3] ; axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time' ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device,/close ; set_plot, 'x' ; ;September 20N Full Longitude Mean ; device, /color, /land, filename='idl_tmt_UT00_Sept_20N_FullLon.ps' ; Title0 = ' Temperature Deviation from Zonal Mean Sept 20N 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev0 = twzonalmeandev ; twzonalmeandevorig0 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_Sept_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT03_Sept_20N_FullLon.ps' ; Title3 = ' Temperature Deviation from Zonal Mean September 20N 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev3 = twzonalmeandev ; twzonalmeandevorig3 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev3, twzonalmeandevorig3, file = 'waccm_mtm_zmd_UT03_Sept_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_Sept_20N_FullLon.ps' ; Title6 = ' Temperature Deviation from Zonal Mean September 20N 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_Sept_20N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_Sept_20N_FullLon.ps' ; Title9 = ' Temperature Deviation from Zonal Mean September 20N 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_Sept_20N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_Sept_20N_FullLon.ps' ; Title12 = ' Temperature Deviation from Zonal Mean September 20N 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_Sept_20N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_Sept_20N_FullLon.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_Sept_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_Sept_20N_FullLon.ps' ; Title18 = ' Temperature Deviation from Zonal Mean September 20N 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_Sept_20N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_Sept_20N_FullLon.ps' ; Title21 = ' Temperature Deviation from Zonal Mean September 20N 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_Sept_20N_FullLon.sav' ; ;September Equator Full Longitude Mean ; iLt = 47 ; ; device, /color, /land, filename='idl_tmt_UT00_Sept_Eq_FullLon.ps' ; Title0 = ' Temperature Deviation from Zonal Mean September Equator 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT03_Sept_Eq_FullLon.ps' ; Title3 = ' Temperature Deviation from Zonal Mean September Equator 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT06_Sept_Eq_FullLon.ps' ; Title6 = ' Temperature Deviation from Zonal Mean September Equator 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_Sept_Eq_FullLon.ps' ; Title9 = ' Temperature Deviation from Zonal Mean September Equator 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_Sept_Eq_FullLon.ps' ; Title12 = ' Temperature Deviation from Zonal Mean September Equator 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_Sept_Eq_FullLon.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March Eq 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT18_Sept_Eq_FullLon.ps' ; Title18 = ' Temperature Deviation from Zonal Mean September Equator 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device, /color, /land, filename='idl_tmt_UT21_Sept_Eq_FullLon.ps' ; Title21 = ' Temperature Deviation from Zonal Mean September Equator 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; ; device,/close ; set_plot, 'x' ; ;March 20N ; device, /color, /land, filename='idl_tmt_UT00.ps' ; Title0 = ' Temperature Deviation from Zonal Mean March 20N 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0,iLt], xrange=[-90.0,90.0],yrange=[-60,100],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt,iLt] ; ; device, /color, /land, filename='idl_tmt_UT03.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March 20N 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0,iLt], xrange=[-135.0,45.0],yrange=[-60,100],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt,iLt] ; ; device, /color, /land, filename='idl_tmt_UT06.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March 20N 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0,iLt], xrange=[-180.0,0.0],yrange=[-60,100],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt,iLt] ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March 20N 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0,iLt],xrange=[135.0,315.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt,iLt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_March_20N.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March 20N 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0,iLt],xrange=[90.0,270.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt,iLt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_March_20N.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0,iLt],xrange=[45.0,225.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt,iLt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_March_20N.sav' ; ; device, /color, /land, filename='idl_tmt_UT15.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15' ; plot, lonw, twzonalmeandev[*,iLt,0,iLt], xrange=[0.0,180.0],yrange=[-60,100],xstyle=1,ystyle=1,title=Title15,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt,iLt] ; ; device, /color, /land, filename='idl_tmt_UT18.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March 20N 285km UT18' ; plot, lonw, twzonalmeandev[*,0,iLt], xrange=[0.0,180.0],yrange=[-60,100],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iPlt,iLt] ; ; device, /color, /land, filename='idl_tmt_UT21.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March 20N 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0,iLt], xrange=[-45.0,135.0],yrange=[-60,100],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt,iLt] ; ;Dec 30S ; iLt = 32 ;30 degrees south latitude ; set_plot, 'ps' ; device, /color, /land, filename='idl_tmt_UT00_Dec_30S.ps' ; Title0 = ' Temperature Deviation from Zonal Mean December 30S 285km UT00' ; plot, lonw, twzonalmeandev[*,0], xrange=[-90.0,90.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iPlt] ; twzonalmeandev0 = twzonalmeandev ; twzonalmeandevorig0 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_Dec_30S.sav' ; ; device, /color, /land, filename='idl_tmt_UT03_Dec_30S.ps' ; Title3 = ' Temperature Deviation from Zonal Mean December 30S 285km UT03' ; plot, lonw, twzonalmeandev[*,0], xrange=[-135.0,45.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iPlt] ; twzonalmeandev3 = twzonalmeandev ; twzonalmeandevorig3 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev3, twzonalmeandevorig3, file = 'waccm_mtm_zmd_UT03_Dec_30S.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_Dec_30S.ps' ; Title6 = ' Temperature Deviation from Zonal Mean December 30S 285km UT06' ; plot, lonw, twzonalmeandev[*,0], xrange=[-180.0,0.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iPlt] ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_Dec_30S.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_Dec_30S.ps' ; Title9 = ' Temperature Deviation from Zonal December 30S 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,0],xrange=[135.0,315.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iPlt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_Dec_30S.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_Dec_30S.ps' ; Title12 = ' Temperature Deviation from Zonal Mean December 30S 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,0],xrange=[90.0,270.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iPlt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_Dec_30S.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_Dec_30S.ps' ; Title15 = ' Temperature Deviation from Zonal Mean December 30S 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,0],xrange=[45.0,225.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iPlt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_Dec_30S.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_Dec_30S.ps' ; Title18 = ' Temperature Deviation from Zonal Mean December 30S 285km UT18' ; plot, lonw, twzonalmeandev[*,0], xrange=[0.0,180.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iPlt] ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_Dec_30S.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_Dec_30S.ps' ; Title21 = ' Temperature Deviation from Zonal Mean December 30S 285km UT21' ; plot, lonw, twzonalmeandev[*,0], xrange=[-45.0,135.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iPlt] ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_Dec_30S.sav' ; ; Dec 30S Full Longitude range mean ; ; set_plot, 'ps' ; iLt = 32 ;30 degrees south latitude ; ; device, /color, /land, filename='idl_tmt_UT00_Dec_30S_FullLon.ps' ; Title0 = ' Temperature Deviation from Zonal Mean December 30S 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev0 = twzonalmeandev ; twzonalmeandevorig0 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_Dec_30S_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT03_Dec_30S_FullLon.ps' ; Title3 = ' Temperature Deviation from Zonal Mean December 30S 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev3 = twzonalmeandev ; twzonalmeandevorig3 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev3, twzonalmeandevorig3, file = 'waccm_mtm_zmd_UT03_Dec_30S_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_Dec_30S_FullLon.ps' ; Title6 = ' Temperature Deviation from Zonal Mean December 30S 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_Dec_30S_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_Dec_30S_FullLon.ps' ; Title9 = ' Temperature Deviation from Zonal December 30S 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_Dec_30S_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_Dec_30S_FullLon.ps' ; Title12 = ' Temperature Deviation from Zonal Mean December 30S 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_Dec_30S_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_Dec_30S_FullLon.ps' ; Title15 = ' Temperature Deviation from Zonal Mean December 30S 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_Dec_30S_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_Dec_30S_FullLon.ps' ; Title18 = ' Temperature Deviation from Zonal Mean December 30S 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_Dec_30S_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_Dec_30S_FullLon.ps' ; Title21 = ' Temperature Deviation from Zonal Mean December 30S 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-120,80],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_Dec_30S_FullLon.sav' ; ;March 20N Gravity Wave parameterization High Resolution ; device, /color, /land, filename='idl_tmt_UT00_MarchGWHR_20N.ps' ; Title0 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev0 = twzonalmeandev ; twzonalmeandevorig0 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_MarchGWHR_20N.sav' ; ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_20N.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev3 = twzonalmeandev ; twzonalmeandevorig3 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev3, twzonalmeandevorig3, file = 'waccm_mtm_zmd_UT03_MarchGWHR_20N.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_20N.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_MarchGWHR_20N.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_20N.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_MarchGWHR_20N.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_20N.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_MarchGWHR_20N.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_20N.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 20N 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_MarchGWHR_20N.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_20N.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_MarchGWHR_20N.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_20N.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March GWHR 20N 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-80,100],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_MarchGWHR_20N.sav' ; ;March 31N Gravity Wave parameterization High Resolution Full Longitude Mean ; iLt = 64 ;31 degrees north latitude ; device, /color, /land, filename='idl_tmt_UT00_MarchGWHR_31N_FullLon.ps' ; Title0 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT00' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-90.0,90.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title0,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev0 = twzonalmeandev ; twzonalmeandevorig0 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev0, twzonalmeandevorig0, file = 'waccm_mtm_zmd_UT00_MarchGWHR_31N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT03_MarchGWHR_31N_FullLon.ps' ; Title3 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT03' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-135.0,45.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title3,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev3 = twzonalmeandev ; twzonalmeandevorig3 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev3, twzonalmeandevorig3, file = 'waccm_mtm_zmd_UT03_MarchGWHR_31N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT06_MarchGWHR_31N_FullLon.ps' ; Title6 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT06' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-180.0,0.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title6,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev6 = twzonalmeandev ; twzonalmeandevorig6 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev6, twzonalmeandevorig6, file = 'waccm_mtm_zmd_UT06_MarchGWHR_31N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT09_MarchGWHR_31N_FullLon.ps' ; Title9 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT09' ; plot, lonworig, twzonalmeandevorig[*,iLt,0],xrange=[135.0,315.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title9,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev9 = twzonalmeandev ; twzonalmeandevorig9 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev9, twzonalmeandevorig9, file = 'waccm_mtm_zmd_UT09_MarchGWHR_31N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT12_MarchGWHR_31N_FullLon.ps' ; Title12 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT12' ; plot, lonworig,twzonalmeandevorig[*,iLt,0],xrange=[90.0,270.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title12,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev12 = twzonalmeandev ; twzonalmeandevorig12 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev12, twzonalmeandevorig12, file = 'waccm_mtm_zmd_UT12_MarchGWHR_31N_FullLon.sav' ; ; xTitleOrig = ' Longitude (Degrees East)' ; device, /color, /land, filename='idl_tmt_UT15_MarchGWHR_31N_FullLon.ps' ; Title15 = ' Temperature Deviation from Zonal Mean March 31N 285km UT15' ; plot,lonworig,twzonalmeandevorig[*,iLt,0],xrange=[45.0,225.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title15,xtitle=xTitleorig,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonworig, twzonalmeandevorig[*,iLt,iPlt] ; twzonalmeandev15 = twzonalmeandev ; twzonalmeandevorig15 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev15, twzonalmeandevorig15, file = 'waccm_mtm_zmd_UT15_MarchGWHR_31N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT18_MarchGWHR_31N_FullLon.ps' ; Title18 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT18' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[0.0,180.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title18,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev18 = twzonalmeandev ; twzonalmeandevorig18 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev18, twzonalmeandevorig18, file = 'waccm_mtm_zmd_UT18_MarchGWHR_31N_FullLon.sav' ; ; device, /color, /land, filename='idl_tmt_UT21_MarchGWHR_31N_FullLon.ps' ; Title21 = ' Temperature Deviation from Zonal Mean March GWHR 31N 285km UT21' ; plot, lonw, twzonalmeandev[*,iLt,0], xrange=[-45.0,135.0],yrange=[-100,80],xstyle=1,ystyle=1,title=Title21,xtitle=xTitle,ytitle=yTitle ; for iPlt = 1,nTAllW-1 DO oplot, lonw, twzonalmeandev[*,iLt,iPlt] ; twzonalmeandev21 = twzonalmeandev ; twzonalmeandevorig21 = twzonalmeandevorig ; save, lonw, lonworig, twzonalmeandev21, twzonalmeandevorig21, file = 'waccm_mtm_zmd_UT21_MarchGWHR_31N_FullLon.sav' ; ; device,/close ; set_plot, 'x' END ;TopTitleA = 'FFT Amplitude of Geopotential Height Lat=10 Level=10,Time=' ;PLOT, /ylog, zAmp[0,*,*,0],xrange=[0,10],xtitle=XTitle,ytitle=YTitle,Title=TopTitleA ; DEVICE, DECOMPOSED = 0 ; LOADCT, 39 ; contour, zamp[1:10,*,20,0],findgen(10),latw,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[-90,90],ystyle=1 ;contour, tamp[1:10,*,20,0],findgen(10),latw,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[-90,90],ystyle=1 set_plot, 'ps' device, /color, filename='idl_tamp.ps' loadct, 39 !p.multi=[2,2] fgen90 = findgen(90) fgen40 = findgen(40) xTitle = 'Time(3hr)' yTitle = 'Frequency' Title = ' Lat=-88 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,1,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-80 Level=100 hPa' contour, reform(transpose(tamp[1:40,5,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-70 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,11,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-60 Level=100 hPa' contour, reform(transpose(tamp[1:40,16,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-50 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,21,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-40 Level=100 hPa' contour, reform(transpose(tamp[1:40,26,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-30 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,32,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-20 Level=100 hPa' contour, reform(transpose(tamp[1:40,37,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-10 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,42,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=0 Level=100 hPa' contour, reform(transpose(tamp[1:40,47,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=10 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,53,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=20 Level=100 hPa' contour, reform(transpose(tamp[1:40,58,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=30 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,63,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=40 Level=100 hPa' contour, reform(transpose(tamp[1:40,69,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=50 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,74,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=60 Level=100 hPa' contour, reform(transpose(tamp[1:40,79,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=70 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,84,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=80 Level=100 hPa' contour, reform(transpose(tamp[1:40,90,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=90 Level=100 hPa FFT Amp. Temperature' contour, reform(transpose(tamp[1:40,94,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle device,/close device, /color, filename='idl_zamp.ps' loadct, 39 !p.multi=[2,2] fgen90 = findgen(90) fgen40 = findgen(40) xTitle = 'Time(3hr)' yTitle = 'Frequency' Title = ' Lat=-88 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,1,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-80 Level=100 hPa' contour, reform(transpose(zamp[1:40,5,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-70 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,11,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-60 Level=100 hPa' contour, reform(transpose(zamp[1:40,16,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-50 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,21,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-40 Level=100 hPa' contour, reform(transpose(zamp[1:40,26,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-30 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,32,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-20 Level=100 hPa' contour, reform(transpose(zamp[1:40,37,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=-10 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,42,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=0 Level=100 hPa' contour, reform(transpose(zamp[1:40,47,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=10 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,53,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=20 Level=100 hPa' contour, reform(transpose(zamp[1:40,58,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=30 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,63,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=40 Level=100 hPa' contour, reform(transpose(zamp[1:40,69,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=50 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,74,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=60 Level=100 hPa' contour, reform(transpose(zamp[1:40,79,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=70 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,84,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=80 Level=100 hPa' contour, reform(transpose(zamp[1:40,90,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle Title = ' Lat=90 Level=100 hPa FFT Amp. Geopot. Ht.' contour, reform(transpose(zamp[1:40,94,1,*])),fgen90,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle device,/close set_plot,'x' ; nTAllW2 = nTAllW/2+1 ; nShiftW = (nXW+1)/2-1 ; zonalfreq,uuNSubW,nTAllW,nXW,nYW,nP,nShiftW,nTAllW2,uuSpcW,uuDTW,uuSDTW,uuP10W,uuP20W,uuP12W,uuP21W ; zonalfreq,vuNSubW,nTAllW,nXW,nYW,nP,nShiftW,nTAllW2,vuSpcW,vuDTW,vuSDTW,vuP10W,vuP20W,vuP12W,vuP21W ; zonalfreq,tuNSubW,nTAllW,nXW,nYW,nP,nShiftW,nTAllW2,tuSpcW,tuDTW,tuSDTW,tuP10W,tuP20W,tuP12W,tuP21W ; xwn = indgen(nx)-nshift ; freq = indgen(nT2) ; xwnw = nshiftw-indgen(nxw) ; freqw = indgen(nTAllW2) ; DEVICE, DECOMPOSED = 0 ; LOADCT, 39 ; nlevels = 20 ; clevels = fltarr(nlevels) ; fmin = min(alog10(uuspc[*,*,36,2])) ; fmax = max(alog10(uuspc[*,*,36,2])) ; conint = (fmax-fmin) / (nlevels-1) ; for i=0,nlevels-1 do clevels(i) = fmin+i*conint ; clevelsw = fltarr(nlevels) ; fminw = min(alog10(uuspcw[*,*,36,2])) ; fmaxw = max(alog10(uuspcw[*,*,36,2])) ; conintw = (fmaxw-fminw) / (nlevels-1) ; for i=0,nlevels-1 do clevelsw(i) = fminw+i*conintw ; window, ret=2 ;contour, alog10(uuspc[*,*,36,0]), xwn,freq, nlev=21, /cell_fill, xrange=[-10,10], yrange=[0,20] ;contour, alog10(uuspc[*,*,36,2]), xwn,freq, nlev=21, /cell_fill, xrange=[-10,10], yrange=[0,10] ;contour, alog10(uuspcw[*,*,36,2]), xwnw,freqw, nlev=21, /cell_fill, xrange=[-10,10], yrange=[0,6] ;plot, lat, uuspc[70,11,*,0] ;plot, lat, uuspc[70,11,*,1] ;plot, lat, uuspc[70,11,*,2] ;contour, alog10(uuspc[*,*,36,2]), xwn,freq,/noerase,/follow,/fill,xstyle=5,ystyle=5,levels=clevels,c_colors=BYTSCL(indgen(20)),xrange=[-10,10], yrange=[0,10] ;contour, alog10(uuspc[*,*,36,2]), xwn,freq,/noerase,/follow,levels=clevels,xstyle=1, xrange=[-10,10], yrange=[0,10] ;contour, alog10(uuspcw[*,*,36,2]), xwnw,freqw,/noerase,/follow,/fill,xstyle=5,ystyle=5,levels=clevelsw,c_colors=BYTSCL(indgen(20)),xrange=[-10,10], yrange=[0,10] ;contour, alog10(uuspcw[*,*,36,2]), xwnw,freqw,/noerase,/follow,levels=clevelsw,xstyle=1, xrange=[-10,10], yrange=[0,10] ; ;set_plot, 'ps' ;device, /color ;loadct, 39 ;!p.multi=[2,2] ;contour, alog10(uuspc[*,*,36,0]), xwn,freq, nlev=21, /cell_fill, xrange=[-10,10],Title='ECMWF 1 hPA' ;contour, alog10(uuspcw[*,*,47,0]), xwnw,freqw, nlev=21, /cell_fill, xrange=[-10,10], yrange=[0,20],Title='WACCM 1 hPA' ;contour, alog10(uuspc[*,*,36,1]), xwn,freq, nlev=21, /cell_fill, xrange=[-10,10],Title='ECMWF 10 hPA' ;contour, alog10(uuspcw[*,*,47,1]), xwnw,freqw, nlev=21, /cell_fill, xrange=[-10,10], yrange=[0,20],Title='WACCM 10 hPA' ;contour, alog10(uuspc[*,*,36,2]), xwn,freq, nlev=21, /cell_fill, xrange=[-10,10],Title='ECMWF 100 hPA' ;contour, alog10(uuspcw[*,*,47,2]), xwnw,freqw, nlev=21, /cell_fill, xrange=[-10,10], yrange=[0,20],Title='WACCM 10 hPA' ;device,/close END