PRO plotting_wam_v, 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/vamp_March_Mean_ZM_WAM.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
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
;tmpa = REFORM(ucomp[72,*,*,1])
tmpa1 = REFORM(uspc[95,0,*,0:149])

zplev1 = zplev[0:149]
;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(16)*2, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional 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
tmpa1 = REFORM(uspc[95,0,*,16:149])

zplev1 = zplev[16:149]

contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(16)*2, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-5],ystyle = 1,title = 'WAM Meridional 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(14)+1
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpa1 = REFORM(uspc[95,0,*,107:149])

zplev1 = zplev[107:149]
;zplev1 = REVERSE(zplevtemp)

;tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(14)*2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-10],ystyle = 1,title = 'WAM Temperature 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(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=8, /fill,c_colors=BYTSCL(indgen(8)),$
;            xrange=[-90.0,90.0],yrange=[1e3,1e-5], xstyle=1,ystyle=1,/ylog,Title=Title,$
;	    xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;CONTOUR, tmpa1,lat, zplev1, nlev=8,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-5],$
;            xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

CONTOUR, tmpa1,lat, zplev1, levels=indgen(16)*2,/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 = 'WAM Merdional Wind Zonal Mean March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(23)+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/vamp_March_Mean_MDT1_WAM.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
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
;tmpa = REFORM(ucomp[72,*,*,1])
tmpa1 = REFORM(uspc[96,1,*,0:149])

zplev1 = zplev[0:149]
;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(19)*8, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-9],ystyle = 1,title = 'WAM Meridional 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
tmpa1 = REFORM(uspc[96,1,*,16:149])

zplev1 = zplev[16:149]

contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(12)*4, xrange=[-90,90], xstyle = 1,$
         yrange = [1e3,1e-5],ystyle = 1,title = 'WAM Meridional 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(14)+1
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpa1 = REFORM(uspc[96,1,*,107:149])

zplev1 = zplev[107:149]
;zplev1 = REVERSE(zplevtemp)

;tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=indgen(14)*2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-10],ystyle = 1,title = 'WAM Temperature 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(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=8, /fill,c_colors=BYTSCL(indgen(8)),$
;            xrange=[-90.0,90.0],yrange=[1e3,1e-5], xstyle=1,ystyle=1,/ylog,Title=Title,$
;	    xtitle=xTitle,ytitle=yTitle,ymargin=[7,3]
;CONTOUR, tmpa1,lat, zplev1, nlev=8,/follow,xrange=[-90.0,90.0], yrange=[1e3,1e-5],$
;            xstyle=1,ystyle=1,/ylog,/overplot,c_charsize=1.2

CONTOUR, tmpa1,lat, zplev1, levels=indgen(18)*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 = 'WAM Merdional Wind Migrating Diurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(23)+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/vamp_March_Mean_MSDT1_WAM.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
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
;tmpa = REFORM(uspc[72,2,*,*])
tmpa1 = REFORM(uspc[97,2,*,0:149])

zplev1 = zplev[0:149]
;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

;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, levels=indgen(14)*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 = 'WAM Meridional Wind Migrating Semidiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(12)+1
tmpa1 = REFORM(uspc[97,2,*,16:149])

zplev1 = zplev[16:149]

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(16)/2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-5],ystyle = 1,title = 'WAM Temperature Migrating Semi-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(5)+1, /follow

CONTOUR, tmpa1,lat, zplev1, levels=indgen(14)*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 = 'WAM Meridional Wind Migrating Semidiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(12)+1	   
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpa1 = REFORM(uspc[97,2,*,107:149])

zplev1 = zplev[107:149]
;zplev1 = REVERSE(zplevtemp)

;tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(10)/2, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-10],ystyle = 1,title = 'WAM Temperature Migrating Semi-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(5)+1, /follow
;, c_labels=intarr(5)+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 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, levels=indgen(14)*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 = 'WAM Merdional Wind Migrating SemiDiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)'
;, c_labels=intarr(12)+1
;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/vamp_March_Mean_MTDT1_WAM.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
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
;tmpa = REFORM(uspc[72,3,*,*])
tmpa1 = REFORM(uspc[98,3,*,0:149])

zplev1 = zplev[0:149]
;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'
print, 'min/max tmpa1 ', min(tmpa1), max(tmpa1)

CONTOUR, tmpa1,lat, zplev1, levels=indgen(7)*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 = 'WAM Merdional Wind Migrating Terdiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(7)+1

tmpa1 = REFORM(uspc[98,3,*,20:149])

zplev1 = zplev[20:149]

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(16)/5, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-5],ystyle = 1,title = 'WAM Temperature Migrating Ter-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(5)+1, /follow
print, 'min/max tmpa1 ', min(tmpa1), max(tmpa1)

CONTOUR, tmpa1,lat, zplev1, levels=indgen(10)*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 = 'WAM Merdional Wind Migrating Terdiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(10)+1
;
; WAM is bottom up as opposed to WACCM top down so no reverse needed
;
tmpa1 = REFORM(uspc[98,3,*,107:149])

zplev1 = zplev[107:149]
;zplev1 = REVERSE(zplevtemp)

;tmpa1 = REVERSE(tmpa, 2)

;contour, tmpa1, lat, zplev1,  /ylog, levels=findgen(10)/5, xrange=[-90,90], xstyle = 1,$
;         yrange = [1e3,1e-10],ystyle = 1,title = 'WAM Temperature Migrating Ter-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(5)+1, /follow
;, c_labels=intarr(5)+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 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, levels=indgen(7)*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 = 'WAM Merdional Wind Migrating TerDiurnal Tide March Mean ',$
	    xtitle = 'Latitude', ytitle = 'Pressure (hPa)', c_labels=intarr(7)+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

SET_PLOT, 'x'

END