nx=144
ny=72
nz=57

projcen = fltarr(2)
projcen(0) = -90.
projcen(1) = 0.
perimlat = -45.
xll = fltarr(2)
xur = xll
yll = xll
yur = xll

xll(0) = 0.1
xll(1) = 0.6
xur(0) = 0.4
xur(1) = 0.9
yll(0) = 0.55
yll(1) = 0.05
yur(0) = 0.95
yur(1) = 0.45
pos1 = [xll(0),yll(0),xur(0),yur(0)]
pos2 = [xll(1),yll(0),xur(1),yur(0)]
pos3 = [xll(0),yll(1),xur(0),yur(1)]
pos4 = [xll(1),yll(1),xur(1),yur(1)]

poscb=[0.9,0.25,0.95,0.75]
posplt = [0.1,0.1,0.8,0.9]

;cst = clevels
;cst(where(clevels lt 0)) = 2
;cst(where(clevels ge 0)) = 0

if (projcen(0) ne -90 and projcen(0) ne 90) then begin
  map_set,projcen(0),projcen(1),/stereo,position=posplt,/noborder,$
;  /continents,/grid
  /continents,/grid,/noerase
endif else begin
  if (projcen(0) eq -90.) then begin
    map_set,projcen(0),projcen(1),/stereo,$
      position=posplt,/noborder,limit=[-90.,-180.,perimlat,180.],$
;  /continents,/grid
  /continents,/grid,/noerase
  endif else begin
    map_set,projcen(0),projcen(1),/stereo,$
      position=posplt,/noborder,limit=[perimlat,-180.,90.,180.],$
;  /grid
  /grid,/noerase
  endelse
endelse

nlev=26
clev=findgen(nlev)/(nlev-1)*.2

nlats = 24
iskip = ny-nlats
genpol = fltarr(nx,nlats)
glat = fltarr(nlats)
for j=0,nlats-1 do glat(j) = lat(j)
for j=0,nlats-1 do genpol(*,j) = qjk(*,j,52,3)
contour,genpol,lon,glat,/norm,/noerase,position=posplt,/fill,$
xstyle=5,ystyle=5,levels=clev,/overplot
;xstyle=5,ystyle=5,levels=clevels,c_linest=cst,/overplot

map_set,projcen(0),projcen(1),/stereo,$
position=posplt,/noborder,limit=[perimlat,-180.,90.,180.],$
/grid,/noerase

labpol,posplt,'S',0.5

mycolorbar,nlev,clev,1,poscb,'K/s'

;xyouts,pos2(0)+.5*(pos2(2)-pos2(0)),pos2(1)-.05, $
;'!17Day 25, !7h!17=1500K',$
;align=.5,/norm,charsize = 0.5
;xyouts,pos2(0)+.5*(pos2(2)-pos2(0)),pos2(3)+.03, '!17Potential Vorticity',$
;align=.5,/norm,charsize = 0.5

end