Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

graphwrf.def @ 937

Last change on this file since 937 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
23	extract='yes' ;; EXECUTE ARWPOST.EXE
24	;; WITH THE FOLLOWING PARAMETERS :
25
26	;-------------------------------
27	;-------------------------------
28
29	nam1=" fields = 'tk,HGT,PSFC,U,V,W' " ;; CHANGE HERE FIELDS TO BE EXTRACTED
30	nam1=" fields = 'tk,HGT' "
31	;nam1=" fields = 'HGT,U,V' "
32	nam1=" 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' "
41	nam1=" fields = 'U,V,HGT,SWDOWN' "
42	nam1=" fields = 'U,V,HGT' "
43	nam1=" fields = 'U,V,W' "
44	nam1=" fields = 'TSURF,HGT,U,V,T' "
45	nam1=" fields = 'PSFC,tk,HGT' "
46	nam1=" 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
53	nam2=" start_date = '2024-07-17_00:00:00' " ;; CHANGE HERE START DATE
54	nam3=" 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 "
61	nam4=" interval_seconds = 22200 "
62	;nam4=" interval_seconds = 44400 "
63	nam5=" interp_method = 0 " ;; SET VERT. COORD. 0: eta levels, 1: height (inc) or pressure (dec)
64	;nam5=" interp_method = 1 "
65	tabnam=[0.,6.,0.1] ;; SET HEIGHT/PRESSURE LEVELS WITH [first alt, last alt, step]
66	;tabnam=[-3.,5.,0.1]
67	tabnam=[-5.,10.,0.2]
68
69	nest='d03' ;; SET NEST
70	nest='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
81	num=1 ;; how many plots ?
82	num=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'
126	field1='tk'
127
128	;----------------------------------
129	;----------------------------------
130	; what kind of plot do you want ?
131	;----------------------------------
132	;----------------------------------
133
134	plot='map'
135	level=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
152	level=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
163	nlon=60
164	nlon=150
165	nlon=250
166	nlon=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
188	minfield_init=200. ;; grepALL
189	maxfield_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=lon59.-lon(0)59. ;; grepMAP ;;attention plots successifs
235	;lat=lat59.-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
254	subtitle_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
282	backup_data='no' ;; save ARWpost-generated text files ? default is no
283	already_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/'
289	datafolder='/d5/aslmd/WRF_Mars/DATA/VM/'
290	plotfolder=datafolder
291
292	textitle='generic' ;; one word. default is 'plot'
293	texcomments='sol' ;; one word. default is 'plot'
294
295
296

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format