source: trunk/MESOSCALE/PLOT/SPEC/POLAR/polar_topo.pro.intercomp @ 134

pro polar_topo
;
;
;
what_I_plot=0.
overcontour=0.
@polar_inc_topo.pro
SPAWN, '\rm param_plot.idl ; cp polar_inc_topo.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=lonfield, invar=longi
getcdf, file=filename, charvar=latfield, 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, -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(*,*,0))
overcontour = reform(field2(*,*,0))
lon = reform(longi(*,*,0))
lat = reform(lati(*,*,0))
overvector_x=0.
overvector_y=0.

;contour, what_I_plot, lon, lat
;PS_End, /PNG
;stop

;
;
;
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
                                        ;PS_End, /PNG
                                        ;end
;
;
;
loadct, 0
MAP_GRID, $
          CHARSIZE = 1., $
          COLOR    = 0,   $
        ;  LABEL    = 2,   $   ;; /LABEL or LABEL=2 (one label any 2 grid lines)
        ;          LATDEL   = 10., $   ;;5
        ;          LONDEL   = 30., $   ;;15
        ;lats=-60, $
        ;          LONLAB   = -70., $ ;(latmin+latmax)/2., $
        ;          LATLAB   = (lonmin+lonmax)/2., $
;LATLAB   = lonmax, $
;LABEL    = 1,   $   
        LABEL    = 2,   $   
        LATDEL   = 10., $   
        LONDEL   = 15., $   
        LONLAB   = latmax, $ 
        LATLAB   = -0.001, $
          GLINESTYLE = 2, $
          GLINETHICK = 0.3
          ;LONALIGN = 0., $
          ;LATALIGN = 1.

;;
;;
PS_End, /PNG
end