print, ' '

IdWAMDec = NCDF_OPEN('/net/aim/d/waccm/wa319_TIPHYS/WAM/wam_t_jun.nc', /NOWRITE)

;; Coordinates

ncdf_varget, IdWAMDec, 'lon', lonWAM         ; longitude
ncdf_varget, IdWAMDec, 'lat', latWAM         ; latitude
ncdf_varget, IdWAMDec, 'lev', zplevWAM          ; log pressure, -ln(5.e-7/p) (p in hPa)
ncdf_varget, IdWAMDec, 'time', time_minWAM      ; time in minutes

lonWAMOrig = lonWAM
lonWAM = lonWAMOrig - 180.0

;; Field variables

;ncdf_varget, Id, 'Z3', zg     ;geopotential height
;zgm = zg * (1. + zg/6370000.0) /1000.

nxWAM = n_elements(lonWAM)
nyWAM = n_elements(latWAM)
nzWAM = n_elements(zplevWAM)
nTWAM = n_elements(time_minWAM)

iZW = 136

; Dimensions lon,lat,lev,time

ncdf_varget, IdWAMDec, 'temp', tNWAM      ;neutral temperature

tZSubWAM = FLTARR(nXWAM,nYWAM,nTWAM)
tZSubWAM(*,*,*) = tNWAM(*,*,iZW,*)

tZSubWAMOrig = FLTARR(nXWAM,nYWAM,nTWAM)

nLonTemp = nXWAM+nXWAM/2
tZSubWTemp = FLTARR(nLonTemp,nYWAM,nTWAM)

FOR iLt=0,nYWAM-1 DO BEGIN

  tZSubWAMOrig[*,iLt,*] = tZSubWAM[*,iLt,*]
  tZSubWTemp[0:nXWAM-1,iLT,*] = tZSubWAM[0:nXWAM-1,iLt,*]
  tZSubWTemp[nXWAM:nLonTemp-1,iLt,*] = tZSubWAM[0:nXWAM/2-1,iLt,*]
  tZSubWAM[0:nXWAM-1,iLt,*] = tZSubWTemp[nXWAM/2:nLonTemp-1,iLt,*]

ENDFOR

xTitle = 'Longitude (Degrees)'
yTitle = 'Latitude (Degrees)'
LTimes = ['18','20','22','0','2','4','6']
!P.CHARSIZE = 1.2
!P.CHARTHICK = 1.5

set_plot, 'ps'
device, /color, /land, filename='wam_mtm_latlonT_JunMean_UT00-21.ps'
loadct, 39

Title = ' Temperature June Mean 7E-08 hPa UT00'
CONTOUR, REFORM(tZSubWAM[48:144,15:78,0]),lonwam[48:144], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-90.0,90.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[48:144,15:78,0]),lonwam[48:144], latwam[15:78], nlev=24,/follow,xrange=[-90.0,90.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT01'
CONTOUR, REFORM(tZSubWAM[40:136,15:78,1]),lonwam[40:136], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-105.0,75.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[40:136,15:78,1]),lonwam[40:136], latwam[15:78], nlev=24,/follow,xrange=[-105.0,75.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT02'
CONTOUR, REFORM(tZSubWAM[32:128,15:78,2]),lonwam[32:128], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-120.0,60.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[32:128,15:78,2]),lonwam[32:128], latwam[15:78], nlev=24,/follow,xrange=[-120.0,60.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT03'
CONTOUR, REFORM(TZSUBWAM[24:120,15:78,3]),lonwam[24:120], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-135.0,45.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAM[24:120,15:78,3]),lonwam[24:120], latwam[15:78], nlev=24,/follow,xrange=[-135.0,45.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT04'
CONTOUR, REFORM(tZSubWAM[16:112,15:78,4]),lonwam[16:112], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-150.0,30.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[16:112,15:78,4]),lonwam[16:112], latwam[15:78], nlev=24,/follow,xrange=[-150.0,30.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT05'
CONTOUR, REFORM(tZSubWAM[8:104,15:78,5]),lonwam[8:104], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-165.0,15.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[8:104,15:78,5]),lonwam[8:104], latwam[15:78], nlev=24,/follow,xrange=[-165.0,15.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT06'
CONTOUR, REFORM(TZSUBWAM[0:96,15:78,6]),lonwam[0:96], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-180.0,0.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAM[0:96,15:78,6]),lonwam[0:96], latwam[15:78], nlev=24,/follow,xrange=[-180.0,0.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT07'
CONTOUR, REFORM(tZSubWAMOrig[88:184,15:78,7]),lonwamorig[88:184], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[165.0,345.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[88:184,15:78,7]),lonwamorig[88:184], latwam[15:78], nlev=24,/follow,xrange=[165.0,345.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT08'
CONTOUR, REFORM(tZSubWAMOrig[80:176,15:78,8]),lonwamorig[80:176], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[150.0,330.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[80:176,15:78,8]),lonwamorig[80:176], latwam[15:78], nlev=24,/follow,xrange=[-90.0,90.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT09'
CONTOUR, REFORM(TZSUBWAMOrig[72:168,15:78,9]),lonwamorig[72:168], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[135.0,315.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAMOrig[72:168,15:78,9]),lonwamorig[72:168], latwam[15:78],nlev=24,/follow,xrange=[135.0,315.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT10'
CONTOUR, REFORM(tZSubWAMOrig[64:160,15:78,10]),lonwamorig[64:160], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[120.0,300.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[64:160,15:78,10]),lonwamorig[64:160], latwam[15:78], nlev=24,/follow,xrange=[120.0,300.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT11'
CONTOUR, REFORM(tZSubWAMOrig[56:152,15:78,11]),lonwamorig[56:152], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[105.0,285.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[56:152,15:78,11]),lonwamorig[56:152], latwam[15:78], nlev=24,/follow,xrange=[105.0,285.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT12'
CONTOUR, REFORM(TZSUBWAMOrig[48:144,15:78,12]),lonwamorig[48:144], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[90.0,270.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAMOrig[48:144,15:78,12]),lonwamorig[48:144], latwam[15:78],nlev=24,/follow,xrange=[90.0,270.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT13'
CONTOUR, REFORM(tZSubWAMOrig[40:136,15:78,13]),lonwamorig[40:136], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[75.0,255.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[40:136,15:78,13]),lonwamorig[40:136], latwam[15:78], nlev=24,/follow,xrange=[75.0,255.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT14'
CONTOUR, REFORM(tZSubWAMOrig[32:128,15:78,14]),lonwamorig[32:128], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[60.0,240.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[32:128,15:78,14]),lonwamorig[32:128], latwam[15:78], nlev=24,/follow,xrange=[60.0,240.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT15'
CONTOUR, REFORM(TZSUBWAMOrig[24:120,15:78,15]),lonwamorig[24:120], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[45.0,225.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAMOrig[24:120,15:78,15]),lonwamorig[24:120], latwam[15:78],nlev=24,/follow,xrange=[45.0,225.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT16'
CONTOUR, REFORM(tZSubWAMOrig[16:112,15:78,16]),lonwamorig[16:112], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[30.0,210.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[16:112,15:78,16]),lonwamorig[16:112], latwam[15:78], nlev=24,/follow,xrange=[30.0,210.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT17'
CONTOUR, REFORM(tZSubWAMOrig[8:104,15:78,17]),lonwamorig[8:104], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[15.0,195.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAMOrig[8:104,15:78,17]),lonwamorig[8:104], latwam[15:78], nlev=24,/follow,xrange=[15.0,195.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT18'
CONTOUR, REFORM(TZSUBWAMOrig[0:96,15:78,18]),lonwamorig[0:96], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[0.0,180.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAMOrig[0:96,15:78,18]),lonwamorig[0:96], latwam[15:78],nlev=24,/follow,xrange=[0.0,180.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT19'
CONTOUR, REFORM(tZSubWAM[88:184,15:78,19]),lonwam[88:184], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-15.0,165.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[88:184,15:78,19]),lonwam[88:184], latwam[15:78], nlev=24,/follow,xrange=[-15.0,165.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT20'
CONTOUR, REFORM(tZSubWAM[80:176,15:78,20]),lonwam[80:176], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-30.0,150.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[80:176,15:78,20]),lonwam[80:176], latwam[15:78], nlev=24,/follow,xrange=[-30.0,150.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT21'
CONTOUR, REFORM(TZSUBWAM[72:168,15:78,21]),lonwam[72:168], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-45.0,135.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBWAM[72:168,15:78,21]),lonwam[72:168], latwam[15:78],nlev=24,/follow,xrange=[-45.0,135.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT022'
CONTOUR, REFORM(tZSubWAM[64:160,15:78,22]),lonwam[64:160], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-60.0,120.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[64:160,15:78,22]),lonwam[64:160], latwam[15:78], nlev=24,/follow,xrange=[-60.0,120.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

Title = ' Temperature June Mean 7E-08 hPa UT23'
CONTOUR, REFORM(tZSubWAM[56:152,15:78,23]),lonwam[56:152], latwam[15:78], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-75.0,105.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSubWAM[56:152,15:78,23]),lonwam[56:152], latwam[15:78], nlev=24,/follow,xrange=[-75.0,105.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,/overplot,c_charsize=1.2
axis, xaxis=1,/save, xticks=6, xtickn=LTimes,xtitle='Local Time'

device,/close
set_plot, 'x'

END