print, ' '

IdJun = NCDF_OPEN('/net/aim/d/waccm/sminyear/06/wa319_2x_sminyear.cam2.h1.1951-06-Mean.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 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_JunMean_UT00-21.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 UT03'
CONTOUR, REFORM(TZSUB[18:90,16:79,1]),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,1]),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 UT06'
CONTOUR, REFORM(TZSUB[0:72,16:79,2]),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,2]),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 UT09'
CONTOUR, REFORM(TZSUBOrig[54:126,16:79,3]),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,3]),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 UT12'
CONTOUR, REFORM(TZSUBOrig[36:108,16:79,4]),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,4]),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 UT15'
CONTOUR, REFORM(TZSUBOrig[18:90,16:79,5]),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,5]),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 UT18'
CONTOUR, REFORM(TZSUBOrig[0:72,16:79,6]),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,6]),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 UT21'
CONTOUR, REFORM(TZSUB[54:126,16:79,7]),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,7]),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'

device,/close
set_plot, 'x'

END