source: trunk/MESOSCALE_DEV/PLOT/RESERVE/graphwrf.def @ 603

Last change on this file since 603 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

  • Property svn:executable set to *
File size: 8.6 KB
Line 
1;;---------------------------------------
2;
3;  fill here plot parameters then type:
4;
5;  @trace.idl
6;  -or- @report.idl
7;  -or- @movie.idl
8;
9;  this file will be saved !
10;
11;;---------------------------------------
12
13
14;-------------------------------
15;-------------------------------
16; namelist.ARWPOST parameters
17;-------------------------------
18;-------------------------------
19
20;; DEFAULT: ARWPOST.EXE ALREADY RUN (i.e. input.ctl and input.dat already there)
21;; UNCOMMENT FOLLOWING LINE TO EXTRACT FILES
22
23extract='yes'   ;;   EXECUTE ARWPOST.EXE
24                ;;   WITH THE FOLLOWING PARAMETERS :
25
26;-------------------------------
27;-------------------------------
28
29nam1=" fields = 'tk,HGT,PSFC,U,V,W' "           ;; CHANGE HERE FIELDS TO BE EXTRACTED
30nam1=" fields = 'tk,HGT' "
31;nam1=" fields = 'HGT,U,V' "
32nam1=" fields = 'V,W' "
33;nam1=" fields = 'V,W,HGT' "
34;nam1=" fields = 'tk,HGT' "
35;nam1=" fields = 'HGT,ICETOT,RAVE,QH2O_ICE,QH2O' "
36;nam1=" fields = 'U,V,W,QH2O' "
37;nam1=" fields = 'HGT,ICETOT,RAVE,tk,QH2O_ICE' "
38;nam1=" fields = 'HGT,ICETOT,RAVE' "
39;nam1=" fields = 'TSK,HGT' "
40;nam1=" fields = 'U,V,PSFC,tk' "
41nam1=" fields = 'U,V,HGT,SWDOWN' "
42nam1=" fields = 'U,V,HGT' "
43nam1=" fields = 'U,V,W' "
44nam1=" fields = 'TSURF,HGT,U,V,T' "
45nam1=" fields = 'PSFC,tk,HGT' "
46nam1=" fields = 'tk,W' "
47
48;nam2=" start_date = '2024-05-03_06:00:00' "    ;; CHANGE HERE START DATE
49;nam3=" end_date   = '2024-05-04_06:00:00' "    ;; CHANGE HERE END DATE
50;nam2=" start_date = '2024-06-53_06:00:00' "    ;; CHANGE HERE START DATE
51;nam3=" end_date = '2024-06-53_12:00:00' "      ;; CHANGE HERE END DATE
52
53nam2=" start_date = '2024-07-17_00:00:00' "     ;; CHANGE HERE START DATE
54nam3=" end_date = '2024-07-18_00:00:00' "       ;; CHANGE HERE END DATE
55
56
57;nam4=" interval_seconds = 14800 "              ;; SET INTERVAL   2h: 7400s, 4h: 14800s, 6h: 22200s     
58;nam4=" interval_seconds = 7400 "               ;; SET INTERVAL   2h: 7400s, 4h: 14800s, 6h: 22200s
59;nam4=" interval_seconds = 3700 "
60;nam4=" interval_seconds = 88800 "
61nam4=" interval_seconds = 22200 "
62;nam4=" interval_seconds = 44400 "
63nam5=" interp_method = 0 "                      ;; SET VERT. COORD. 0: eta levels, 1: height (inc) or pressure (dec)
64;nam5=" interp_method = 1 "
65tabnam=[0.,6.,0.1]      ;; SET HEIGHT/PRESSURE LEVELS WITH [first alt, last alt, step]
66;tabnam=[-3.,5.,0.1]
67tabnam=[-5.,10.,0.2]
68
69nest='d03'  ;; SET NEST
70nest='d01'  ;; SET NEST
71
72;----------------------------------------------------
73;----------------------------------------------------
74; what do you want to plot ?
75; THIS MUST BE SET ACCORDINGLY TO PREVIOUS SECTION
76;   - do not plot a field not extracted
77;   - do not ask for more plot than data available
78;----------------------------------------------------
79;----------------------------------------------------
80
81num=1                   ;; how many plots ?
82num=6                   ;; how many plots ?
83;num=2                   ;; how many plots ?
84;num=3                  ;; how many plots ?
85;num=12
86;num=13
87;num=4
88;num=7
89
90;@map_sea_level.idl
91;@map_ertel_vorticity.idl
92;field1='HGT'
93;field2='U'
94;winds=['U','V']
95;topo=1
96
97;winds=['U','V'] & field1='SWDOWN' & topo=1
98;winds=['V','W'] & field1='V'
99;field1='ICETOT' & topo=1
100;field1='RAVE' & topo=1
101;field1='QH2O_ICE' & topo=1
102;field1='QH2O' & winds=['U','W']
103;field1='QH2O' & field2='W'
104;field1='QH2O' & winds=['U','W']
105;field1='QH2O_ICE' & field2='tk'
106;field1='tk' & topo=1
107;field1='V' & field2='V'
108;field1='W' & winds=['V','W']
109;field1='tk' & winds=['V','W']
110;field1='tk'
111;field1='W' & topo=1
112;field1='TSK' & topo=1
113;field1='U' & winds=['U','V'] & topo=1
114;winds=['V','W'] & field1='V'
115;winds=['V','W'] & field1='W'
116;field1='U' & winds=['U','V'] & topo=1
117;field1='HGT' & field2='HGT'
118;field1='W' ;& winds=['U','V'] ;& topo=1
119;field1='PSFC' & winds=['U','V'] & field2='W'
120;field1='W' & winds=['U','V'] & field2='PSFC'
121;field1='TSURF' & topo=1 ;& winds=['U','V']
122;field1='T' & topo=1 & winds=['U','V']
123;field1='PSFC' & winds=['tk','HGT'] & topo=1
124;field1='W' & field2='U'
125;field2='tk' & field1='W'
126field1='tk'
127
128;----------------------------------
129;----------------------------------
130; what kind of plot do you want ?
131;----------------------------------
132;----------------------------------
133
134plot='map'     
135level=1         ; first level is 1
136;level=20
137;;level=41
138;level=10
139;level=11
140;;level=1
141;level=9
142;level=8
143;level=3
144;level=49
145;level=5
146;level=15
147;level=17
148;level=5
149;level=10
150;level=7  ;;SEA LEVEL
151;level=30
152level=7
153
154;plot='meridional'     
155;nlon=75        ; first is 1
156;;;;nlon=87
157;nlon=72
158;nlon=70
159;nlon=69
160;nlon=60 
161;nlon=162
162;nlon=100
163nlon=60
164nlon=150
165nlon=250
166nlon=30
167
168;plot='zonal'           
169;;;nlat=145     ; first is 1
170;;nlat=24
171;;nlat=25
172;nlat=100.
173;nlat=165
174
175;@map_wind_stress.idl
176
177
178;----------------------------
179;----------------------------
180; change plot settings
181;----------------------------
182;----------------------------
183
184;; uncomment to modify default plot settings
185;; | | | | | | | | | | | | | | | | | | | | |
186;; v v v v v v v v v v v v v v v v v v v v v
187
188minfield_init=200.              ;; grepALL
189maxfield_init=260.              ;; grepALL
190;minfield_init=0.              ;; grepALL
191;maxfield_init=48.               ;; grepALL
192;intervalx=5.                   ;; grepALL
193;intervaly=5.                    ;; grepALL
194;title_user='Vertical Wind (!Nm!N.s!U-1!N), LT 13, 980m ABG'    ;; grepALL
195;title_axis=['x coord (km)','y coord (km)']     ;; grepALL
196;subtitle_user='LT 13, 980m ABG, Contours: surface pressure anomaly each 0.5 Pa'        ;; grepALL
197;colors=64                      ;; grepALL
198;pal=33                         ;; grepALL
199;pal=0                          ;; grepALL
200;pal=22                         ;; grepALL
201;format='(I0)'                  ;; grepALL
202;windex=30.                     ;; grepALL
203;stride=3                       ;; grepALL
204;flag_cb='false'                ;; grepALL
205;what_I_plot=abs(what_I_plot)    ;; grepALL ??! minfield
206;what_I_plot=what_I_plot+220.
207;isotropic='false'              ;; grepALL
208;overcontour=overcontour-mean(overcontour)      ;; grepALL
209;lev=[-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.]           ;; grepALL
210;lev=[-5.,0.,5.,10.,15.,20.,25.] ;; grepALL
211;lev=[-1.,5.,10.,15.,20.,25.] ;; grepALL
212
213
214;minalt=-3.                     ;; grepSEC
215;maxalt=5.                      ;; grepSEC
216;;;altitude=altitude*1000.      ;; grepSEC
217;minspace=-18.0                 ;; grepSEC
218;maxspace=-9.0                  ;; grepSEC
219;minspace=-11.0                 ;; grepSEC
220;maxspace=-7.0                  ;; grepSEC
221;factor=20.  ;10.                          ;; grepSEC
222;space=space*59. & spacekm='true'   ;; grepSEC
223;minspace=minspace*59./factor       ;; grepSEC
224;maxspace=maxspace*59./factor       ;; grepSEC
225;intervalx=round(intervalx*59./factor)       ;; grepSEC
226
227
228;windowy=[10.,13.]        ;; grepMAP
229;windowx=[14.,17.]      ;; grepMAP
230;windowy=[-12.,-2.]        ;; grepMAP
231;windowx=[-78.,-68.]     ;; grepMAP
232;windowy=[45.,55.]        ;; grepMAP
233;windowx=[-45.,-35.]     ;; grepMAP
234;lon=lon*59.-lon(0)*59.         ;; grepMAP ;;attention plots successifs
235;lat=lat*59.-lat(0)*59.         ;; grepMAP
236;windowy=[0.,14.]       ;; grepMAP
237;windowx=[0.,14.]   ;; grepMAP
238
239
240;lim_max = maxfield_init                                ;; 1/3 grepALL
241;w=where(what_I_plot ge lim_max)                ;; 2/3 grepALL
242;if (w[0] ne -1) then what_I_plot[w]=lim_max    ;; 3/3 grepALL
243
244;lim_min = minfield_init                         ;; 1/3 grepALL
245;w=where(what_I_plot le lim_min)                 ;; 2/3 grepALL
246;if (w[0] ne -1) then what_I_plot[w]=lim_min     ;; 3/3 grepALL
247
248
249;lim_blank = 10                                         ;; 1/3 grepALL
250;w=where(abs(what_I_plot) le lim_blank)                 ;; 2/3 grepALL
251;if (w[0] ne -1) then what_I_plot[w]=missing_value      ;; 3/3 grepALL
252
253;title_user='Sea-level surface pressure'+title_user     ;; grepALL
254subtitle_user='LMD mesoscale model'                     ;; grepALL
255
256
257;; pour tracer juste les vecteurs sur un fond uni ou vide
258;pal=0                                          ;; 1/4 grepMAP
259;what_I_plot(*,*)=what_I_plot(*,*)*0.+0.2       ;; 2/4 grepMAP
260;what_I_plot(0,0)=0.                            ;; 3/4 grepMAP
261;flag_cb='false'                                        ;; 4/4 grepMAP
262
263
264;;; truc pour tracer juste les contours et une zone grisee de topo
265;pal=0                                           ;; 1/9 grepALL
266;w=where(abs(what_I_plot) lt missing_value)     ;; 2/9 grepALL
267;what_I_plot[w]=0.                              ;; 3/9 grepALL
268;w=where(abs(what_I_plot) gt missing_value)     ;; 4/9 grepALL
269;what_I_plot[w]=0.75                            ;; 5/9 grepALL
270;w=where(what_I_plot eq 0.)                     ;; 6/9 grepALL
271;what_I_plot[w]=missing_value                           ;; 7/9 grepALL
272;what_I_plot(0,0)=1. & what_I_plot(1,0)=0.      ;; 8/9 grepALL
273;flag_cb='false'                                ;; 9/9 grepALL
274
275
276;------------------------------
277;------------------------------
278; some last possible settings
279;------------------------------
280;------------------------------
281
282backup_data='no'        ;; save ARWpost-generated text files ? default is no
283already_data='no'       ;; load ARWpost-generated text files ? default is no
284;datafolder='./WHERE_RESULTS_ARE/'      ;; default is ./
285;plotfolder='./WHERE_RESULTS_ARE/'       ;; default is ./
286;plotfolder='./WHERE_PLOTS_ARE/'                ;; default is ./
287;datafolder='/d5/aslmd/WRF_Mars/EN_COURS/newTH20WC/month5/'
288;datafolder='/d5/aslmd/WRF_Mars/ARTICLE/VM/'
289datafolder='/d5/aslmd/WRF_Mars/DATA/VM/'
290plotfolder=datafolder
291
292textitle='generic'              ;; one word.  default is 'plot'
293texcomments='sol'       ;; one word.  default is 'plot'
294
295
296
Note: See TracBrowser for help on using the repository browser.