PRO plot_waccm_wam_tuv_wvn_recon_dec, tnWACCM, tnWAM, unWACCM, unWAM, vnWACCM, vnWAM, zplev, zplevWAM, nz, nzWAM

; This procedure was created on Tuesday, August 31st, 2010.
; This procedure will produce line plots of temperature or wind amplitudes (y-axis) and wavenumber (x-axis)

SET_PLOT, 'PS'
;
;Temperature Reconstructed WNs 1-12
;
;psfile = '/local/i/joemci/SW/WAM/twvn_Dec_Mean_40N_WACCM_WAM_recon.ps'
psfile = 'twvn_Dec_Mean_40N_WACCM_WAM_recon.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmp1WACCMMin = FLTARR(nz)
tmp2WACCMMin = FLTARR(nz)
tmp3WACCMMin = FLTARR(nz)
tmp4WACCMMin = FLTARR(nz)
tmp5WACCMMin = FLTARR(nz)
tmp6WACCMMin = FLTARR(nz)
tmp7WACCMMin = FLTARR(nz)
tmp8WACCMMin = FLTARR(nz)
tmp9WACCMMin = FLTARR(nz)
tmp10WACCMMin = FLTARR(nz)
tmp11WACCMMin = FLTARR(nz)
tmp12WACCMMin = FLTARR(nz)

tmp1WACCMMax = FLTARR(nz)
tmp2WACCMMax = FLTARR(nz)
tmp3WACCMMax = FLTARR(nz)
tmp4WACCMMax = FLTARR(nz)
tmp5WACCMMax = FLTARR(nz)
tmp6WACCMMax = FLTARR(nz)
tmp7WACCMMax = FLTARR(nz)
tmp8WACCMMax = FLTARR(nz)
tmp9WACCMMax = FLTARR(nz)
tmp10WACCMMax = FLTARR(nz)
tmp11WACCMMax = FLTARR(nz)
tmp12WACCMMax = FLTARR(nz)

FOR iz = nz-1,0,-1 DO BEGIN

  newiz = nz - iz - 1

  tmp1WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,1]))
  tmp2WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,2]))
  tmp3WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,3]))
  tmp4WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,4]))
  tmp5WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,5]))
  tmp6WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,6]))
  tmp7WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,7]))
  tmp8WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,8]))
  tmp9WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,9]))
  tmp10WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,10]))
  tmp11WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,11]))
  tmp12WACCMMin(newiz) = REFORM(MIN(tnWACCM[*,69,iz,12]))

  tmp1WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,1]))
  tmp2WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,2]))
  tmp3WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,3]))
  tmp4WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,4]))
  tmp5WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,5]))
  tmp6WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,6]))
  tmp7WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,7]))
  tmp8WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,8]))
  tmp9WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,9]))
  tmp10WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,10]))
  tmp11WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,11]))
  tmp12WACCMMax(newiz) = REFORM(MAX(tnWACCM[*,69,iz,12]))

ENDFOR

zplevT = REVERSE(zplev)

tmp1WACCMMin = interpol(tmp1WACCMMin, zplevT, zplevWAM)
tmp2WACCMMin = interpol(tmp2WACCMMin, zplevT, zplevWAM)
tmp3WACCMMin = interpol(tmp3WACCMMin, zplevT, zplevWAM)
tmp4WACCMMin = interpol(tmp4WACCMMin, zplevT, zplevWAM)
tmp5WACCMMin = interpol(tmp5WACCMMin, zplevT, zplevWAM)
tmp6WACCMMin = interpol(tmp6WACCMMin, zplevT, zplevWAM)
tmp7WACCMMin = interpol(tmp7WACCMMin, zplevT, zplevWAM)
tmp8WACCMMin = interpol(tmp8WACCMMin, zplevT, zplevWAM)
tmp9WACCMMin = interpol(tmp9WACCMMin, zplevT, zplevWAM)
tmp10WACCMMin = interpol(tmp10WACCMMin, zplevT, zplevWAM)
tmp11WACCMMin = interpol(tmp11WACCMMin, zplevT, zplevWAM)
tmp12WACCMMin = interpol(tmp12WACCMMin, zplevT, zplevWAM)

tmp1WACCMMax = interpol(tmp1WACCMMax, zplevT, zplevWAM)
tmp2WACCMMax = interpol(tmp2WACCMMax, zplevT, zplevWAM)
tmp3WACCMMax = interpol(tmp3WACCMMax, zplevT, zplevWAM)
tmp4WACCMMax = interpol(tmp4WACCMMax, zplevT, zplevWAM)
tmp5WACCMMax = interpol(tmp5WACCMMax, zplevT, zplevWAM)
tmp6WACCMMax = interpol(tmp6WACCMMax, zplevT, zplevWAM)
tmp7WACCMMax = interpol(tmp7WACCMMax, zplevT, zplevWAM)
tmp8WACCMMax = interpol(tmp8WACCMMax, zplevT, zplevWAM)
tmp9WACCMMax = interpol(tmp9WACCMMax, zplevT, zplevWAM)
tmp10WACCMMax = interpol(tmp10WACCMMax, zplevT, zplevWAM)
tmp11WACCMMax = interpol(tmp11WACCMMax, zplevT, zplevWAM)
tmp12WACCMMax = interpol(tmp12WACCMMax, zplevT, zplevWAM)

;tmpa1 = ABS(FFT(tmp1,-1))
;tmpa2 = ABS(FFT(tmp2,-1))
;tmpa3 = ABS(FFT(tmp3,-1))
;tmpa4 = ABS(FFT(tmp4,-1))
;tmpa5 = ABS(FFT(tmp5,-1))
;tmpa6 = ABS(FFT(tmp6,-1))
;tmpa7 = ABS(FFT(tmp7,-1))
;tmpa8 = ABS(FFT(tmp8,-1))
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmp1WAMMin = FLTARR(nzWAM)
tmp2WAMMin = FLTARR(nzWAM)
tmp3WAMMin = FLTARR(nzWAM)
tmp4WAMMin = FLTARR(nzWAM)
tmp5WAMMin = FLTARR(nzWAM)
tmp6WAMMin = FLTARR(nzWAM)
tmp7WAMMin = FLTARR(nzWAM)
tmp8WAMMin = FLTARR(nzWAM)
tmp9WAMMin = FLTARR(nzWAM)
tmp10WAMMin = FLTARR(nzWAM)
tmp11WAMMin = FLTARR(nzWAM)
tmp12WAMMin = FLTARR(nzWAM)

tmp1WAMMax = FLTARR(nzWAM)
tmp2WAMMax = FLTARR(nzWAM)
tmp3WAMMax = FLTARR(nzWAM)
tmp4WAMMax = FLTARR(nzWAM)
tmp5WAMMax = FLTARR(nzWAM)
tmp6WAMMax = FLTARR(nzWAM)
tmp7WAMMax = FLTARR(nzWAM)
tmp8WAMMax = FLTARR(nzWAM)
tmp9WAMMax = FLTARR(nzWAM)
tmp10WAMMax = FLTARR(nzWAM)
tmp11WAMMax = FLTARR(nzWAM)
tmp12WAMMax = FLTARR(nzWAM)

FOR iz = 0, nzWAM-1 DO BEGIN

  tmp1WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,1]))
  tmp2WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,2]))
  tmp3WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,3]))
  tmp4WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,4]))
  tmp5WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,5]))
  tmp6WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,6]))
  tmp7WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,7]))
  tmp8WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,8]))
  tmp9WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,9]))
  tmp10WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,10]))
  tmp11WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,11]))
  tmp12WAMMin(iz) = REFORM(MIN(tnWAM[*,68,iz,12]))

  tmp1WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,1]))
  tmp2WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,2]))
  tmp3WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,3]))
  tmp4WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,4]))
  tmp5WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,5]))
  tmp6WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,6]))
  tmp7WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,7]))
  tmp8WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,8]))
  tmp9WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,9]))
  tmp10WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,10]))
  tmp11WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,11]))
  tmp12WAMMax(iz) = REFORM(MAX(tnWAM[*,68,iz,12]))

ENDFOR

plot, tmp1WACCMMax,zplevWAM,xrange=[0.0,70.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 1 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp1WAMMin,zplevWAM,linestyle=1
;oplot, tmp1WACCMMin,zplevWAM
oplot, tmp1WAMMax,zplevWAM,linestyle=1

plot, tmp2WACCMMax,zplevWAM,xrange=[0.0,40.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 2 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp2WAMMin,zplevWAM,linestyle=1
;oplot, tmp2WACCMMin,zplevWAM
oplot, tmp2WAMMax,zplevWAM,linestyle=1

plot, tmp3WACCMMax,zplevWAM,xrange=[0.0,30.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 3 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp3WAMMin,zplevWAM,linestyle=1
;oplot, tmp3WACCMMin,zplevWAM
oplot, tmp3WAMMax,zplevWAM,linestyle=1

plot, tmp4WACCMMax,zplevWAM,xrange=[0.0,15.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 4 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp4WAMMin,zplevWAM,linestyle=1
;oplot, tmp4WACCMMin,zplevWAM
oplot, tmp4WAMMax,zplevWAM,linestyle=1

plot, tmp5WACCMMax,zplevWAM,xrange=[0.0,15.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 5 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp5WAMMin,zplevWAM,linestyle=1
;oplot, tmp5WACCMMin,zplevWAM
oplot, tmp5WAMMax,zplevWAM,linestyle=1

plot, tmp6WACCMMax,zplevWAM,xrange=[0.0,15.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 6 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp6WAMMin,zplevWAM,linestyle=1
;oplot, tmp6WACCMMin,zplevWAM
oplot, tmp6WAMMax,zplevWAM,linestyle=1

plot, tmp7WACCMMax,zplevWAM,xrange=[0.0,10.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 7 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp7WAMMin,zplevWAM,linestyle=1
;oplot, tmp7WACCMMin,zplevWAM
oplot, tmp7WAMMax,zplevWAM,linestyle=1

plot, tmp8WACCMMax,zplevWAM,xrange=[0.0,5.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 8 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp8WAMMin,zplevWAM,linestyle=1
;oplot, tmp8WACCMMin,zplevWAM
oplot, tmp8WAMMax,zplevWAM,linestyle=1

plot, tmp9WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 9 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp9WAMMin,zplevWAM,linestyle=1
;oplot, tmp9WACCMMin,zplevWAM
oplot, tmp9WAMMax,zplevWAM,linestyle=1

plot, tmp10WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 10 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp10WAMMin,zplevWAM,linestyle=1
;oplot, tmp10WACCMMin,zplevWAM
oplot, tmp10WAMMax,zplevWAM,linestyle=1

;plot, tmp11WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
;	  yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 11 UT00 40N',$
;	  ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
;	  charsize = 1.2 
;
;;oplot, tmp11WAMMin,zplevWAM,linestyle=1
;;oplot, tmp11WACCMMin,zplevWAM
;oplot, tmp11WAMMax,zplevWAM,linestyle=1
;
;plot, tmp12WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
;	  yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Temperature December Mean WN 12 UT00 40N',$
;	  ytitle = 'Pressure', xtitle = 'Temperature ',ycharsize = 1.0, xcharsize = 1.2,$
;	  charsize = 1.2 
;
;oplot, tmp12WAMMin,zplevWAM,linestyle=1
;oplot, tmp12WACCMMin,zplevWAM
;oplot, tmp12WAMMax,zplevWAM,linestyle=1

DEVICE, /CLOSE

;
;Zonal Wind Reconstructed WNs 1-12
;
;psfile = '/local/i/joemci/SW/WAM/uwvn_Dec_Mean_40N_WACCM_WAM_recon.ps'
psfile = 'uwvn_Dec_Mean_40N_WACCM_WAM_recon.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmp1WACCMMin = FLTARR(nz)
tmp2WACCMMin = FLTARR(nz)
tmp3WACCMMin = FLTARR(nz)
tmp4WACCMMin = FLTARR(nz)
tmp5WACCMMin = FLTARR(nz)
tmp6WACCMMin = FLTARR(nz)
tmp7WACCMMin = FLTARR(nz)
tmp8WACCMMin = FLTARR(nz)
tmp9WACCMMin = FLTARR(nz)
tmp10WACCMMin = FLTARR(nz)
tmp11WACCMMin = FLTARR(nz)
tmp12WACCMMin = FLTARR(nz)

tmp1WACCMMax = FLTARR(nz)
tmp2WACCMMax = FLTARR(nz)
tmp3WACCMMax = FLTARR(nz)
tmp4WACCMMax = FLTARR(nz)
tmp5WACCMMax = FLTARR(nz)
tmp6WACCMMax = FLTARR(nz)
tmp7WACCMMax = FLTARR(nz)
tmp8WACCMMax = FLTARR(nz)
tmp9WACCMMax = FLTARR(nz)
tmp10WACCMMax = FLTARR(nz)
tmp11WACCMMax = FLTARR(nz)
tmp12WACCMMax = FLTARR(nz)

FOR iz = nz-1,0,-1 DO BEGIN

  newiz = nz - iz - 1

  tmp1WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,1]))
  tmp2WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,2]))
  tmp3WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,3]))
  tmp4WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,4]))
  tmp5WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,5]))
  tmp6WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,6]))
  tmp7WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,7]))
  tmp8WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,8]))
  tmp9WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,9]))
  tmp10WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,10]))
  tmp11WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,11]))
  tmp12WACCMMin(newiz) = REFORM(MIN(unWACCM[*,69,iz,12]))

  tmp1WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,1]))
  tmp2WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,2]))
  tmp3WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,3]))
  tmp4WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,4]))
  tmp5WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,5]))
  tmp6WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,6]))
  tmp7WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,7]))
  tmp8WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,8]))
  tmp9WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,9]))
  tmp10WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,10]))
  tmp11WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,11]))
  tmp12WACCMMax(newiz) = REFORM(MAX(unWACCM[*,69,iz,12]))

ENDFOR

zplevU = REVERSE(zplev)

tmp1WACCMMin = interpol(tmp1WACCMMin, zplevU, zplevWAM)
tmp2WACCMMin = interpol(tmp2WACCMMin, zplevU, zplevWAM)
tmp3WACCMMin = interpol(tmp3WACCMMin, zplevU, zplevWAM)
tmp4WACCMMin = interpol(tmp4WACCMMin, zplevU, zplevWAM)
tmp5WACCMMin = interpol(tmp5WACCMMin, zplevU, zplevWAM)
tmp6WACCMMin = interpol(tmp6WACCMMin, zplevU, zplevWAM)
tmp7WACCMMin = interpol(tmp7WACCMMin, zplevU, zplevWAM)
tmp8WACCMMin = interpol(tmp8WACCMMin, zplevU, zplevWAM)
tmp9WACCMMin = interpol(tmp9WACCMMin, zplevU, zplevWAM)
tmp10WACCMMin = interpol(tmp10WACCMMin, zplevU, zplevWAM)
tmp11WACCMMin = interpol(tmp11WACCMMin, zplevU, zplevWAM)
tmp12WACCMMin = interpol(tmp12WACCMMin, zplevU, zplevWAM)

tmp1WACCMMax = interpol(tmp1WACCMMax, zplevU, zplevWAM)
tmp2WACCMMax = interpol(tmp2WACCMMax, zplevU, zplevWAM)
tmp3WACCMMax = interpol(tmp3WACCMMax, zplevU, zplevWAM)
tmp4WACCMMax = interpol(tmp4WACCMMax, zplevU, zplevWAM)
tmp5WACCMMax = interpol(tmp5WACCMMax, zplevU, zplevWAM)
tmp6WACCMMax = interpol(tmp6WACCMMax, zplevU, zplevWAM)
tmp7WACCMMax = interpol(tmp7WACCMMax, zplevU, zplevWAM)
tmp8WACCMMax = interpol(tmp8WACCMMax, zplevU, zplevWAM)
tmp9WACCMMax = interpol(tmp9WACCMMax, zplevU, zplevWAM)
tmp10WACCMMax = interpol(tmp10WACCMMax, zplevU, zplevWAM)
tmp11WACCMMax = interpol(tmp11WACCMMax, zplevU, zplevWAM)
tmp12WACCMMax = interpol(tmp12WACCMMax, zplevU, zplevWAM)
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmp1WAMMin = FLTARR(nzWAM)
tmp2WAMMin = FLTARR(nzWAM)
tmp3WAMMin = FLTARR(nzWAM)
tmp4WAMMin = FLTARR(nzWAM)
tmp5WAMMin = FLTARR(nzWAM)
tmp6WAMMin = FLTARR(nzWAM)
tmp7WAMMin = FLTARR(nzWAM)
tmp8WAMMin = FLTARR(nzWAM)
tmp9WAMMin = FLTARR(nzWAM)
tmp10WAMMin = FLTARR(nzWAM)
tmp11WAMMin = FLTARR(nzWAM)
tmp12WAMMin = FLTARR(nzWAM)

tmp1WAMMax = FLTARR(nzWAM)
tmp2WAMMax = FLTARR(nzWAM)
tmp3WAMMax = FLTARR(nzWAM)
tmp4WAMMax = FLTARR(nzWAM)
tmp5WAMMax = FLTARR(nzWAM)
tmp6WAMMax = FLTARR(nzWAM)
tmp7WAMMax = FLTARR(nzWAM)
tmp8WAMMax = FLTARR(nzWAM)
tmp9WAMMax = FLTARR(nzWAM)
tmp10WAMMax = FLTARR(nzWAM)
tmp11WAMMax = FLTARR(nzWAM)
tmp12WAMMax = FLTARR(nzWAM)

FOR iz = 0, nzWAM-1 DO BEGIN

  tmp1WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,1]))
  tmp2WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,2]))
  tmp3WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,3]))
  tmp4WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,4]))
  tmp5WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,5]))
  tmp6WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,6]))
  tmp7WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,7]))
  tmp8WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,8]))
  tmp9WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,9]))
  tmp10WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,10]))
  tmp11WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,11]))
  tmp12WAMMin(iz) = REFORM(MIN(unWAM[*,68,iz,12]))

  tmp1WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,1]))
  tmp2WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,2]))
  tmp3WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,3]))
  tmp4WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,4]))
  tmp5WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,5]))
  tmp6WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,6]))
  tmp7WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,7]))
  tmp8WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,8]))
  tmp9WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,9]))
  tmp10WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,10]))
  tmp11WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,11]))
  tmp12WAMMax(iz) = REFORM(MAX(unWAM[*,68,iz,12]))

ENDFOR

plot, tmp1WACCMMax,zplevWAM,xrange=[0.0,80.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 1 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp1WAMMin,zplevWAM,linestyle=1
;oplot, tmp1WACCMMin,zplevWAM
oplot, tmp1WAMMax,zplevWAM,linestyle=1

plot, tmp2WACCMMax,zplevWAM,xrange=[0.0,40.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 2 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp2WAMMin,zplevWAM,linestyle=1
;oplot, tmp2WACCMMin,zplevWAM
oplot, tmp2WAMMax,zplevWAM,linestyle=1

plot, tmp3WACCMMax,zplevWAM,xrange=[0.0,30.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 3 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp3WAMMin,zplevWAM,linestyle=1
;oplot, tmp3WACCMMin,zplevWAM
oplot, tmp3WAMMax,zplevWAM,linestyle=1

plot, tmp4WACCMMax,zplevWAM,xrange=[0.0,20.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 4 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp4WAMMin,zplevWAM,linestyle=1
;oplot, tmp4WACCMMin,zplevWAM
oplot, tmp4WAMMax,zplevWAM,linestyle=1

plot, tmp5WACCMMax,zplevWAM,xrange=[0.0,30.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 5 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp5WAMMin,zplevWAM,linestyle=1
;oplot, tmp5WACCMMin,zplevWAM
oplot, tmp5WAMMax,zplevWAM,linestyle=1

plot, tmp6WACCMMax,zplevWAM,xrange=[0.0,10.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 6 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp6WAMMin,zplevWAM,linestyle=1
;oplot, tmp6WACCMMin,zplevWAM
oplot, tmp6WAMMax,zplevWAM,linestyle=1

plot, tmp7WACCMMax,zplevWAM,xrange=[0.0,10.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 7 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp7WAMMin,zplevWAM,linestyle=1
;oplot, tmp7WACCMMin,zplevWAM
oplot, tmp7WAMMax,zplevWAM,linestyle=1

plot, tmp8WACCMMax,zplevWAM,xrange=[0.0,10.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 8 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp8WAMMin,zplevWAM,linestyle=1
;oplot, tmp8WACCMMin,zplevWAM
oplot, tmp8WAMMax,zplevWAM,linestyle=1

plot, tmp9WACCMMax,zplevWAM,xrange=[0.0,5.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 9 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp9WAMMin,zplevWAM,linestyle=1
;oplot, tmp9WACCMMin,zplevWAM
oplot, tmp9WAMMax,zplevWAM,linestyle=1

plot, tmp10WACCMMax,zplevWAM,xrange=[0.0,4.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 10 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp10WAMMin,zplevWAM,linestyle=1
;oplot, tmp10WACCMMin,zplevWAM
oplot, tmp10WAMMax,zplevWAM,linestyle=1

;plot, tmp11WACCMMax,zplevWAM,xrange=[-3.0,3.0],xstyle=1,$
;	  yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 11 UT00 40N',$
;	  ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
;	  charsize = 1.2 
;
;;oplot, tmp11WAMMin,zplevWAM,linestyle=1
;;oplot, tmp11WACCMMin,zplevWAM
;oplot, tmp11WAMMax,zplevWAM,linestyle=1
;
;plot, tmp12WACCMMax,zplevWAM,xrange=[-4.0,4.0],xstyle=1,$
;	  yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Zonal Wind December Mean WN 12 UT00 40N',$
;	  ytitle = 'Pressure', xtitle = 'Zonal Wind ',ycharsize = 1.0, xcharsize = 1.2,$
;	  charsize = 1.2 
;
;oplot, tmp12WAMMin,zplevWAM,linestyle=1
;oplot, tmp12WACCMMin,zplevWAM
;oplot, tmp12WAMMax,zplevWAM,linestyle=1

DEVICE, /CLOSE
;
;Meridional Wind Reconstructed WNs 1-12
;
;psfile = '/local/i/joemci/SW/WAM/WAM/vwvn_Dec_Mean_40N_WACCM_WAM_recon.ps'
psfile = 'vwvn_Dec_Mean_40N_WACCM_WAM_recon.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile
;
; WACCM is top down so reverse needed
;
tmp1WACCMMin = FLTARR(nz)
tmp2WACCMMin = FLTARR(nz)
tmp3WACCMMin = FLTARR(nz)
tmp4WACCMMin = FLTARR(nz)
tmp5WACCMMin = FLTARR(nz)
tmp6WACCMMin = FLTARR(nz)
tmp7WACCMMin = FLTARR(nz)
tmp8WACCMMin = FLTARR(nz)
tmp9WACCMMin = FLTARR(nz)
tmp10WACCMMin = FLTARR(nz)
tmp11WACCMMin = FLTARR(nz)
tmp12WACCMMin = FLTARR(nz)

tmp1WACCMMax = FLTARR(nz)
tmp2WACCMMax = FLTARR(nz)
tmp3WACCMMax = FLTARR(nz)
tmp4WACCMMax = FLTARR(nz)
tmp5WACCMMax = FLTARR(nz)
tmp6WACCMMax = FLTARR(nz)
tmp7WACCMMax = FLTARR(nz)
tmp8WACCMMax = FLTARR(nz)
tmp9WACCMMax = FLTARR(nz)
tmp10WACCMMax = FLTARR(nz)
tmp11WACCMMax = FLTARR(nz)
tmp12WACCMMax = FLTARR(nz)

FOR iz = nz-1,0,-1 DO BEGIN

  newiz = nz - iz - 1

  tmp1WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,1]))
  tmp2WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,2]))
  tmp3WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,3]))
  tmp4WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,4]))
  tmp5WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,5]))
  tmp6WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,6]))
  tmp7WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,7]))
  tmp8WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,8]))
  tmp9WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,9]))
  tmp10WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,10]))
  tmp11WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,11]))
  tmp12WACCMMin(newiz) = REFORM(MIN(vnWACCM[*,69,iz,12]))

  tmp1WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,1]))
  tmp2WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,2]))
  tmp3WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,3]))
  tmp4WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,4]))
  tmp5WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,5]))
  tmp6WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,6]))
  tmp7WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,7]))
  tmp8WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,8]))
  tmp9WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,9]))
  tmp10WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,10]))
  tmp11WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,11]))
  tmp12WACCMMax(newiz) = REFORM(MAX(vnWACCM[*,69,iz,12]))

ENDFOR

zplevV = REVERSE(zplev)

tmp1WACCMMin = interpol(tmp1WACCMMin, zplevV, zplevWAM)
tmp2WACCMMin = interpol(tmp2WACCMMin, zplevV, zplevWAM)
tmp3WACCMMin = interpol(tmp3WACCMMin, zplevV, zplevWAM)
tmp4WACCMMin = interpol(tmp4WACCMMin, zplevV, zplevWAM)
tmp5WACCMMin = interpol(tmp5WACCMMin, zplevV, zplevWAM)
tmp6WACCMMin = interpol(tmp6WACCMMin, zplevV, zplevWAM)
tmp7WACCMMin = interpol(tmp7WACCMMin, zplevV, zplevWAM)
tmp8WACCMMin = interpol(tmp8WACCMMin, zplevV, zplevWAM)
tmp9WACCMMin = interpol(tmp9WACCMMin, zplevV, zplevWAM)
tmp10WACCMMin = interpol(tmp10WACCMMin, zplevV, zplevWAM)
tmp11WACCMMin = interpol(tmp11WACCMMin, zplevV, zplevWAM)
tmp12WACCMMin = interpol(tmp12WACCMMin, zplevV, zplevWAM)

tmp1WACCMMax = interpol(tmp1WACCMMax, zplevV, zplevWAM)
tmp2WACCMMax = interpol(tmp2WACCMMax, zplevV, zplevWAM)
tmp3WACCMMax = interpol(tmp3WACCMMax, zplevV, zplevWAM)
tmp4WACCMMax = interpol(tmp4WACCMMax, zplevV, zplevWAM)
tmp5WACCMMax = interpol(tmp5WACCMMax, zplevV, zplevWAM)
tmp6WACCMMax = interpol(tmp6WACCMMax, zplevV, zplevWAM)
tmp7WACCMMax = interpol(tmp7WACCMMax, zplevV, zplevWAM)
tmp8WACCMMax = interpol(tmp8WACCMMax, zplevV, zplevWAM)
tmp9WACCMMax = interpol(tmp9WACCMMax, zplevV, zplevWAM)
tmp10WACCMMax = interpol(tmp10WACCMMax, zplevV, zplevWAM)
tmp11WACCMMax = interpol(tmp11WACCMMax, zplevV, zplevWAM)
tmp12WACCMMax = interpol(tmp12WACCMMax, zplevV, zplevWAM)

;tmpa1 = ABS(FFT(tmp1,-1))
;tmpa2 = ABS(FFT(tmp2,-1))
;tmpa3 = ABS(FFT(tmp3,-1))
;tmpa4 = ABS(FFT(tmp4,-1))
;tmpa5 = ABS(FFT(tmp5,-1))
;tmpa6 = ABS(FFT(tmp6,-1))
;tmpa7 = ABS(FFT(tmp7,-1))
;tmpa8 = ABS(FFT(tmp8,-1))
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmp1WAMMin = FLTARR(nzWAM)
tmp2WAMMin = FLTARR(nzWAM)
tmp3WAMMin = FLTARR(nzWAM)
tmp4WAMMin = FLTARR(nzWAM)
tmp5WAMMin = FLTARR(nzWAM)
tmp6WAMMin = FLTARR(nzWAM)
tmp7WAMMin = FLTARR(nzWAM)
tmp8WAMMin = FLTARR(nzWAM)
tmp9WAMMin = FLTARR(nzWAM)
tmp10WAMMin = FLTARR(nzWAM)
tmp11WAMMin = FLTARR(nzWAM)
tmp12WAMMin = FLTARR(nzWAM)

tmp1WAMMax = FLTARR(nzWAM)
tmp2WAMMax = FLTARR(nzWAM)
tmp3WAMMax = FLTARR(nzWAM)
tmp4WAMMax = FLTARR(nzWAM)
tmp5WAMMax = FLTARR(nzWAM)
tmp6WAMMax = FLTARR(nzWAM)
tmp7WAMMax = FLTARR(nzWAM)
tmp8WAMMax = FLTARR(nzWAM)
tmp9WAMMax = FLTARR(nzWAM)
tmp10WAMMax = FLTARR(nzWAM)
tmp11WAMMax = FLTARR(nzWAM)
tmp12WAMMax = FLTARR(nzWAM)

FOR iz = 0, nzWAM-1 DO BEGIN

  tmp1WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,1]))
  tmp2WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,2]))
  tmp3WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,3]))
  tmp4WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,4]))
  tmp5WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,5]))
  tmp6WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,6]))
  tmp7WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,7]))
  tmp8WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,8]))
  tmp9WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,9]))
  tmp10WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,10]))
  tmp11WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,11]))
  tmp12WAMMin(iz) = REFORM(MIN(vnWAM[*,68,iz,12]))

  tmp1WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,1]))
  tmp2WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,2]))
  tmp3WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,3]))
  tmp4WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,4]))
  tmp5WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,5]))
  tmp6WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,6]))
  tmp7WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,7]))
  tmp8WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,8]))
  tmp9WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,9]))
  tmp10WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,10]))
  tmp11WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,11]))
  tmp12WAMMax(iz) = REFORM(MAX(vnWAM[*,68,iz,12]))

ENDFOR

plot, tmp1WACCMMax,zplevWAM,xrange=[0.0,15.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 1 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp1WAMMin,zplevWAM,linestyle=1
;oplot, tmp1WACCMMin,zplevWAM
oplot, tmp1WAMMax,zplevWAM,linestyle=1

plot, tmp2WACCMMax,zplevWAM,xrange=[0.0,25.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 2 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp2WAMMin,zplevWAM,linestyle=1
;oplot, tmp2WACCMMin,zplevWAM
oplot, tmp2WAMMax,zplevWAM,linestyle=1

plot, tmp3WACCMMax,zplevWAM,xrange=[0.0,25.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 3 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp3WAMMin,zplevWAM,linestyle=1
;oplot, tmp3WACCMMin,zplevWAM
oplot, tmp3WAMMax,zplevWAM,linestyle=1

plot, tmp4WACCMMax,zplevWAM,xrange=[0.0,10.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 4 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp4WAMMin,zplevWAM,linestyle=1
;oplot, tmp4WACCMMin,zplevWAM
oplot, tmp4WAMMax,zplevWAM,linestyle=1

plot, tmp5WACCMMax,zplevWAM,xrange=[0.0,10.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 5 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp5WAMMin,zplevWAM,linestyle=1
;oplot, tmp5WACCMMin,zplevWAM
oplot, tmp5WAMMax,zplevWAM,linestyle=1

plot, tmp6WACCMMax,zplevWAM,xrange=[0.0,5.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 6 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp6WAMMin,zplevWAM,linestyle=1
;oplot, tmp6WACCMMin,zplevWAM
oplot, tmp6WAMMax,zplevWAM,linestyle=1

plot, tmp7WACCMMax,zplevWAM,xrange=[0.0,5.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 7 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp7WAMMin,zplevWAM,linestyle=1
;oplot, tmp7WACCMMin,zplevWAM
oplot, tmp7WAMMax,zplevWAM,linestyle=1

plot, tmp8WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 8 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp8WAMMin,zplevWAM,linestyle=1
;oplot, tmp8WACCMMin,zplevWAM
oplot, tmp8WAMMax,zplevWAM,linestyle=1

plot, tmp9WACCMMax,zplevWAM,xrange=[0.0,5.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 9 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp9WAMMin,zplevWAM,linestyle=1
;oplot, tmp9WACCMMin,zplevWAM
oplot, tmp9WAMMax,zplevWAM,linestyle=1

plot, tmp10WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 10 UT00 40N',$
	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
	 charsize = 1.2 

;oplot, tmp10WAMMin,zplevWAM,linestyle=1
;oplot, tmp10WACCMMin,zplevWAM
oplot, tmp10WAMMax,zplevWAM,linestyle=1

;plot, tmp11WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
;         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 11 UT00 40N',$
;	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
;	 charsize = 1.2 
;
;oplot, tmp11WAMMin,zplevWAM,linestyle=1
;oplot, tmp11WACCMMin,zplevWAM
;oplot, tmp11WAMMax,zplevWAM,linestyle=1
;
;plot, tmp12WACCMMax,zplevWAM,xrange=[0.0,3.0],xstyle=1,$
;         yrange = [1000.0,1.E-09],/ylog,ystyle = 1,title = 'WACCM/WAM Maximum Meridional Wind December Mean WN 12 UT00 40N',$
;	 ytitle = 'Pressure', xtitle = 'Meridional Wind ',ycharsize = 1.0, xcharsize = 1.2,$
;	 charsize = 1.2 
;
;oplot, tmp12WAMMin,zplevWAM,linestyle=1
;oplot, tmp12WACCMMin,zplevWAM
;oplot, tmp12WAMMax,zplevWAM,linestyle=1

DEVICE, /CLOSE

SET_PLOT, 'x'

RETURN

END