source: trunk/mesoscale/PLOT/SPEC/POLAR/polar_w.pro @ 85

pro polar_w
;
;
;
what_I_plot=0.
overcontour=0.
hem = -1.
@polar_inc_w.pro
SPAWN, '\rm param_plot.idl ; cp polar_inc_w.pro param_plot.idl'
;
;
;
  zefile=save_ps
  PS_Start, filename=zefile+'.ps'
  print, zefile+'.ps'
  ;!p.multi=[0,2,2]
  ;!P.Charsize = 0.6
;
;
;
getcdf, file=filename, charvar=cfield1, invar=field1
getcdf, file=filename, charvar=cfield2, invar=field2
;getcdf, file=filename, charvar=cfield3, invar=field3
getcdf, file=filename, charvar='XLONG', invar=longi
getcdf, file=filename, charvar='XLAT', invar=lati
;
;
;

;zefile=save_ps+string(100+ntime,'(I0)')
;PS_Start, filename=zefile+'.ps'
;!P.Charsize = 1.2
;!p.charthick = 2.0
;!p.thick = 2.0
;!x.thick = 2.0
;!y.thick = 2.0
;
;
;
        ;latmin = -90. & latmax = 10.0 & lonmin = -180. & lonmax = 180.
map_set, hem*90., 0., /isotropic, /azimuthal, /noborder, limit=[latmin,lonmin,latmax,lonmax],title=title_user,/advance ;, position=[0.10, 0.12, 0.90, 0.92]
;
;
;
what_I_plot = reform(field1(*,*,nlevel,ntime))
overcontour = reform(field2(*,*,ntime))
;overcontour = 0.
lon = reform(longi(*,*,ntime))
lat = reform(lati(*,*,ntime))
overvector_x=0.
overvector_y=0.
;
;
;
map_latlon, $
        what_I_plot, $                          ; 2D field
        lon, $                                  ; 1D latitude
        lat, $                                  ; 1D longitude
        minfield=minfield_init, $               ; minimum value of plotted field (=0: calculate)
        maxfield=maxfield_init, $               ; maximum value of plotted field (=0: calculate)
        overcontour=overcontour, $              ; another 2D field to overplot with contour lines (=0: no)
        overvector_x=overvector_x, $            ; wind vector - x component (=0: no)
        overvector_y=overvector_y, $            ; wind vector - y component (=0: no)
        ct=pal, $                               ; color table (33-rainbow is default)
        colors=colors, $                        ; number of colors/levels (32 is default)
        title=title_user, $                     ; title of the plot ('' is default)
        format=format                           ; format of colorbar annotations ('(F6.2)' is default)
;;
;;
;;
;loadct, 0
;contour, reform(field3(*,*,ntime)), $
;       lon,lat, $
;        /overplot, $ 
;       levels=[1,100,200], $
;       c_labels=[0,0,0], $
;        c_thick=1.5, $
;       color=255

;
;
;
loadct, 0
MAP_GRID, CHARSIZE = 0.8, $
          COLOR    = 0,   $
          LABEL    = 1,   $   ;(one label any 1 grid lines)
          LATDEL   = 02., $
          LONDEL   = 45., $
          ;lats=-60, $
          GLINESTYLE = 1, $
          GLINETHICK = 0.1, $
          LONLAB   = -77., $
          LATLAB   = 180.
;;
;;
PS_End, /PNG
end