| 1 | pro polar_topo |
|---|
| 2 | ; |
|---|
| 3 | ; |
|---|
| 4 | ; |
|---|
| 5 | what_I_plot=0. |
|---|
| 6 | overcontour=0. |
|---|
| 7 | @polar_inc_topo.pro |
|---|
| 8 | SPAWN, '\rm param_plot.idl ; cp polar_inc_topo.pro param_plot.idl' |
|---|
| 9 | ; |
|---|
| 10 | ; |
|---|
| 11 | ; |
|---|
| 12 | zefile=save_ps |
|---|
| 13 | PS_Start, filename=zefile+'.ps' |
|---|
| 14 | print, zefile+'.ps' |
|---|
| 15 | ;!p.multi=[0,2,2] |
|---|
| 16 | ;!P.Charsize = 0.6 |
|---|
| 17 | ; |
|---|
| 18 | ; |
|---|
| 19 | ; |
|---|
| 20 | getcdf, file=filename, charvar=cfield1, invar=field1 |
|---|
| 21 | getcdf, file=filename, charvar=cfield2, invar=field2 |
|---|
| 22 | getcdf, file=filename, charvar=cfield3, invar=field3 |
|---|
| 23 | getcdf, file=filename, charvar='XLONG_M', invar=longi |
|---|
| 24 | getcdf, file=filename, charvar='XLAT_M', invar=lati |
|---|
| 25 | ; |
|---|
| 26 | ; |
|---|
| 27 | ; |
|---|
| 28 | |
|---|
| 29 | ;zefile=save_ps+string(100+ntime,'(I0)') |
|---|
| 30 | ;PS_Start, filename=zefile+'.ps' |
|---|
| 31 | ;!P.Charsize = 1.2 |
|---|
| 32 | ;!p.charthick = 2.0 |
|---|
| 33 | ;!p.thick = 2.0 |
|---|
| 34 | ;!x.thick = 2.0 |
|---|
| 35 | ;!y.thick = 2.0 |
|---|
| 36 | ; |
|---|
| 37 | ; |
|---|
| 38 | ; |
|---|
| 39 | ;latmin = -90. & latmax = 10.0 & lonmin = -180. & lonmax = 180. |
|---|
| 40 | 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] |
|---|
| 41 | ; |
|---|
| 42 | ; |
|---|
| 43 | ; |
|---|
| 44 | what_I_plot = field1 |
|---|
| 45 | overcontour = field2 |
|---|
| 46 | lon = longi |
|---|
| 47 | lat = lati |
|---|
| 48 | overvector_x=0. |
|---|
| 49 | overvector_y=0. |
|---|
| 50 | ; |
|---|
| 51 | ; |
|---|
| 52 | ; |
|---|
| 53 | map_latlon, $ |
|---|
| 54 | what_I_plot, $ ; 2D field |
|---|
| 55 | lon, $ ; 1D latitude |
|---|
| 56 | lat, $ ; 1D longitude |
|---|
| 57 | minfield=minfield_init, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 58 | maxfield=maxfield_init, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 59 | overcontour=overcontour, $ ; another 2D field to overplot with contour lines (=0: no) |
|---|
| 60 | overvector_x=overvector_x, $ ; wind vector - x component (=0: no) |
|---|
| 61 | overvector_y=overvector_y, $ ; wind vector - y component (=0: no) |
|---|
| 62 | ct=pal, $ ; color table (33-rainbow is default) |
|---|
| 63 | colors=colors, $ ; number of colors/levels (32 is default) |
|---|
| 64 | title=title_user, $ ; title of the plot ('' is default) |
|---|
| 65 | format=format ; format of colorbar annotations ('(F6.2)' is default) |
|---|
| 66 | ;; |
|---|
| 67 | ;; |
|---|
| 68 | ;; |
|---|
| 69 | ;loadct, 0 |
|---|
| 70 | ;contour, reform(field3(*,*,ntime)), $ |
|---|
| 71 | ; lon,lat, $ |
|---|
| 72 | ; /overplot, $ |
|---|
| 73 | ; levels=[1,100,200], $ |
|---|
| 74 | ; c_labels=[0,0,0], $ |
|---|
| 75 | ; c_thick=1.5, $ |
|---|
| 76 | ; color=255 |
|---|
| 77 | |
|---|
| 78 | ; |
|---|
| 79 | ; |
|---|
| 80 | ; |
|---|
| 81 | loadct, 0 |
|---|
| 82 | MAP_GRID, $ |
|---|
| 83 | CHARSIZE = 1., $ |
|---|
| 84 | COLOR = 0, $ |
|---|
| 85 | ; LABEL = 2, $ ;; /LABEL or LABEL=2 (one label any 2 grid lines) |
|---|
| 86 | ; LATDEL = 10., $ ;;5 |
|---|
| 87 | ; LONDEL = 30., $ ;;15 |
|---|
| 88 | ;lats=-60, $ |
|---|
| 89 | ; LONLAB = -70., $ ;(latmin+latmax)/2., $ |
|---|
| 90 | ; LATLAB = (lonmin+lonmax)/2., $ |
|---|
| 91 | LABEL = 1, $ ;; /LABEL or LABEL=2 (one label any 2 grid lines) |
|---|
| 92 | LATDEL = 10., $ ;;5 |
|---|
| 93 | LONDEL = 15., $ ;;15 |
|---|
| 94 | LONLAB = -0., $ ;(latmin+latmax)/2., $ |
|---|
| 95 | LATLAB = -0.001, $ |
|---|
| 96 | GLINESTYLE = 2, $ |
|---|
| 97 | GLINETHICK = 0.3 |
|---|
| 98 | ;LONALIGN = 0., $ |
|---|
| 99 | ;LATALIGN = 1. |
|---|
| 100 | |
|---|
| 101 | ;; |
|---|
| 102 | ;; |
|---|
| 103 | PS_End, /PNG |
|---|
| 104 | end |
|---|
