print, 'Running waccm_waveamp'

;  waccmFilename = "/ptmp/joemci/WACCMOut/wa319_4x5_std.cam2.h1.1999-02-25-00000.nc"
;  waccmFilename = "/waccm/liuh/ensmhr/djf_ensm_03.cam2.h1.1997-12-21-00000.nc"

  nFilesW = 90

  waccmFiles = ["/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-01-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-02-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-03-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-04-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-05-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-06-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-07-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-08-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-09-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-10-00000.nc", $ 
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-11-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-12-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-13-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-14-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-15-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-16-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-17-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-18-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-19-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-20-00000.nc", $ 
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-21-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-22-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-23-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-24-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-25-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-26-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-27-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-28-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-29-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-30-00000.nc", $
                "/aim/d/waccm/smedyear/12/wa319_2x_smedyear.cam2.h1.1951-12-31-00000.nc", $
		"/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-01-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-02-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-03-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-04-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-05-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-06-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-07-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-08-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-09-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-10-00000.nc", $ 
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-11-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-12-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-13-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-14-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-15-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-16-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-17-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-18-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-19-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-20-00000.nc", $ 
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-21-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-22-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-23-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-24-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-25-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-26-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-27-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-28-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-29-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-30-00000.nc", $
                "/aim/d/waccm/smedyear/01/wa319_2x_smedyear.cam2.h1.1951-01-31-00000.nc", $
		"/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-01-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-02-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-03-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-04-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-05-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-06-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-07-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-08-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-09-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-10-00000.nc", $ 
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-11-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-12-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-13-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-14-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-15-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-16-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-17-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-18-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-19-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-20-00000.nc", $ 
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-21-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-22-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-23-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-24-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-25-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-26-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-27-00000.nc", $
                "/aim/d/waccm/smedyear/02/wa319_2x_smedyear.cam2.h1.1951-02-28-00000.nc"]

;
;  Read first WACCM file to get dimensions of data time read
;					       
;  nFirstDayW = 26
  waccmFilename = waccmFiles[0]
  ncdf_waccm_read,waccmFilename,lonW,latW,lev3D,timeOutW,nXW,nYW,nZW,nTW,nameCnst,$
                                     GEOPTH=zW,XWind=uNW,YWind=vNW,ZWind=wNW,TEMP=tNW
  nTInAllW = nTW * nFilesW
;  nFirstDayW = 26
  nTAllW = 0
;
; Loop over input WACCM files to accumlate in time
;  
  tNAllW = FLTARR(nXW,nYW,nZW,nTInAllW)
  uNAllW = FLTARR(nXW,nYW,nZW,nTInAllW)
  vNAllW = FLTARR(nXW,nYW,nZW,nTInAllW)
  zAllW = FLTARR(nXW,nYW,nZW,nTInAllW)
  
  FOR iFW = 0, nFilesW-1 DO BEGIN

;    fileN = iFW + nFirstDayW
;    fileNC = STRTRIM(fileN)
    waccmFilename = waccmFiles[iFW]

    ncdf_waccm_read,waccmFilename,lonW,latW,lev3D,timeOutW,nXW,nYW,nZW,nTW,nameCnst,$
                                     GEOPTH=zW,XWind=uNW,YWind=vNW,ZWind=wNW,TEMP=tNW

    nTInW = 0
    FOR iTW = nTAllW,nTAllW+nTW-1 DO BEGIN

      tNAllW[*,*,*,iTW] = tNW[*,*,*,nTInW] 
      uNAllW[*,*,*,iTW] = uNW[*,*,*,nTInW] 
      vNAllW[*,*,*,iTW] = vNW[*,*,*,nTInW] 
      zAllW[*,*,*,iTW]  = zW[*,*,*,nTInW] 
      nTInW = nTInW + 1

    ENDFOR

    nTAllW = nTAllW + nTW
    
  ENDFOR

;lon,lat,level,time zAmp
 
  zAmp    = FLTARR(nXW,nYW,nZW,nTAllW)
  tAmp    = FLTARR(nXW,nYW,nZW,nTAllW)
  
  zWLonMean    = FLTARR(nYW,nZW,nTAllW)
  zWLatMean    = FLTARR(nZW,nTAllW)
;
;  Find pressure levels to process in WACCM data
;
  FOR iT=0,nTAllW-1 DO BEGIN

    FOR iLt=0,nYW-1 DO BEGIN
;
;  For WACCM data need to find data for selected vertical levels using 3D vertical array
;
      FOR iZ=0,nZW-1 DO BEGIN

        zAmp[*,iLt,iZ,iT] = ABS(FFT(zAllW[*,iLt,iZ,iT]/1000.0, -1))
        tAmp[*,iLt,iZ,iT] = ABS(FFT(tNAllW[*,iLt,iZ,iT], -1))
	zWLonMean[iLt,iZ,iT] = MEAN(zAllW[*,iLt,iZ,iT])
	
      ENDFOR

    ENDFOR

    FOR iZZ=0,nZW-1 DO zWLatMean[iZZ,iT] = MEAN(zWLonMean[*,iZZ,iT])

  ENDFOR

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
loadct, 39
!p.multi=[2,2]
fgen80 = findgen(80)
fgen40 = findgen(40)
xTitle = 'Time(3hr)'
yTitle = 'Frequency'

Title = ' Lat=-88 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,1,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-80 Level=20 km'
contour, reform(transpose(tamp[1:40,5,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-70 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,11,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-60 Level=20 km'
contour, reform(transpose(tamp[1:40,16,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-50 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,21,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-40 Level=20 km'
contour, reform(transpose(tamp[1:40,26,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-30 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,32,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-20 Level=20 km'
contour, reform(transpose(tamp[1:40,37,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=-10 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,42,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=0 Level=20 km'
contour, reform(transpose(tamp[1:40,47,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=10 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,53,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=20 Level=20 km'
contour, reform(transpose(tamp[1:40,58,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=30 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,63,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=40 Level=20 km'
contour, reform(transpose(tamp[1:40,69,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=50 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,74,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=60 Level=20 km'
contour, reform(transpose(tamp[1:40,79,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=70 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,84,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=80 Level=20 km'
contour, reform(transpose(tamp[1:40,90,58,*])),fgen80,fgen40,nlev=140,c_colors=BYTSCL(indgen(140)),/fill,yrange=[0,10],ystyle=1,title=Title, xtitle=xTitle, ytitle=yTitle
Title = ' Lat=90 Level=20 km FFT Amp. Temperature'
contour, reform(transpose(tamp[1:40,94,58,*])),fgen80,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