source: trunk/mesoscale/PLOT/SPEC/GCM/psurf_param_plot.idl @ 85

Last change on this file since 85 was 85, checked in by aslmd, 14 years ago

LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT

File size: 2.5 KB
Line 
1;;;;;*************************************** PLOT TITLES
2        subtitle_user   =       ''
3        title_axis      =       ['East longitude','North latitude']
4        title_user      =       'Surface pressure (Pa)'
5
6;;;;;*************************************** COLOR TABLES
7        flag_cb         =       'true' 
8        format          =       '(I0)'
9        colors          =       64
10        pal             =       16              ;; GOOD: 4, 18, 22, 16, 37, 33, 39
11
12;;;;;*************************************** FILL LIMITS
13        minfield_init   =       0.     
14        maxfield_init   =       800.
15
16;;;;;*************************************** LIMIT TRICKS
17                lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max
18                lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
19;               lim_blank = 2. & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
20
21;;;;;*************************************** WINDS
22        windex          =       10.
23        stride          =       2.
24
25;;;;;*************************************** CONTOUR
26        lev             =       130. + 10.*findgen(20)
27
28;;;;;*************************************** AXIS
29        isotropic       =       'true' 
30        intervalx       =       5.
31        intervaly       =       5.
32
33;;;;;*************************************** MAP LIMITS
34        windowx         =       [-145.,-125.]
35        windowy         =       [5.,30.]
36
37
38
39
40;;;;;***************************************
41;;;;;***************************************
42;;;;;***************************************
43;;;;;***************************************
44
45
46
47
48;;;
49;;; VERTICAL SECTION LIMITS
50;;;
51;minalt=-5.             ;; grepSEC
52;maxalt=40.             ;; grepSEC
53;minspace=0.0           ;; grepSEC
54;maxspace=35.0          ;; grepSEC
55
56;;;
57;;; METRIC UNITS FOR VERTICAL SECTION
58;;;
59;factor=10.        ;; grepSEC
60;space=space*60. & spacekm='true'       ;; grepSEC
61;minspace=minspace*60./factor           ;; grepSEC
62;maxspace=maxspace*60./factor           ;; grepSEC
63;intervalx=round(intervalx*60./factor)   ;; grepSEC
64
65
66;;;
67;;; TRICKS
68;;;
69
70;; pour tracer juste les vecteurs sur un fond uni ou vide
71;pal=0                                          ;; 1/4 grepMAP
72;what_I_plot(*,*)=what_I_plot(*,*)*0.+0.2       ;; 2/4 grepMAP
73;what_I_plot(0,0)=0.                            ;; 3/4 grepMAP
74;flag_cb='false'                                        ;; 4/4 grepMAP
75
76;;; truc pour tracer juste les contours et une zone grisee de topo
77;pal=0                                           ;; 1/9 grepALL
78;w=where(abs(what_I_plot) lt missing_value)     ;; 2/9 grepALL
79;what_I_plot[w]=0.                              ;; 3/9 grepALL
80;w=where(abs(what_I_plot) gt missing_value)     ;; 4/9 grepALL
81;what_I_plot[w]=0.75                            ;; 5/9 grepALL
82;w=where(what_I_plot eq 0.)                     ;; 6/9 grepALL
83;what_I_plot[w]=missing_value                           ;; 7/9 grepALL
84;what_I_plot(0,0)=1. & what_I_plot(1,0)=0.      ;; 8/9 grepALL
85;flag_cb='false'                                ;; 9/9 grepALL
86
87
88
Note: See TracBrowser for help on using the repository browser.