print, ' '

IdJun = NCDF_OPEN('/net/aim/d/waccm/wa319_TIPHYS/WACCM/WA319JanSMinHourly_2x.cam2.h1.1993-06_MeanTZ3.nc', /NOWRITE)

;; Coordinates

ncdf_varget, IdJun, 'lon', lon         ; longitude
ncdf_varget, IdJun, 'lat', lat         ; latitude
ncdf_varget, IdJun, 'lev', zplev          ; log pressure, -ln(5.e-7/p) (p in hPa)
ncdf_varget, IdJun, 'time', time_min      ; time in minutes

lonOrig = lon
lon = lonOrig - 180.0

;; Field variables

; Dimensions lon,lat,lev,time

ncdf_varget, IdJun, 'T', tN      ;neutral temperature

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

nx = n_elements(lon)
ny = n_elements(lat)
nz = n_elements(zplev)
nT = n_elements(time_min)

tZSub = FLTARR(nX,nY,nT)
;
;  Find altitude index for this time
;    
FOR iLn=0,nX-1 DO BEGIN

  FOR iLt=0,nY-1 DO BEGIN

    FOR iT = 0, nT-1 DO BEGIN

      zWLnLtT = zgm[iLn,iLt,*,iT]
      zInd = WHERE(zWLnLtT GE 280.0 AND zWLnLtT LE 315.0)
;      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

      tZSub(iLn,iLt,iT) = tN(iLn,iLt,iZ,iT)

    ENDFOR

  ENDFOR

ENDFOR

tZSubOrig = FLTARR(nX,nY,nT)

nLonTemp = nX+nX/2
tZSubWTemp = FLTARR(nLonTemp,nY,nT)

FOR iLt=0,nY-1 DO BEGIN

  tZSubOrig[*,iLt,*] = tZSub[*,iLt,*]
  tZSubWTemp[0:nX-1,iLT,*] = tZSub[0:nX-1,iLt,*]
  tZSubWTemp[nX:nLonTemp-1,iLt,*] = tZSub[0:nX/2-1,iLt,*]
  tZSub[0:nX-1,iLt,*] = tZSubWTemp[nX/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='waccm_mtm_latlonT_JunMean1Hrly_UT00-23_295km.ps'
loadct, 39

Title = ' WACCM Temperature June Mean 285 km UT00'
CONTOUR, REFORM(tZSub[36:108,16:79,0]),lon[36:108], lat[16:79], 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(tZSub[36:108,16:79,0]),lon[36:108], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT01'
CONTOUR, REFORM(tZSub[30:102,16:79,1]),lon[30:102], lat[16:79], 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(tZSub[30:102,16:79,1]),lon[30:102], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT02'
CONTOUR, REFORM(tZSub[24:96,16:79,2]),lon[24:96], lat[16:79], 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(tZSub[24:96,16:79,2]),lon[24:96], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT03'
CONTOUR, REFORM(TZSUB[18:90,16:79,3]),lon[18:90], lat[16:79], 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(TZSUB[18:90,16:79,3]),lon[18:90], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT04'
CONTOUR, REFORM(TZSUB[12:84,16:79,4]),lon[12:84], lat[16:79], 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(TZSUB[12:84,16:79,4]),lon[12:84], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT05'
CONTOUR, REFORM(TZSUB[6:78,16:79,5]),lon[6:78], lat[16:79], 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(TZSUB[6:78,16:79,5]),lon[6:78], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT06'
CONTOUR, REFORM(TZSUB[0:72,16:79,6]),lon[0:72], lat[16:79], 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(TZSUB[0:72,16:79,6]),lon[0:72], lat[16:79], 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 = ' WACCM Temperature June Mean 285 km UT07'
CONTOUR, REFORM(TZSUBOrig[66:138,16:79,7]),lonorig[66:138], lat[16:79], 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(TZSUBOrig[66:138,16:79,7]),lonorig[66:138], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT08'
CONTOUR, REFORM(TZSUBOrig[60:132,16:79,8]),lonorig[60:132], lat[16:79], 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(TZSUBOrig[60:132,16:79,8]),lonorig[60:132], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT09'
CONTOUR, REFORM(TZSUBOrig[54:126,16:79,9]),lonorig[54:126], lat[16:79], 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(TZSUBOrig[54:126,16:79,9]),lonorig[54:126], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT10'
CONTOUR, REFORM(TZSUBOrig[48:120,16:79,10]),lonorig[48:120], lat[16:79], 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(TZSUBOrig[48:120,16:79,10]),lonorig[48:120], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT11'
CONTOUR, REFORM(TZSUBOrig[42:114,16:79,11]),lonorig[42:114], lat[16:79], nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[105.0,315.0], yrange=[-60.0,60.0], xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[42:114,16:79,11]),lonorig[42:114], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT12'
CONTOUR, REFORM(TZSUBOrig[36:108,16:79,12]),lonorig[36:108], lat[16:79], 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(TZSUBOrig[36:108,16:79,12]),lonorig[36:108], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT13'
CONTOUR, REFORM(TZSUBOrig[30:102,16:79,13]),lonorig[30:102], lat[16:79], 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(TZSUBOrig[30:102,16:79,13]),lonorig[30:102], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT14'
CONTOUR, REFORM(TZSUBOrig[24:96,16:79,14]),lonorig[24:96], lat[16:79], 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(TZSUBOrig[24:96,16:79,14]),lonorig[24:96], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT15'
CONTOUR, REFORM(TZSUBOrig[18:90,16:79,15]),lonorig[18:90], lat[16:79], 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(TZSUBOrig[18:90,16:79,15]),lonorig[18:90], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT16'
CONTOUR, REFORM(TZSUBOrig[12:84,16:79,16]),lonorig[12:84], lat[16:79], 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(TZSUBOrig[12:84,16:79,16]),lonorig[12:84], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT17'
CONTOUR, REFORM(TZSUBOrig[6:78,16:79,17]),lonorig[6:78], lat[16:79], 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(TZSUBOrig[6:78,16:79,17]),lonorig[6:78], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT18'
CONTOUR, REFORM(TZSUBOrig[0:72,16:79,18]),lonorig[0:72], lat[16:79], 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(TZSUBOrig[0:72,16:79,18]),lonorig[0:72], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT19'
CONTOUR, REFORM(TZSUB[66:138,16:79,19]),lon[66:138], lat[16:79], 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(TZSUB[66:138,16:79,19]),lon[66:138], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT20'
CONTOUR, REFORM(TZSUB[60:132,16:79,20]),lon[60:132], lat[16:79], 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(TZSUB[60:132,16:79,20]),lon[60:132], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT21'
CONTOUR, REFORM(TZSUB[54:126,16:79,21]),lon[54:126], lat[16:79], 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(TZSUB[54:126,16:79,21]),lon[54:126], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT22'
CONTOUR, REFORM(TZSUB[48:120,16:79,22]),lon[48:120], lat[16:79], 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(TZSUB[48:120,16:79,22]),lon[48:120], lat[16:79],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 = ' WACCM Temperature June Mean 285 km UT23'
CONTOUR, REFORM(TZSUB[42:114,16:79,23]),lon[42:114], lat[16:79], 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(TZSUB[42:114,16:79,23]),lon[42:114], lat[16:79],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'

device,/close
set_plot, 'x'

END