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,nZ,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 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

iLt = 58 ;20N	
;iLt = 56 ;15N	
;iLt = 53 ;10N

tZSub[*,*,*] = tN(*,iLt,*,*)

tZSubOrig = FLTARR(nX,nZ,nT)

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

FOR iZ=0,nZ-1 DO BEGIN

  tZSubOrig[*,iZ,*] = tZSub[*,iZ,*]
  tZSubWTemp[0:nX-1,iZ,*] = tZSub[0:nX-1,iZ,*]
  tZSubWTemp[nX:nLonTemp-1,iZ,*] = tZSub[0:nX/2-1,iZ,*]
  tZSub[0:nX-1,iZ,*] = tZSubWTemp[nX/2:nLonTemp-1,iZ,*]

ENDFOR

xTitle = 'Longitude (Degrees)'
yTitle = 'Pressure (hPa)'
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_latlevT_JunMean1Hrly_UT00-23_20N.ps'
loadct, 39

Title = ' WACCM Temperature June Mean 285 km UT00'
CONTOUR, REFORM(tZSub[36:108,*,0]),lon[36:108],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)),xrange=[-90.0,90.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSub[36:108,*,0]),lon[36:108],zplev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,1]),lon[30:102],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-105.0,75.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSub[30:102,*,1]),lon[30:102],zplev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,2]),lon[24:96],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-120.0,60.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(tZSub[24:96,*,2]),lon[24:96],zplev, nlev=24,/follow,xrange=[-90.0,90.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,3]),lon[18:90],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-135.0,45.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[18:90,*,3]),lon[18:90],zplev, nlev=24,/follow,xrange=[-135.0,45.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,4]),lon[12:84],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-150.0,30.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[12:84,*,4]),lon[12:84],zplev, nlev=24,/follow,xrange=[-135.0,45.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,5]),lon[6:78],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-165.0,15.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[6:78,*,5]),lon[6:78],zplev, nlev=24,/follow,xrange=[-135.0,45.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,6]),lon[0:72],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-180.0,0.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[0:72,*,6]),lon[0:72],zplev, nlev=24,/follow,xrange=[-180.0,0.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,7]),lonorig[66:138],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[165.0,345.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[66:138,*,7]),lonorig[66:138],zplev,nlev=24,/follow,xrange=[135.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,8]),lonorig[60:132],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[150.0,330.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[60:132,*,8]),lonorig[60:132],zplev,nlev=24,/follow,xrange=[135.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,9]),lonorig[54:126],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[135.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[54:126,*,9]),lonorig[54:126],zplev,nlev=24,/follow,xrange=[135.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,10]),lonorig[48:120],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[120.0,300.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[48:120,*,10]),lonorig[48:120],zplev,nlev=24,/follow,xrange=[135.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,11]),lonorig[42:114],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[105.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[42:114,*,11]),lonorig[42:114],zplev,nlev=24,/follow,xrange=[135.0,315.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,12]),lonorig[36:108],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[90.0,270.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[36:108,*,12]),lonorig[36:108],zplev,nlev=24,/follow,xrange=[90.0,270.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,13]),lonorig[30:102],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[75.0,255.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[30:102,*,13]),lonorig[30:102],zplev,nlev=24,/follow,xrange=[90.0,270.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,14]),lonorig[24:96],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[60.0,240.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[24:96,*,14]),lonorig[24:96],zplev,nlev=24,/follow,xrange=[90.0,270.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,15]),lonorig[18:90],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[45.0,225.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[18:90,*,15]),lonorig[18:90],zplev,nlev=24,/follow,xrange=[45.0,225.0], yrange=[10.0E3,10.0E-9], /ylog,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]),lonorig[12:84],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[30.0,210.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[12:84,*,16]),lonorig[12:84],zplev,nlev=24,/follow,xrange=[45.0,225.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,17]),lonorig[6:78],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[15.0,195.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[6:78,*,17]),lonorig[6:78],zplev,nlev=24,/follow,xrange=[45.0,225.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,18]),lonorig[0:72],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[0.0,180.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUBOrig[0:72,*,18]),lonorig[0:72],zplev,nlev=24,/follow,xrange=[0.0,180.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,19]),lon[66:138],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-15.0,165.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[66:138,*,19]),lon[66:138],zplev,nlev=24,/follow,xrange=[-45.0,135.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,20]),lon[60:132],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-30.0,150.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[60:132,*,20]),lon[60:132],zplev,nlev=24,/follow,xrange=[-45.0,135.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,21]),lon[54:126],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-45.0,135.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[54:126,*,21]),lon[54:126],zplev,nlev=24,/follow,xrange=[-45.0,135.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,22]),lon[48:120],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-60.0,120.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[48:120,*,22]),lon[48:120],zplev,nlev=24,/follow,xrange=[-45.0,135.0], yrange=[10.0E3,10.0E-9], /ylog,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,*,23]),lon[42:114],zplev, nlev=24, /fill, c_colors=BYTSCL(indgen(24)), xrange=[-75.0,105.0], yrange=[10.0E3,10.0E-9], /ylog,xstyle=9,ystyle=1,subTitle=Title,xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
CONTOUR, REFORM(TZSUB[42:114,*,23]),lon[42:114],zplev,nlev=24,/follow,xrange=[-45.0,135.0], yrange=[10.0E3,10.0E-9], /ylog,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