PRO plotting_waccm_u, un, lat, zplev, uspc,ucomp

; This procedure was created on Tuesday, June 9th, 2009.
; This procedure was created by SOARS Protege McArthur (Mack) Jones Jr.
; This procedure will produce a contour w/ lat. on the x-axis and altitude on the y-axis

;newcolor

SET_PLOT, 'PS'

;DEVICE, SET_FONT = 'Times Bold', /TT_FONT, /COLOR, bits_per_pixel=8
;psfile = '/home/joemci/MackStuff/MyVersion/tcomp_March_21.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_21_MSDTT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MSTDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MSTDT1.ps'

psfile = '/home/joemci/MackStuff/MyVersion/WAM/uamp_March_Mean_ZM_WACCM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_Mean_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MDT1.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile

;tmpa = REFORM(ucomp[72,*,*,1])
tmpa = REFORM(uspc[71,0,*,0:80])

zplevtemp = zplev[0:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;levels = findgen(26)

;levels = findgen(13)

;levels1 = REVERSE(findgen(26))
;
;levels1 = levels1 * (-1)
;
;levels1 = levels1[0:24]
;
;levels = [levels1, levels]

;help, levels
;print, levels
;print, levels1

;theMin = levels[0]

;theMax = levels[50]
;theMax = levels[12]

;nlevels = n_elements(levels)

;ncolors = nlevels + 1

;c_colors = findgen(ncolors)

;bottom = 1 

LOADCT, 39
;LOADCT, 39, ncolors = ncolors, bottom = bottom
;LOADCT, 40, ncolors = ncolors, bottom = bottom

 ;tmpa1_d11 = tmpa1

;SAVE, tmpa1_d11, lat, zplev1, filename='tamp_diurnal_migrating_MarchMeanGWHR.sav'
;SAVE, tmpa1_d11, lat, zplev1, filename='tamp_diurnal_migrating_March11.sav'
;SAVE, tmpa1_d11, lat, zplev1, filename='tamp_diurnal_migrating_March27.sav'

contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(12)*5, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM-X 3.19 Zonal Wind Zonal Mean March Mean ',$
	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], /follow
;, c_labels=intarr(20)+1
tmpa = REFORM(uspc[71,0,*,16:80])

zplevtemp = zplev[16:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(12)*5, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-5],ystyle = 1,title = 'WACCM-X 3.19 Zonal Wind Zonal Mean March Mean ',$
	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(14)+1, /follow

;CONTOUR, tmpa1,lat, zplev1, nlev=16,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-5],$
;            xstyle=1,ystyle=1,/ylog,c_charsize=1.2

tmpa = REFORM(uspc[71,0,*,0:21])

zplevtemp = zplev[0:21]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(30)*2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM-X Migrating Diurnal Tide Temperature March Mean GWHR',$
;	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(7)+1, /follow


;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_colors = c_colors, xrange=[-90,90], xstyle = 1, yrange = [1e3,1e-5],  ystyle = 1, $
;	 title = 'WACCM-X Diurnal Tide s=-1 Temperature March 21', xtitle = 'Latitude', ytitle = 'Pressure (hPa)', /ylog,  $
;	 xtickinterval = 30, /fill, ymargin = [7,2], xmargin = [10,8], ycharsize = 1.02, xcharsize = 1.02, charsize = 1.27

;CONTOUR, tmpa1,lat, zplev1, nlev=16, /fill,c_colors=BYTSCL(indgen(16)),$
;;            xrange=[-90.0,90.0],yrange=[1e-4,1e-9], xstyle=1,ystyle=1,/ylog,$
;	    ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
;	    title = 'WACCM-X 3.19 Migrating Diurnal Tide Temperature March Mean ',$
;	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;	    Title=Title,$
;	    xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;CONTOUR, tmpa1,lat, zplev1, nlev=16,/follow,xrange=[-90.0,90.0], yrange=[1e-4,1e-9],$
;            xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

CONTOUR, tmpa1,lat, zplev1,levels=indgen(12)*5,/follow,xrange=[-90.0,90.0], yrange=[1e-4,1e-9],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Zonal Mean March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(20)+1
;height = ['0','16','32','48','64','80','96','112']

;AXIS, yaxis = 1, ytitle = 'Approx. Altitude (km)', yticks = 7, ytickname = height, ycharsize = 1.3

levels = findgen(6)

levels1 = REVERSE(findgen(6))

levels1 = levels1 * (-1)

levels1 = levels1[0:4]

levels = [levels1, levels]

levels = levels * 5

c_labels = [replicate(1,n_elements(levels))]

;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_labels = c_labels, /overplot

;colorbar, ncolors, theMin, theMax

DEVICE, /CLOSE

;psfile = '/home/joemci/MackStuff/MyVersion/tcomp_March_21.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_21_MSDTT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MSTDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MSTDT1.ps'

psfile = '/home/joemci/MackStuff/MyVersion/uamp_March_Mean_MDT1_WACCM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_Mean_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MDT1.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile

;tmpa = REFORM(ucomp[72,*,*,1])
tmpa = REFORM(uspc[72,1,*,0:80])

zplevtemp = zplev[0:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;levels = findgen(26)

;levels = findgen(13)

;levels1 = REVERSE(findgen(26))
;
;levels1 = levels1 * (-1)
;
;levels1 = levels1[0:24]
;
;levels = [levels1, levels]

;help, levels
;print, levels
;print, levels1

;theMin = levels[0]

;theMax = levels[50]
;theMax = levels[12]

;nlevels = n_elements(levels)

;ncolors = nlevels + 1

;c_colors = findgen(ncolors)

;bottom = 1 

LOADCT, 39
;LOADCT, 39, ncolors = ncolors, bottom = bottom
;LOADCT, 40, ncolors = ncolors, bottom = bottom

 ;tmpa1_d11 = tmpa1

;SAVE, tmpa1_d11, lat, zplev1, filename='tamp_diurnal_migrating_MarchMeanGWHR.sav'
;SAVE, tmpa1_d11, lat, zplev1, filename='tamp_diurnal_migrating_March11.sav'
;SAVE, tmpa1_d11, lat, zplev1, filename='tamp_diurnal_migrating_March27.sav'

contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(21)*8, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM-X 3.19 Zonal Wind Migrating Diurnal Tide March Mean ',$
	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], /follow
;, c_labels=intarr(20)+1
tmpa = REFORM(uspc[72,1,*,16:80])

zplevtemp = zplev[16:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(19)*3, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-5],ystyle = 1,title = 'WACCM-X 3.19 Zonal Wind Migrating Diurnal Tide March Mean ',$
	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(14)+1, /follow

;CONTOUR, tmpa1,lat, zplev1, nlev=16,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-5],$
;            xstyle=1,ystyle=1,/ylog,c_charsize=1.2

tmpa = REFORM(uspc[72,1,*,0:21])

zplevtemp = zplev[0:21]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(30)*2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM-X Migrating Diurnal Tide Temperature March Mean GWHR',$
;	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(7)+1, /follow


;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_colors = c_colors, xrange=[-90,90], xstyle = 1, yrange = [1e3,1e-5],  ystyle = 1, $
;	 title = 'WACCM-X Diurnal Tide s=-1 Temperature March 21', xtitle = 'Latitude', ytitle = 'Pressure (hPa)', /ylog,  $
;	 xtickinterval = 30, /fill, ymargin = [7,2], xmargin = [10,8], ycharsize = 1.02, xcharsize = 1.02, charsize = 1.27

;CONTOUR, tmpa1,lat, zplev1, nlev=16, /fill,c_colors=BYTSCL(indgen(16)),$
;;            xrange=[-90.0,90.0],yrange=[1e-4,1e-9], xstyle=1,ystyle=1,/ylog,$
;	    ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
;	    title = 'WACCM-X 3.19 Migrating Diurnal Tide Temperature March Mean ',$
;	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;	    Title=Title,$
;	    xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;CONTOUR, tmpa1,lat, zplev1, nlev=16,/follow,xrange=[-90.0,90.0], yrange=[1e-4,1e-9],$
;            xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

CONTOUR, tmpa1,lat, zplev1,levels=indgen(25)*8,/follow,xrange=[-90.0,90.0], yrange=[1e-4,1e-9],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Diurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(20)+1
;height = ['0','16','32','48','64','80','96','112']

;AXIS, yaxis = 1, ytitle = 'Approx. Altitude (km)', yticks = 7, ytickname = height, ycharsize = 1.3

levels = findgen(6)

levels1 = REVERSE(findgen(6))

levels1 = levels1 * (-1)

levels1 = levels1[0:4]

levels = [levels1, levels]

levels = levels * 5

c_labels = [replicate(1,n_elements(levels))]

;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_labels = c_labels, /overplot

;colorbar, ncolors, theMin, theMax

DEVICE, /CLOSE
;
;  MSDT
;
;DEVICE, SET_FONT = 'Times Bold', /TT_FONT, /COLOR, bits_per_pixel=8

psfile = '/home/joemci/MackStuff/MyVersion/uamp_March_Mean_MSDT1_WACCM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MSDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MSDT1.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile

;tmpa = REFORM(uspc[73,2,*,*])
tmpa = REFORM(uspc[73,2,*,0:80])

zplevtemp = zplev[0:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;levels = findgen(26)

;levels = findgen(21)/5.

;levels1 = REVERSE(findgen(26))
;
;levels1 = levels1 * (-1)
;
;levels1 = levels1[0:24]
;
;levels = [levels1, levels]

;help, levels
;print, levels
;print, levels1

;theMin = levels[0]

;theMax = levels[50]
;theMax = levels[20]

;nlevels = n_elements(levels)

;ncolors = nlevels + 1

;c_colors = findgen(ncolors)

;bottom = 1 

LOADCT, 39, ncolors = ncolors, bottom = bottom
;LOADCT, 40, ncolors = ncolors, bottom = bottom

 tmpa1_sd11 = tmpa1

CONTOUR, tmpa1,lat, zplev1, levels=findgen(10)*4,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-9],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Semidiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(9)+1
tmpa = REFORM(uspc[73,2,*,16:80])

zplevtemp = zplev[16:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;SAVE, tmpa1_sd11, lat, zplev1, filename='tamp_semidiurnal_migrating_MarchMeanGWHR.sav'
;SAVE, tmpa1_sd11, lat, zplev1, filename='tamp_semidiurnal_migrating_March11.sav'
;SAVE, tmpa1_sd11, lat, zplev1, filename='tamp_semidiurnal_migrating_March27.sav'

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(14)/2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-5],ystyle = 1,title = 'WACCM-X Migrating Semidiurnal Tide Temperature March Mean GWHR',$
;	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(14)+1, /follow

CONTOUR, tmpa1,lat, zplev1, levels=findgen(7)*4,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-5],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Semidiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(7)+1

tmpa = REFORM(uspc[73,2,*,1:21])

zplevtemp = zplev[1:21]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(30)/2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM-X Migrating Semidiurnal Tide Temperature March Mean GWHR',$
;	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(5)+1, /follow
;, c_labels=intarr(5)+1

CONTOUR, tmpa1,lat, zplev1, levels=indgen(10)*4,/follow,xrange=[-90.0,90.0], yrange=[1e-4,1e-9],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Semidiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(10)+1
;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_colors = c_colors, xrange=[-90,90], xstyle = 1, yrange = [1e3,1e-4],  ystyle = 1, $
;	 title = 'WACCM-X Semidiurnal Tide s=-1 Temperature March 21', xtitle = 'Latitude', ytitle = 'Pressure (hPa)', /ylog,  $
;	 xtickinterval = 30, /fill, ymargin = [7,2], xmargin = [10,8], ycharsize = 1.02, xcharsize = 1.02, charsize = 1.27

;height = ['0','16','32','48','64','80','96','112']

;AXIS, yaxis = 1, ytitle = 'Approx. Altitude (km)', yticks = 7, ytickname = height, ycharsize = 1.3

;height = ['0','22','44','66','88','110','132']

;AXIS, yaxis = 1, ytitle = 'Approx. Altitude (km)', yticks = 6, ytickname = height, ycharsize = 1.3

levels = findgen(6)

levels1 = REVERSE(findgen(6))

levels1 = levels1 * (-1)

levels1 = levels1[0:4]

levels = [levels1, levels]

levels = levels * 5

c_labels = [replicate(1,n_elements(levels))]

;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_labels = c_labels, /overplot

;colorbar, ncolors, theMin, theMax

DEVICE, /CLOSE
;
; MTDT 
;
;DEVICE, SET_FONT = 'Times Bold', /TT_FONT, /COLOR, bits_per_pixel=8

psfile = '/home/joemci/MackStuff/MyVersion/uamp_March_Mean_MTDT1_WACCM.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_11_MTDT1.ps'
;psfile = '/home/joemci/MackStuff/MyVersion/tamp_March_27_MTDT1.ps'

DEVICE, /COLOR, /LAND, FILENAME = psfile

;tmpa = REFORM(uspc[72,3,*,*])
tmpa = REFORM(uspc[74,3,*,0:80])

zplevtemp = zplev[0:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;levels = findgen(26)

levels = findgen(13)

;levels1 = REVERSE(findgen(26))
;
;levels1 = levels1 * (-1)
;
;levels1 = levels1[0:24]
;
;levels = [levels1, levels]

help, levels
print, levels
;print, levels1

theMin = levels[0]

;theMax = levels[50]
theMax = levels[12]

nlevels = n_elements(levels)

ncolors = nlevels + 1

c_colors = findgen(ncolors)

bottom = 1 

LOADCT, 39, ncolors = ncolors, bottom = bottom
;LOADCT, 40, ncolors = ncolors, bottom = bottom

 tmpa1_td11 = tmpa1

; SAVE, tmpa1_td11, lat, zplev1, filename='tamp_terdiurnal_migrating_MarchMeanGWHR.sav'
; SAVE, tmpa1_td11, lat, zplev1, filename='tamp_terdiurnal_migrating_March11.sav'
; SAVE, tmpa1_td11, lat, zplev1, filename='tamp_terdiurnal_migrating_March27.sav'

CONTOUR, tmpa1,lat, zplev1, levels=indgen(13)*3,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-9],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Terdiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(13)+1

tmpa = REFORM(uspc[74,3,*,16:80])

zplevtemp = zplev[16:80]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(18)/5, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-3],ystyle = 1,title = 'WACCM-X 3.19 Migrating Terdiurnal Tide March Mean ',$
;	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(5)+1, /follow
;
CONTOUR, tmpa1,lat, zplev1, levels=indgen(5)*2,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-5],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Terdiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(5)+1

tmpa = REFORM(uspc[74,3,*,1:21])

zplevtemp = zplev[1:21]
zplev1 = REVERSE(zplevtemp)

tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(30)/5, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-9],ystyle = 1,title = 'WACCM-X Migrating Terdiurnal Tide Temperature March Mean GWHR',$
;	 xtitle = 'Latitude', ytitle = 'Pressure (hPa)',ycharsize = 1.2, xcharsize = 1.2,$
;	 charsize = 1.3, xticks=9, c_charsize=1.0, xmargin = [10,8], c_labels=intarr(5)+1, /follow
;, c_labels=intarr(5)+1

CONTOUR, tmpa1,lat, zplev1, levels=indgen(13)*3,/follow,xrange=[-90.0,90.0], yrange=[1e-4,1e-9],$
            xstyle=1,ystyle=1,/ylog,ycharsize=1.2,xcharsize=1.2,c_charsize=1.2,xmargin=[10,8],$
	    title = 'WACCM-X 3.19 Zonal Wind Migrating Terdiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'

;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_colors = c_colors, xrange=[-90,90], xstyle = 1, yrange = [1e3,1e-4],  ystyle = 1, $
;	 title = 'WACCM-X Terdiurnal Tide s=-1 Temperature March 21', xtitle = 'Latitude', ytitle = 'Pressure (hPa)', /ylog,  $
;	 xtickinterval = 30, /fill, ymargin = [7,2], xmargin = [10,8], ycharsize = 1.02, xcharsize = 1.02, charsize = 1.27

;height = ['0','16','32','48','64','80','96','112']

;AXIS, yaxis = 1, ytitle = 'Approx. Altitude (km)', yticks = 7, ytickname = height, ycharsize = 1.3

levels = findgen(6)

levels1 = REVERSE(findgen(6))

levels1 = levels1 * (-1)

levels1 = levels1[0:4]

levels = [levels1, levels]

levels = levels * 5

c_labels = [replicate(1,n_elements(levels))]

;CONTOUR, tmpa1, lat, zplev1, levels = levels, c_labels = c_labels, /overplot

;colorbar, ncolors, theMin, theMax

DEVICE, /CLOSE

SET_PLOT, 'x'

END