1 | load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" |
---|
2 | load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" |
---|
3 | load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" |
---|
4 | |
---|
5 | begin |
---|
6 | |
---|
7 | ; *************************************************** |
---|
8 | ; read parameter_list |
---|
9 | ; *************************************************** |
---|
10 | |
---|
11 | if (isfilepresent("~/.ncl_preferences")) then |
---|
12 | parameter = asciiread("~/.ncl_preferences",75,"string") |
---|
13 | delete(parameter@_FillValue) |
---|
14 | else |
---|
15 | print(" ") |
---|
16 | print("Please copy '.ncl_preferences' into your $home dircetory") |
---|
17 | print(" ") |
---|
18 | exit |
---|
19 | end if |
---|
20 | |
---|
21 | ; *************************************************** |
---|
22 | ; set up default parameter values and strings if not assigned in prompt or parameter list |
---|
23 | ; *************************************************** |
---|
24 | |
---|
25 | if ( .not. isvar("file_in") ) then ; path+name of input file |
---|
26 | if (parameter(7) .EQ. "input file") then |
---|
27 | print(" ") |
---|
28 | print("Please provide input file 'file_in = ' either in prompt or parameter_list") |
---|
29 | print(" ") |
---|
30 | exit |
---|
31 | else |
---|
32 | file_in = parameter(7) |
---|
33 | end if |
---|
34 | end if |
---|
35 | if ( .not. isvar("format_out") ) then ; format of output file |
---|
36 | format_out = "x11" |
---|
37 | if (parameter(9) .NE. "x11") then |
---|
38 | format_out = parameter(9) |
---|
39 | end if |
---|
40 | end if |
---|
41 | if ( .not. isvar("file_out") ) then ; path+name of output file |
---|
42 | file_out = "test" |
---|
43 | if (parameter(11) .NE. "test_ts") then |
---|
44 | file_out = parameter(11) |
---|
45 | end if |
---|
46 | end if |
---|
47 | if ( .not. isvar("no_columns") ) then ; number of plots in one row |
---|
48 | no_columns = 1 |
---|
49 | if (parameter(17) .NE. "1") then |
---|
50 | no_columns = stringtointeger(parameter(17)) |
---|
51 | end if |
---|
52 | end if |
---|
53 | if ( .not. isvar("no_lines") ) then ; number of plot-lines on one sheet |
---|
54 | no_lines = 2 |
---|
55 | if (parameter(19) .NE. "2") then |
---|
56 | no_lines = stringtointeger(parameter(19)) |
---|
57 | end if |
---|
58 | end if |
---|
59 | if ( .not. isvar("combine") ) then ; color of lines |
---|
60 | combine = 0 |
---|
61 | if (parameter(23) .NE. "0") then |
---|
62 | combine = stringtointeger(parameter(23)) |
---|
63 | if (stringtointeger(parameter(23)) .NE. 1) then |
---|
64 | print(" ") |
---|
65 | print("Please set 'combine' to 0 or 1") |
---|
66 | print(" ") |
---|
67 | exit |
---|
68 | end if |
---|
69 | end if |
---|
70 | end if |
---|
71 | if (combine .EQ. 1) then |
---|
72 | if( .not. isvar("c_var") ) then |
---|
73 | if (parameter(27) .EQ. "c_variables") then |
---|
74 | print(" ") |
---|
75 | print("Please select variables for overlaying ('c_var')") |
---|
76 | print(" ") |
---|
77 | exit |
---|
78 | else |
---|
79 | c_var=parameter(27) |
---|
80 | end if |
---|
81 | end if |
---|
82 | end if |
---|
83 | if ( .not. isvar("black") ) then ; color of lines |
---|
84 | black = 0 |
---|
85 | if (parameter(31) .NE. "0") then |
---|
86 | black = stringtointeger(parameter(31)) |
---|
87 | if (stringtointeger(parameter(31)) .NE. 1) then |
---|
88 | print(" ") |
---|
89 | print("Please set 'black' to 0 or 1") |
---|
90 | print(" ") |
---|
91 | exit |
---|
92 | end if |
---|
93 | end if |
---|
94 | end if |
---|
95 | if ( .not. isvar("dash") ) then ; pattern of lines |
---|
96 | dash = 0 |
---|
97 | if (parameter(29) .NE. "0") then |
---|
98 | dash = stringtointeger(parameter(29)) |
---|
99 | if (stringtointeger(parameter(29)) .NE. 1) then |
---|
100 | print(" ") |
---|
101 | print("Please set 'dash' to 0 or 1") |
---|
102 | print(" ") |
---|
103 | exit |
---|
104 | end if |
---|
105 | end if |
---|
106 | end if |
---|
107 | if ( .not. isvar("over") ) then ; switches overlaying plots on |
---|
108 | over = 0 |
---|
109 | if (stringtointeger(parameter(73)) .NE. 0) then |
---|
110 | over = stringtointeger(parameter(73)) |
---|
111 | if (stringtointeger(parameter(73)) .NE. 1) then |
---|
112 | print(" ") |
---|
113 | print("Please set 'over' to 0 or 1") |
---|
114 | print(" ") |
---|
115 | exit |
---|
116 | end if |
---|
117 | end if |
---|
118 | end if |
---|
119 | |
---|
120 | ; *************************************************** |
---|
121 | ; open input file |
---|
122 | ; *************************************************** |
---|
123 | |
---|
124 | f=addfile( file_in,"r") |
---|
125 | |
---|
126 | vNam = getfilevarnames(f) |
---|
127 | print(" ") |
---|
128 | print("Variable on netCDF file: " + vNam) |
---|
129 | print(" ") |
---|
130 | dim = dimsizes(vNam) |
---|
131 | if (dim .EQ. 0) then |
---|
132 | print(" ") |
---|
133 | print("There are no data on file") |
---|
134 | print(" ") |
---|
135 | end if |
---|
136 | |
---|
137 | do varn = dim-1,0,1 |
---|
138 | if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then |
---|
139 | varn=varn-1 |
---|
140 | continue |
---|
141 | end if |
---|
142 | if (vNam(varn) .EQ. "u") |
---|
143 | z_u = f->zu |
---|
144 | break |
---|
145 | else |
---|
146 | if (vNam(varn) .EQ. "v")then |
---|
147 | z_u = f->zv |
---|
148 | break |
---|
149 | else |
---|
150 | if(vNam(varn) .EQ. "pt")then |
---|
151 | z_u = f->zpt |
---|
152 | break |
---|
153 | else |
---|
154 | if(vNam(varn) .EQ. "vpt")then |
---|
155 | z_u = f->zvpt |
---|
156 | break |
---|
157 | else |
---|
158 | if(vNam(varn) .EQ. "lpt")then |
---|
159 | z_u = f->zlpt |
---|
160 | break |
---|
161 | else |
---|
162 | if(vNam(varn) .EQ. "q")then |
---|
163 | z_u = f->zq |
---|
164 | break |
---|
165 | else |
---|
166 | if(vNam(varn) .EQ. "qv")then |
---|
167 | z_u = f->zqv |
---|
168 | break |
---|
169 | else |
---|
170 | if(vNam(varn) .EQ. "ql")then |
---|
171 | z_u = f->zql |
---|
172 | break |
---|
173 | else |
---|
174 | if(vNam(varn) .EQ. "rho")then |
---|
175 | z_u = f->zrho |
---|
176 | break |
---|
177 | else |
---|
178 | if(vNam(varn) .EQ. "s")then |
---|
179 | z_u = f->zs |
---|
180 | break |
---|
181 | else |
---|
182 | if(vNam(varn) .EQ. "sa")then |
---|
183 | z_u = f->zsa |
---|
184 | break |
---|
185 | else |
---|
186 | if(vNam(varn) .EQ. "e")then |
---|
187 | z_u = f->ze |
---|
188 | break |
---|
189 | else |
---|
190 | if(vNam(varn) .EQ. "es")then |
---|
191 | z_u = f->zes |
---|
192 | break |
---|
193 | else |
---|
194 | if(vNam(varn) .EQ. "km")then |
---|
195 | z_u = f->zkm |
---|
196 | break |
---|
197 | else |
---|
198 | if(vNam(varn) .EQ. "kh")then |
---|
199 | z_u = f->zkh |
---|
200 | break |
---|
201 | else |
---|
202 | if(vNam(varn) .EQ. "l")then |
---|
203 | z_u = f->zl |
---|
204 | break |
---|
205 | else |
---|
206 | if(vNam(varn) .EQ. "us2")then |
---|
207 | z_u = f->zus2 |
---|
208 | break |
---|
209 | else |
---|
210 | if(vNam(varn) .EQ. "vs2")then |
---|
211 | z_u = f->zvs2 |
---|
212 | break |
---|
213 | else |
---|
214 | if(vNam(varn) .EQ. "pts2")then |
---|
215 | z_u = f->zpts2 |
---|
216 | break |
---|
217 | else |
---|
218 | if(vNam(varn) .EQ. "wsususodz")then |
---|
219 | z_u = f->zwsususodz |
---|
220 | break |
---|
221 | else |
---|
222 | if(vNam(varn) .EQ. "wspsodz")then |
---|
223 | z_u = f->zwspsodz |
---|
224 | break |
---|
225 | else |
---|
226 | if(vNam(varn) .EQ. "wpeodz")then |
---|
227 | z_u = f->zwpeodz |
---|
228 | break |
---|
229 | end if |
---|
230 | end if |
---|
231 | end if |
---|
232 | end if |
---|
233 | end if |
---|
234 | end if |
---|
235 | end if |
---|
236 | end if |
---|
237 | end if |
---|
238 | end if |
---|
239 | end if |
---|
240 | end if |
---|
241 | end if |
---|
242 | end if |
---|
243 | end if |
---|
244 | end if |
---|
245 | end if |
---|
246 | end if |
---|
247 | end if |
---|
248 | end if |
---|
249 | end if |
---|
250 | end if |
---|
251 | varn=varn-1 |
---|
252 | end do |
---|
253 | |
---|
254 | do varn=dim-1,0,1 |
---|
255 | if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then |
---|
256 | varn=varn-1 |
---|
257 | continue |
---|
258 | end if |
---|
259 | if (vNam(varn) .EQ. "w") |
---|
260 | z_w = f->zw |
---|
261 | break |
---|
262 | else |
---|
263 | if (vNam(varn) .EQ. "wpup")then |
---|
264 | z_w = f->zwpup |
---|
265 | break |
---|
266 | else |
---|
267 | if(vNam(varn) .EQ. "wsus")then |
---|
268 | z_w = f->zwsus |
---|
269 | break |
---|
270 | else |
---|
271 | if(vNam(varn) .EQ. "wu")then |
---|
272 | z_w = f->zwu |
---|
273 | break |
---|
274 | else |
---|
275 | if(vNam(varn) .EQ. "wpvp")then |
---|
276 | z_w = f->zwpvp |
---|
277 | break |
---|
278 | else |
---|
279 | if(vNam(varn) .EQ. "wsvs")then |
---|
280 | z_w = f->zwsvs |
---|
281 | break |
---|
282 | else |
---|
283 | if(vNam(varn) .EQ. "wv")then |
---|
284 | z_w = f->zwv |
---|
285 | break |
---|
286 | else |
---|
287 | if(vNam(varn) .EQ. "wptpp")then |
---|
288 | z_w = f->zwptpp |
---|
289 | break |
---|
290 | else |
---|
291 | if(vNam(varn) .EQ. "wspts")then |
---|
292 | z_w = f->zwspts |
---|
293 | break |
---|
294 | else |
---|
295 | if(vNam(varn) .EQ. "wpt")then |
---|
296 | z_w = f->zwpt |
---|
297 | break |
---|
298 | else |
---|
299 | if(vNam(varn) .EQ. "wsptsBC")then |
---|
300 | z_w = f->zwsptsBC |
---|
301 | break |
---|
302 | else |
---|
303 | if(vNam(varn) .EQ. "wptBC")then |
---|
304 | z_w = f->zwptBC |
---|
305 | break |
---|
306 | else |
---|
307 | if(vNam(varn) .EQ. "wpvptp")then |
---|
308 | z_w = f->zwpvptp |
---|
309 | break |
---|
310 | else |
---|
311 | if(vNam(varn) .EQ. "wsvpts")then |
---|
312 | z_w = f->zwsvpts |
---|
313 | break |
---|
314 | else |
---|
315 | if(vNam(varn) .EQ. "wvpt")then |
---|
316 | z_w = f->zwvpt |
---|
317 | break |
---|
318 | else |
---|
319 | if(vNam(varn) .EQ. "wpqp")then |
---|
320 | z_w = f->zwpqp |
---|
321 | break |
---|
322 | else |
---|
323 | if(vNam(varn) .EQ. "wsqs")then |
---|
324 | z_w = f->zwsqs |
---|
325 | break |
---|
326 | else |
---|
327 | if(vNam(varn) .EQ. "wq")then |
---|
328 | z_w = f->zwq |
---|
329 | break |
---|
330 | else |
---|
331 | if(vNam(varn) .EQ. "wpqvp")then |
---|
332 | z_w = f->zwpqvp |
---|
333 | break |
---|
334 | else |
---|
335 | if(vNam(varn) .EQ. "wsqvs")then |
---|
336 | z_w = f->zwsqvs |
---|
337 | break |
---|
338 | else |
---|
339 | if(vNam(varn) .EQ. "wqv")then |
---|
340 | z_w = f->zwqv |
---|
341 | break |
---|
342 | else |
---|
343 | if(vNam(varn) .EQ. "wpsp")then |
---|
344 | z_w = f->zwpsp |
---|
345 | break |
---|
346 | else |
---|
347 | if(vNam(varn) .EQ. "wsss")then |
---|
348 | z_w = f->zwsss |
---|
349 | break |
---|
350 | else |
---|
351 | if(vNam(varn) .EQ. "ws")then |
---|
352 | z_w = f->zws |
---|
353 | break |
---|
354 | else |
---|
355 | if(vNam(varn) .EQ. "wpsap")then |
---|
356 | z_w = f->zwpsap |
---|
357 | break |
---|
358 | else |
---|
359 | if(vNam(varn) .EQ. "wssas")then |
---|
360 | z_w = f->zwssas |
---|
361 | break |
---|
362 | else |
---|
363 | if(vNam(varn) .EQ. "wsa")then |
---|
364 | z_w = f->zwsa |
---|
365 | break |
---|
366 | else |
---|
367 | if(vNam(varn) .EQ. "wses")then |
---|
368 | z_w = f->zwses |
---|
369 | break |
---|
370 | else |
---|
371 | if(vNam(varn) .EQ. "ws2")then |
---|
372 | z_w = f->zws2 |
---|
373 | break |
---|
374 | else |
---|
375 | if(vNam(varn) .EQ. "ws3")then |
---|
376 | z_w = f->zws3 |
---|
377 | break |
---|
378 | else |
---|
379 | if(vNam(varn) .EQ. "Sw")then |
---|
380 | z_w = f->zSw |
---|
381 | break |
---|
382 | else |
---|
383 | if(vNam(varn) .EQ. "ws2pts")then |
---|
384 | z_w = f->zws2pts |
---|
385 | break |
---|
386 | else |
---|
387 | if(vNam(varn) .EQ. "wspts2")then |
---|
388 | z_w = f->zwspts2 |
---|
389 | break |
---|
390 | end if |
---|
391 | end if |
---|
392 | end if |
---|
393 | end if |
---|
394 | end if |
---|
395 | end if |
---|
396 | end if |
---|
397 | end if |
---|
398 | end if |
---|
399 | end if |
---|
400 | end if |
---|
401 | end if |
---|
402 | end if |
---|
403 | end if |
---|
404 | end if |
---|
405 | end if |
---|
406 | end if |
---|
407 | end if |
---|
408 | end if |
---|
409 | end if |
---|
410 | end if |
---|
411 | end if |
---|
412 | end if |
---|
413 | end if |
---|
414 | end if |
---|
415 | end if |
---|
416 | end if |
---|
417 | end if |
---|
418 | end if |
---|
419 | end if |
---|
420 | end if |
---|
421 | end if |
---|
422 | end if |
---|
423 | varn=varn-1 |
---|
424 | end do |
---|
425 | if ( .not. isvar("z_u") ) then |
---|
426 | zu = 0 |
---|
427 | else |
---|
428 | zu = 1 |
---|
429 | end if |
---|
430 | if ( .not. isvar("z_w") ) then |
---|
431 | zw = 0 |
---|
432 | else |
---|
433 | zw = 1 |
---|
434 | end if |
---|
435 | |
---|
436 | if (zu .EQ. 0 .AND. zw .EQ. 0) then |
---|
437 | do varn=0,dim-1 |
---|
438 | if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then |
---|
439 | check = False |
---|
440 | else |
---|
441 | if (.not. isvar("var")) then |
---|
442 | check = True |
---|
443 | if (parameter(21) .NE. "variables") then |
---|
444 | var=parameter(21) |
---|
445 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
446 | end if |
---|
447 | else |
---|
448 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
449 | end if |
---|
450 | end if |
---|
451 | if (check) |
---|
452 | z = f->$vNam(varn+1)$ |
---|
453 | dimz = dimsizes(z) |
---|
454 | break |
---|
455 | end if |
---|
456 | end do |
---|
457 | end if |
---|
458 | |
---|
459 | dimz = dimsizes(z_w) |
---|
460 | dimz = dimsizes(z_u) |
---|
461 | t_all = f->time |
---|
462 | nt = dimsizes(t_all) |
---|
463 | |
---|
464 | ; **************************************************** |
---|
465 | ; start of time step and different types of mistakes that could be done |
---|
466 | ; **************************************************** |
---|
467 | |
---|
468 | if ( .not. isvar("start_time_step") ) then |
---|
469 | start_time_step = 1 |
---|
470 | if (parameter(13) .NE. "1") then |
---|
471 | if (parameter(13) .LE. "1") |
---|
472 | print(" ") |
---|
473 | print("Begin at least with time step 2") |
---|
474 | print(" ") |
---|
475 | exit |
---|
476 | end if |
---|
477 | if (stringtointeger(parameter(13)) .GE. nt) |
---|
478 | print(" ") |
---|
479 | print("'start_time_step' = "+ parameter(13) +" is greater than available time steps = " + (nt-1)) |
---|
480 | print(" ") |
---|
481 | exit |
---|
482 | end if |
---|
483 | start_time_step = stringtointeger(parameter(13))-1 |
---|
484 | end if |
---|
485 | else |
---|
486 | if (start_time_step .LE. 0) |
---|
487 | print(" ") |
---|
488 | print("Begin at least with time step 1") |
---|
489 | print(" ") |
---|
490 | exit |
---|
491 | end if |
---|
492 | if (start_time_step .GE. nt) |
---|
493 | print(" ") |
---|
494 | print("'start_time_step' = "+ start_time_step +" is greater than available time steps = " + (nt-1)) |
---|
495 | print(" ") |
---|
496 | exit |
---|
497 | end if |
---|
498 | start_time_step = start_time_step - 1 |
---|
499 | end if |
---|
500 | |
---|
501 | ; **************************************************** |
---|
502 | ; end of time step and different types of mistakes that could be done |
---|
503 | ; **************************************************** |
---|
504 | |
---|
505 | if ( .not. isvar("end_time_step") ) then |
---|
506 | end_time_step = nt-1 |
---|
507 | if (parameter(15) .NE. "nt") then |
---|
508 | if (parameter(15) .LE. "0") |
---|
509 | print(" ") |
---|
510 | print("'end_time_step' = "+parameter(15)+ " is too small; 'end_time_step' should be at least 1 ") |
---|
511 | print(" ") |
---|
512 | exit |
---|
513 | end if |
---|
514 | if (stringtointeger(parameter(15)) .GE. nt) |
---|
515 | print(" ") |
---|
516 | print("'end_time_step' = "+ parameter(15) +" is greater than available time steps = " + (nt-1)) |
---|
517 | print(" ") |
---|
518 | exit |
---|
519 | end if |
---|
520 | if (stringtointeger(parameter(15)) .LT. stringtointeger(parameter(13)) ) |
---|
521 | print(" ") |
---|
522 | print("'end_time_step' = "+ parameter(15) +" is lower than 'start_time_step' = "+parameter(13)) |
---|
523 | print(" ") |
---|
524 | exit |
---|
525 | end if |
---|
526 | end_time_step = stringtointeger(parameter(15)) |
---|
527 | end if |
---|
528 | else |
---|
529 | if (end_time_step .LE. 0) |
---|
530 | print(" ") |
---|
531 | print("'end_time_step' = "+end_time_step+ " is too small; 'end_time_step' should be at least 1 ") |
---|
532 | print(" ") |
---|
533 | exit |
---|
534 | end if |
---|
535 | if (end_time_step .GE. nt) |
---|
536 | print(" ") |
---|
537 | print("'end_time_step' = "+ end_time_step +" is greater than available time steps = "+(nt-1)) |
---|
538 | print(" ") |
---|
539 | exit |
---|
540 | end if |
---|
541 | if (end_time_step .LT. start_time_step) |
---|
542 | print(" ") |
---|
543 | print("'end_time_step' = "+end_time_step +" is lower than 'start_time_step' = "+start_time_step) |
---|
544 | print(" ") |
---|
545 | exit |
---|
546 | end if |
---|
547 | end_time_step = end_time_step - 1 |
---|
548 | end if |
---|
549 | |
---|
550 | ; **************************************************** |
---|
551 | ; set up legend and colors |
---|
552 | ; **************************************************** |
---|
553 | |
---|
554 | legend_label=new(nt,double) |
---|
555 | do p=start_time_step,end_time_step |
---|
556 | legend_label(p-start_time_step)=t_all(p) |
---|
557 | end do |
---|
558 | |
---|
559 | np = end_time_step-start_time_step |
---|
560 | if ( black .eq. 0 ) then |
---|
561 | res=True |
---|
562 | res@xyLineColors = ispan(2,237,235/np) |
---|
563 | end if |
---|
564 | |
---|
565 | ; *************************************************** |
---|
566 | ; set up recourses |
---|
567 | ; *************************************************** |
---|
568 | |
---|
569 | res@gsnDraw = False |
---|
570 | res@gsnFrame = False |
---|
571 | res@gsnPaperOrientation = "portrait" |
---|
572 | res@gsnPaperWidth = 8.27 |
---|
573 | res@gsnPaperHeight = 11.69 |
---|
574 | res@gsnPaperMargin = 0.79 |
---|
575 | res@txFont = "helvetica" |
---|
576 | res@tiMainFont = "helvetica" |
---|
577 | res@tiXAxisFont = "helvetica" |
---|
578 | res@tiYAxisFont = "helvetica" |
---|
579 | res@tmXBLabelFont = "helvetica" |
---|
580 | res@tmYLLabelFont = "helvetica" |
---|
581 | res@lgLabelFont = "helvetica" |
---|
582 | res@tmLabelAutoStride = True |
---|
583 | res@pmLegendDisplayMode = "Always" |
---|
584 | res@pmLegendSide = "Top" |
---|
585 | res@xyExplicitLegendLabels = legend_label |
---|
586 | res@pmLegendParallelPosF = 1.15 |
---|
587 | res@pmLegendOrthogonalPosF = -1.0 |
---|
588 | res@pmLegendWidthF = 0.12 |
---|
589 | res@pmLegendHeightF = 0.3 |
---|
590 | res@lgLabelFontHeightF = .02 |
---|
591 | res@txFontHeightF = 0.02 |
---|
592 | res@tiXAxisFontHeightF = 0.02 |
---|
593 | res@tiYAxisFontHeightF = 0.02 |
---|
594 | res@tiXAxisString = " " |
---|
595 | res@tiYAxisString = "Height [z]" |
---|
596 | |
---|
597 | if ( dash .eq. 0 ) then |
---|
598 | res@xyMonoDashPattern = True |
---|
599 | end if |
---|
600 | |
---|
601 | resP = True |
---|
602 | resP@txFont = "helvetica" |
---|
603 | resP@txString = f@title |
---|
604 | resP@txFuncCode = "~" |
---|
605 | resP@txFontHeightF = 0.014 |
---|
606 | |
---|
607 | ; *************************************************** |
---|
608 | ; set up graphics for plot |
---|
609 | ; *************************************************** |
---|
610 | |
---|
611 | plot = new(dim,graphic) |
---|
612 | plot_ = new(dim,graphic) |
---|
613 | |
---|
614 | if (combine .EQ. 1) then |
---|
615 | if ( .not. isvar("number_comb") ) then |
---|
616 | if (parameter(25) .EQ. "0") then |
---|
617 | print(" ") |
---|
618 | print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot") |
---|
619 | print(" ") |
---|
620 | exit |
---|
621 | else |
---|
622 | number_comb=stringtointeger(parameter(25)) |
---|
623 | plot_o = new(number_comb,graphic) |
---|
624 | end if |
---|
625 | else |
---|
626 | if(number_comb .EQ. 2 .OR. number_comb .EQ. 3) then |
---|
627 | plot_o = new(number_comb,graphic) |
---|
628 | else |
---|
629 | print(" ") |
---|
630 | print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot") |
---|
631 | print(" ") |
---|
632 | exit |
---|
633 | end if |
---|
634 | end if |
---|
635 | label=new(number_comb,string) |
---|
636 | color_o=new(number_comb,integer) |
---|
637 | mini=new(number_comb,float) |
---|
638 | maxi=new(number_comb,float) |
---|
639 | end if |
---|
640 | |
---|
641 | wks=gsn_open_wks(format_out,file_out) |
---|
642 | gsn_define_colormap(wks,"rainbow+white") |
---|
643 | |
---|
644 | ; *************************************************** |
---|
645 | ; indicate plot number |
---|
646 | ; *************************************************** |
---|
647 | |
---|
648 | if (combine .EQ. 1) then |
---|
649 | n = 1 |
---|
650 | else |
---|
651 | n = 0 |
---|
652 | end if |
---|
653 | |
---|
654 | ; *************************************************** |
---|
655 | ; set up minimum and maximum height |
---|
656 | ; *************************************************** |
---|
657 | |
---|
658 | if (.not. isvar("min_z")) |
---|
659 | min_z=0 |
---|
660 | if (stringtointeger(parameter(33)) .NE. 0) then |
---|
661 | if (stringtointeger(parameter(33)) .GE. max(z_u) ) then |
---|
662 | print(" ") |
---|
663 | print("Minimum of height ('min_z'="+stringtointeger(parameter(33))+") is greater than available heights (="+max(z_u)+")") |
---|
664 | print(" ") |
---|
665 | exit |
---|
666 | end if |
---|
667 | if (stringtointeger(parameter(33)) .LT. 0 ) then |
---|
668 | print(" ") |
---|
669 | print("Begin minimum of height 'min_z' with 0") |
---|
670 | print(" ") |
---|
671 | exit |
---|
672 | end if |
---|
673 | min_z=stringtointeger(parameter(33)) |
---|
674 | end if |
---|
675 | else |
---|
676 | if (min_z .GE. max(z_u) ) then |
---|
677 | print(" ") |
---|
678 | print("Minimum of height ('min_z'="+min_z+") is greater than available heights (="+max(z_u)+")") |
---|
679 | print(" ") |
---|
680 | exit |
---|
681 | end if |
---|
682 | if (min_z .LT. 0 ) then |
---|
683 | print(" ") |
---|
684 | print("Begin minimum of height 'min_z' with 0") |
---|
685 | print(" ") |
---|
686 | exit |
---|
687 | end if |
---|
688 | end if |
---|
689 | |
---|
690 | if (.not. isvar("max_z")) |
---|
691 | max_z=max(z_u) |
---|
692 | if ((parameter(35)) .NE. "max(z_u)") then |
---|
693 | if (stringtofloat(parameter(35)) .GE. max(z_u) ) then |
---|
694 | print(" ") |
---|
695 | print("Maximum of height ('max_z'="+parameter(35)+") is greater than available heights (="+max(z_u)+")") |
---|
696 | print(" ") |
---|
697 | exit |
---|
698 | end if |
---|
699 | if (stringtointeger(parameter(35)) .LE. 0 ) then |
---|
700 | print(" ") |
---|
701 | print("Maximum of height 'max_z' should be at least 1") |
---|
702 | print(" ") |
---|
703 | exit |
---|
704 | end if |
---|
705 | max_z=stringtointeger(parameter(35)) |
---|
706 | end if |
---|
707 | else |
---|
708 | if (max_z .GE. max(z_pr) ) then |
---|
709 | print(" ") |
---|
710 | print("Maximum of height ('max_z'="+max_z+") is greater than available heights (="+max(z_u)+")") |
---|
711 | print(" ") |
---|
712 | exit |
---|
713 | end if |
---|
714 | if (max_z .LE. 0 ) then |
---|
715 | print(" ") |
---|
716 | print("Maximum of height 'max_z' should be at least 1") |
---|
717 | print(" ") |
---|
718 | exit |
---|
719 | end if |
---|
720 | end if |
---|
721 | |
---|
722 | ; *************************************************** |
---|
723 | ; read data and create plots |
---|
724 | ; *************************************************** |
---|
725 | |
---|
726 | do ti = start_time_step, end_time_step |
---|
727 | if( t_all(ti) .lt. 10^36) then |
---|
728 | start_time_step = ti |
---|
729 | break |
---|
730 | end if |
---|
731 | end do |
---|
732 | |
---|
733 | data = new((/dim,(end_time_step-start_time_step)+1,dimz/),float) |
---|
734 | data_0 = new((/(end_time_step-start_time_step)+1,dimz/),float) |
---|
735 | data_0 = 0.0 |
---|
736 | t = new((/(end_time_step-start_time_step)+1,dimz/),float) |
---|
737 | t = 0.0 |
---|
738 | unit = new(dim,string) |
---|
739 | |
---|
740 | if (over .EQ. 1) then |
---|
741 | plot_u = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
742 | miniu = 1.E27 |
---|
743 | maxiu =-1.E27 |
---|
744 | plot_v = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
745 | miniv = 1.E27 |
---|
746 | maxiv =-1.E27 |
---|
747 | plot_w = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
748 | miniw = 1.E27 |
---|
749 | maxiw =-1.E27 |
---|
750 | plot_pt = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
751 | minipt = 1.E27 |
---|
752 | maxipt =-1.E27 |
---|
753 | plot_vpt = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
754 | minivpt = 1.E27 |
---|
755 | maxivpt =-1.E27 |
---|
756 | plot_lpt = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
757 | minilpt = 1.E27 |
---|
758 | maxilpt =-1.E27 |
---|
759 | plot_q = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
760 | miniq = 1.E27 |
---|
761 | maxiq =-1.E27 |
---|
762 | plot_qv = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
763 | miniqv = 1.E27 |
---|
764 | maxiqv =-1.E27 |
---|
765 | plot_ql = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
766 | miniql = 1.E27 |
---|
767 | maxiql =-1.E27 |
---|
768 | plot_rho = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
769 | plot_s = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
770 | plot_sa = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
771 | plot_e = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
772 | minie = 1.E27 |
---|
773 | maxie =-1.E27 |
---|
774 | plot_es = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
775 | minie = 1.E27 |
---|
776 | maxie =-1.E27 |
---|
777 | plot_km = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
778 | minie = 1.E27 |
---|
779 | maxie =-1.E27 |
---|
780 | plot_kh = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
781 | minie = 1.E27 |
---|
782 | maxie =-1.E27 |
---|
783 | plot_l = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
784 | plot_wpup = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
785 | miniwpup = 1.E27 |
---|
786 | maxiwpup =-1.E27 |
---|
787 | plot_wsus = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
788 | miniwsus = 1.E27 |
---|
789 | maxiwsus =-1.E27 |
---|
790 | plot_wu = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
791 | miniwu = 1.E27 |
---|
792 | maxiwu =-1.E27 |
---|
793 | plot_wpvp = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
794 | miniwpvp = 1.E27 |
---|
795 | maxiwpvp =-1.E27 |
---|
796 | plot_wsvs = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
797 | miniwsvs = 1.E27 |
---|
798 | maxiwsvs =-1.E27 |
---|
799 | plot_wv = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
800 | miniwv = 1.E27 |
---|
801 | maxiwv =-1.E27 |
---|
802 | plot_wpptp = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
803 | miniwpptp = 1.E27 |
---|
804 | maxiwpptp =-1.E27 |
---|
805 | plot_wspts = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
806 | miniwspts = 1.E27 |
---|
807 | maxiwspts =-1.E27 |
---|
808 | plot_wpt = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
809 | miniwpt = 1.E27 |
---|
810 | maxiwpt =-1.E27 |
---|
811 | plot_wsptsBC = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
812 | miniwsptsBC = 1.E27 |
---|
813 | maxiwsptsBC =-1.E27 |
---|
814 | plot_wptBC = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
815 | miniwptBC = 1.E27 |
---|
816 | maxiwptBC =-1.E27 |
---|
817 | plot_wpvptp = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
818 | miniwpvptp = 1.E27 |
---|
819 | maxiwpvptp =-1.E27 |
---|
820 | plot_wsvpts = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
821 | miniwsvpts = 1.E27 |
---|
822 | maxiewsvpts=-1.E27 |
---|
823 | plot_wvpt = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
824 | miniwvpt = 1.E27 |
---|
825 | maxiwvpt =-1.E27 |
---|
826 | plot_wpqp = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
827 | miniwpqp = 1.E27 |
---|
828 | maxiwpqp =-1.E27 |
---|
829 | plot_wsqs = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
830 | miniwsqs = 1.E27 |
---|
831 | maxiwsqs =-1.E27 |
---|
832 | plot_wq = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
833 | miniwq = 1.E27 |
---|
834 | maxiwq =-1.E27 |
---|
835 | plot_wpqvp = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
836 | miniwpqvp = 1.E27 |
---|
837 | maxiwpqvp =-1.E27 |
---|
838 | plot_wsqvs = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
839 | miniwsqvs = 1.E27 |
---|
840 | maxiwsqvs =-1.E27 |
---|
841 | plot_wqv = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
842 | miniwqv = 1.E27 |
---|
843 | maxiwqv =-1.E27 |
---|
844 | plot_wpsp = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
845 | miniwpsp = 1.E27 |
---|
846 | maxiwpsp =-1.E27 |
---|
847 | plot_wsss = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
848 | miniwsss = 1.E27 |
---|
849 | maxiwsss =-1.E27 |
---|
850 | plot_ws = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
851 | miniws = 1.E27 |
---|
852 | maxiws =-1.E27 |
---|
853 | plot_wpsap = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
854 | miniwpsap = 1.E27 |
---|
855 | maxiwpsap =-1.E27 |
---|
856 | plot_wssas = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
857 | miniwssas = 1.E27 |
---|
858 | maxiwssas =-1.E27 |
---|
859 | plot_wsa = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
860 | miniwsa = 1.E27 |
---|
861 | maxiwsa =-1.E27 |
---|
862 | plot_wses = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
863 | plot_us2 = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
864 | minius2 = 1.E27 |
---|
865 | maxius2 =-1.E27 |
---|
866 | plot_vs2 = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
867 | minivs2 = 1.E27 |
---|
868 | maxivs2 =-1.E27 |
---|
869 | plot_ws2 = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
870 | miniws2 = 1.E27 |
---|
871 | maxiws2 =-1.E27 |
---|
872 | plot_pts2 = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
873 | plot_ws3 = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
874 | plot_Sw = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
875 | plot_ws2pts = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
876 | plot_wspts2 = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
877 | plot_wsususodz = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
878 | miniwsususodz = 1.E27 |
---|
879 | maxiwsususodz =-1.E27 |
---|
880 | plot_wspsodz = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
881 | miniwspsodz = 1.E27 |
---|
882 | maxiwspsodz =-1.E27 |
---|
883 | plot_wpeodz = gsn_csm_xy(wks,t,data_0(:,:),res) |
---|
884 | miniwpeodz = 1.E27 |
---|
885 | maxiwpeodz =-1.E27 |
---|
886 | end if |
---|
887 | |
---|
888 | n_o=0 |
---|
889 | do varn = 0,dim-1 |
---|
890 | |
---|
891 | temp = f->$vNam(varn)$ |
---|
892 | |
---|
893 | if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then |
---|
894 | check = False |
---|
895 | else |
---|
896 | if (.not. isvar("var")) then |
---|
897 | check = True |
---|
898 | if (parameter(21) .NE. "variables") then |
---|
899 | var=parameter(21) |
---|
900 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
901 | end if |
---|
902 | else |
---|
903 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
904 | end if |
---|
905 | end if |
---|
906 | |
---|
907 | if (combine .EQ. 1) then |
---|
908 | com=isStrSubset(c_var,","+vNam(varn)+"," ) |
---|
909 | if (com) |
---|
910 | check = False |
---|
911 | data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1) |
---|
912 | unit(varn) = temp@units |
---|
913 | if (n_o .GT. number_comb-1) then |
---|
914 | print(" ") |
---|
915 | print("Please set 'number_comb' to the right number of overlaying variables ('c_var')") |
---|
916 | print(" ") |
---|
917 | exit |
---|
918 | end if |
---|
919 | mini(n_o)=min(data(varn,:,:)) |
---|
920 | maxi(n_o)=max(data(varn,:,:)) |
---|
921 | varn=varn+1 |
---|
922 | n_o=n_o+1 |
---|
923 | end if |
---|
924 | end if |
---|
925 | |
---|
926 | if(check) then |
---|
927 | |
---|
928 | z = f->$vNam(varn+1)$ |
---|
929 | |
---|
930 | unit(varn) = temp@units |
---|
931 | data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1) |
---|
932 | |
---|
933 | if (over .EQ. 0) then |
---|
934 | res@gsnLeftString = vNam(varn) |
---|
935 | res@gsnRightString = unit(varn) |
---|
936 | res@trYMinF = min_z |
---|
937 | res@trYMaxF = max_z |
---|
938 | res@trXMinF = min(data(varn,:,:)) |
---|
939 | res@trXMaxF = max(data(varn,:,:)) |
---|
940 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) |
---|
941 | end if |
---|
942 | |
---|
943 | if (vNam(varn) .EQ. "u") then |
---|
944 | miniu=min(data(varn,:,:)) |
---|
945 | maxiu=max(data(varn,:,:)) |
---|
946 | if (over .EQ. 1) then |
---|
947 | res@xyDashPattern = 0 |
---|
948 | plot_u = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
949 | else |
---|
950 | res@gsnLeftString = vNam(varn) |
---|
951 | res@gsnRightString = unit(varn) |
---|
952 | res@trXMinF = miniu |
---|
953 | res@trXMaxF = maxiu |
---|
954 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
955 | end if |
---|
956 | end if |
---|
957 | if (vNam(varn) .EQ. "v") then |
---|
958 | miniv=min(data(varn,:,:)) |
---|
959 | maxiv=max(data(varn,:,:)) |
---|
960 | if (over .EQ. 1) then |
---|
961 | res@xyDashPattern = 1 |
---|
962 | plot_v = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
963 | else |
---|
964 | res@gsnLeftString = vNam(varn) |
---|
965 | res@gsnRightString = unit(varn) |
---|
966 | res@trXMinF = miniv |
---|
967 | res@trXMaxF = maxiv |
---|
968 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
969 | end if |
---|
970 | end if |
---|
971 | if (vNam(varn) .EQ. "w") then |
---|
972 | miniw=min(data(varn,:,:)) |
---|
973 | maxiw=max(data(varn,:,:)) |
---|
974 | if (over .EQ. 1) then |
---|
975 | res@xyDashPattern = 2 |
---|
976 | plot_w = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
977 | else |
---|
978 | res@gsnLeftString = vNam(varn) |
---|
979 | res@gsnRightString = unit(varn) |
---|
980 | res@trXMinF = miniw |
---|
981 | res@trXMaxF = maxiw |
---|
982 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
983 | end if |
---|
984 | end if |
---|
985 | |
---|
986 | if (vNam(varn) .EQ. "pt") then |
---|
987 | minipt=min(data(varn,:,:)) |
---|
988 | maxipt=max(data(varn,:,:)) |
---|
989 | if (over .EQ. 1) then |
---|
990 | res@xyDashPattern = 0 |
---|
991 | plot_pt = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
992 | else |
---|
993 | res@gsnLeftString = vNam(varn) |
---|
994 | res@gsnRightString = unit(varn) |
---|
995 | res@trXMinF = minipt |
---|
996 | res@trXMaxF = maxipt |
---|
997 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
998 | end if |
---|
999 | end if |
---|
1000 | if (vNam(varn) .EQ. "vpt") then |
---|
1001 | minivpt=min(data(varn,:,:)) |
---|
1002 | maxivpt=max(data(varn,:,:)) |
---|
1003 | if (over .EQ. 1) then |
---|
1004 | res@xyDashPattern = 1 |
---|
1005 | plot_vpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1006 | else |
---|
1007 | res@gsnLeftString = vNam(varn) |
---|
1008 | res@gsnRightString = unit(varn) |
---|
1009 | res@trXMinF = minivpt |
---|
1010 | res@trXMaxF = maxivpt |
---|
1011 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1012 | end if |
---|
1013 | end if |
---|
1014 | if (vNam(varn) .EQ. "lpt") then |
---|
1015 | minilpt=min(data(varn,:,:)) |
---|
1016 | maxilpt=max(data(varn,:,:)) |
---|
1017 | if (over .EQ. 1) then |
---|
1018 | res@xyDashPattern = 2 |
---|
1019 | plot_lpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1020 | else |
---|
1021 | res@gsnLeftString = vNam(varn) |
---|
1022 | res@gsnRightString = unit(varn) |
---|
1023 | res@trXMinF = minilpt |
---|
1024 | res@trXMaxF = maxilpt |
---|
1025 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1026 | end if |
---|
1027 | end if |
---|
1028 | |
---|
1029 | if (vNam(varn) .EQ. "q") then |
---|
1030 | miniq=min(data(varn,:,:)) |
---|
1031 | maxiq=max(data(varn,:,:)) |
---|
1032 | if (over .EQ. 1) then |
---|
1033 | res@xyDashPattern = 0 |
---|
1034 | plot_q = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1035 | else |
---|
1036 | res@gsnLeftString = vNam(varn) |
---|
1037 | res@gsnRightString = unit(varn) |
---|
1038 | res@trXMinF = miniq |
---|
1039 | res@trXMaxF = maxiq |
---|
1040 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1041 | end if |
---|
1042 | end if |
---|
1043 | if (vNam(varn) .EQ. "qv") then |
---|
1044 | miniqv=min(data(varn,:,:)) |
---|
1045 | maxiqv=max(data(varn,:,:)) |
---|
1046 | if (over .EQ. 1) then |
---|
1047 | res@xyDashPattern = 1 |
---|
1048 | plot_qv = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1049 | else |
---|
1050 | res@gsnLeftString = vNam(varn) |
---|
1051 | res@gsnRightString = unit(varn) |
---|
1052 | res@trXMinF = miniqv |
---|
1053 | res@trXMaxF = maxiqv |
---|
1054 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1055 | end if |
---|
1056 | end if |
---|
1057 | if (vNam(varn) .EQ. "ql") then |
---|
1058 | miniql=min(data(varn,:,:)) |
---|
1059 | maxiql=max(data(varn,:,:)) |
---|
1060 | if (over .EQ. 1) then |
---|
1061 | res@xyDashPattern = 2 |
---|
1062 | plot_ql = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1063 | else |
---|
1064 | res@gsnLeftString = vNam(varn) |
---|
1065 | res@gsnRightString = unit(varn) |
---|
1066 | res@trXMinF = miniql |
---|
1067 | res@trXMaxF = maxiql |
---|
1068 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1069 | end if |
---|
1070 | end if |
---|
1071 | |
---|
1072 | if (vNam(varn) .EQ. "rho") then |
---|
1073 | if (over .EQ. 1) then |
---|
1074 | res@gsnLeftString = vNam(varn) |
---|
1075 | res@gsnRightString = unit(varn) |
---|
1076 | res@trXMinF = min(data(varn,:,:)) |
---|
1077 | res@trXMaxF = max(data(varn,:,:)) |
---|
1078 | plot_q = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1079 | else |
---|
1080 | res@gsnLeftString = vNam(varn) |
---|
1081 | res@gsnRightString = unit(varn) |
---|
1082 | res@trXMinF = min(data(varn,:,:)) |
---|
1083 | res@trXMaxF = max(data(varn,:,:)) |
---|
1084 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1085 | end if |
---|
1086 | end if |
---|
1087 | |
---|
1088 | if (vNam(varn) .EQ. "s") then |
---|
1089 | if (over .EQ. 1) then |
---|
1090 | res@gsnLeftString = vNam(varn) |
---|
1091 | res@gsnRightString = unit(varn) |
---|
1092 | res@trXMinF = min(data(varn,:,:)) |
---|
1093 | res@trXMaxF = max(data(varn,:,:)) |
---|
1094 | plot_s = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1095 | else |
---|
1096 | res@gsnLeftString = vNam(varn) |
---|
1097 | res@gsnRightString = unit(varn) |
---|
1098 | res@trXMinF = min(data(varn,:,:)) |
---|
1099 | res@trXMaxF = max(data(varn,:,:)) |
---|
1100 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1101 | end if |
---|
1102 | end if |
---|
1103 | |
---|
1104 | if (vNam(varn) .EQ. "sa") then |
---|
1105 | if (over .EQ. 1) then |
---|
1106 | res@gsnLeftString = vNam(varn) |
---|
1107 | res@gsnRightString = unit(varn) |
---|
1108 | res@trXMinF = min(data(varn,:,:)) |
---|
1109 | res@trXMaxF = max(data(varn,:,:)) |
---|
1110 | plot_sa = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1111 | else |
---|
1112 | res@gsnLeftString = vNam(varn) |
---|
1113 | res@gsnRightString = unit(varn) |
---|
1114 | res@trXMinF = min(data(varn,:,:)) |
---|
1115 | res@trXMaxF = max(data(varn,:,:)) |
---|
1116 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1117 | end if |
---|
1118 | end if |
---|
1119 | |
---|
1120 | if (vNam(varn) .EQ. "e") then |
---|
1121 | minie=min(data(varn,:,:)) |
---|
1122 | maxie=max(data(varn,:,:)) |
---|
1123 | if (over .EQ. 1) then |
---|
1124 | res@xyDashPattern = 0 |
---|
1125 | plot_e = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1126 | else |
---|
1127 | res@gsnLeftString = vNam(varn) |
---|
1128 | res@gsnRightString = unit(varn) |
---|
1129 | res@trXMinF = minie |
---|
1130 | res@trXMaxF = maxie |
---|
1131 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1132 | end if |
---|
1133 | end if |
---|
1134 | if (vNam(varn) .EQ. "es") then |
---|
1135 | minies=min(data(varn,:,:)) |
---|
1136 | maxies=max(data(varn,:,:)) |
---|
1137 | if (over .EQ. 1) then |
---|
1138 | res@xyDashPattern = 1 |
---|
1139 | plot_es = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1140 | else |
---|
1141 | res@gsnLeftString = vNam(varn) |
---|
1142 | res@gsnRightString = unit(varn) |
---|
1143 | res@trXMinF = minies |
---|
1144 | res@trXMaxF = maxies |
---|
1145 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1146 | end if |
---|
1147 | end if |
---|
1148 | |
---|
1149 | if (vNam(varn) .EQ. "km") then |
---|
1150 | minikm=min(data(varn,:,:)) |
---|
1151 | maxikm=max(data(varn,:,:)) |
---|
1152 | if (over .EQ. 1) then |
---|
1153 | res@xyDashPattern = 0 |
---|
1154 | plot_km = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1155 | else |
---|
1156 | res@gsnLeftString = vNam(varn) |
---|
1157 | res@gsnRightString = unit(varn) |
---|
1158 | res@trXMinF = minikm |
---|
1159 | res@trXMaxF = maxikm |
---|
1160 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1161 | end if |
---|
1162 | end if |
---|
1163 | if (vNam(varn) .EQ. "kh") then |
---|
1164 | minikh=min(data(varn,:,:)) |
---|
1165 | maxikh=max(data(varn,:,:)) |
---|
1166 | if (over .EQ. 1) then |
---|
1167 | res@xyDashPattern = 1 |
---|
1168 | plot_kh = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1169 | else |
---|
1170 | res@gsnLeftString = vNam(varn) |
---|
1171 | res@gsnRightString = unit(varn) |
---|
1172 | res@trXMinF = minikh |
---|
1173 | res@trXMaxF = maxikh |
---|
1174 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1175 | end if |
---|
1176 | end if |
---|
1177 | |
---|
1178 | if (vNam(varn) .EQ. "l") then |
---|
1179 | if (over .EQ. 1) then |
---|
1180 | res@gsnLeftString = vNam(varn) |
---|
1181 | res@gsnRightString = unit(varn) |
---|
1182 | res@trXMinF = min(data(varn,:,:)) |
---|
1183 | res@trXMaxF = max(data(varn,:,:)) |
---|
1184 | plot_l = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1185 | else |
---|
1186 | res@gsnLeftString = vNam(varn) |
---|
1187 | res@gsnRightString = unit(varn) |
---|
1188 | res@trXMinF = min(data(varn,:,:)) |
---|
1189 | res@trXMaxF = max(data(varn,:,:)) |
---|
1190 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1191 | end if |
---|
1192 | end if |
---|
1193 | |
---|
1194 | if (vNam(varn) .EQ. "wpup") then |
---|
1195 | miniwpup=min(data(varn,:,:)) |
---|
1196 | maxiwpup=max(data(varn,:,:)) |
---|
1197 | if (over .EQ. 1) then |
---|
1198 | res@xyDashPattern = 0 |
---|
1199 | plot_wpup = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1200 | else |
---|
1201 | res@gsnLeftString = vNam(varn) |
---|
1202 | res@gsnRightString = unit(varn) |
---|
1203 | res@trXMinF = miniwpup |
---|
1204 | res@trXMaxF = maxiwpup |
---|
1205 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1206 | end if |
---|
1207 | end if |
---|
1208 | if (vNam(varn) .EQ. "wsus") then |
---|
1209 | miniwsus=min(data(varn,:,:)) |
---|
1210 | maxiwsus=max(data(varn,:,:)) |
---|
1211 | if (over .EQ. 1) then |
---|
1212 | res@xyDashPattern = 1 |
---|
1213 | plot_wsus = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1214 | else |
---|
1215 | res@gsnLeftString = vNam(varn) |
---|
1216 | res@gsnRightString = unit(varn) |
---|
1217 | res@trXMinF = miniwsus |
---|
1218 | res@trXMaxF = maxiwsus |
---|
1219 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1220 | end if |
---|
1221 | end if |
---|
1222 | if (vNam(varn) .EQ. "wu") then |
---|
1223 | miniwu=min(data(varn,:,:)) |
---|
1224 | maxiwu=max(data(varn,:,:)) |
---|
1225 | if (over .EQ. 1) then |
---|
1226 | res@xyDashPattern = 2 |
---|
1227 | plot_wu = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1228 | else |
---|
1229 | res@gsnLeftString = vNam(varn) |
---|
1230 | res@gsnRightString = unit(varn) |
---|
1231 | res@trXMinF = miniwu |
---|
1232 | res@trXMaxF = maxiwu |
---|
1233 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1234 | end if |
---|
1235 | end if |
---|
1236 | |
---|
1237 | if (vNam(varn) .EQ. "wpvp") then |
---|
1238 | miniwpvp=min(data(varn,:,:)) |
---|
1239 | maxiwpvp=max(data(varn,:,:)) |
---|
1240 | if (over .EQ. 1) then |
---|
1241 | res@xyDashPattern = 0 |
---|
1242 | plot_wpvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1243 | else |
---|
1244 | res@gsnLeftString = vNam(varn) |
---|
1245 | res@gsnRightString = unit(varn) |
---|
1246 | res@trXMinF = miniwpvp |
---|
1247 | res@trXMaxF = maxiwpvp |
---|
1248 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1249 | end if |
---|
1250 | end if |
---|
1251 | if (vNam(varn) .EQ. "wsvs") then |
---|
1252 | miniwsvs=min(data(varn,:,:)) |
---|
1253 | maxiwsvs=max(data(varn,:,:)) |
---|
1254 | if (over .EQ. 1) then |
---|
1255 | res@xyDashPattern = 1 |
---|
1256 | plot_wsvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1257 | else |
---|
1258 | res@gsnLeftString = vNam(varn) |
---|
1259 | res@gsnRightString = unit(varn) |
---|
1260 | res@trXMinF = miniwsvs |
---|
1261 | res@trXMaxF = maxiwsvs |
---|
1262 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1263 | end if |
---|
1264 | end if |
---|
1265 | if (vNam(varn) .EQ. "wv") then |
---|
1266 | miniwv=min(data(varn,:,:)) |
---|
1267 | maxiwv=max(data(varn,:,:)) |
---|
1268 | if (over .EQ. 1) then |
---|
1269 | res@xyDashPattern = 2 |
---|
1270 | plot_wv = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1271 | else |
---|
1272 | res@gsnLeftString = vNam(varn) |
---|
1273 | res@gsnRightString = unit(varn) |
---|
1274 | res@trXMinF = miniwv |
---|
1275 | res@trXMaxF = maxiwv |
---|
1276 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1277 | end if |
---|
1278 | end if |
---|
1279 | |
---|
1280 | if (vNam(varn) .EQ. "wpptp") then |
---|
1281 | miniwpptp=min(data(varn,:,:)) |
---|
1282 | maxiwpptp=max(data(varn,:,:)) |
---|
1283 | if (over .EQ. 1) then |
---|
1284 | res@xyDashPattern = 0 |
---|
1285 | plot_wpptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1286 | else |
---|
1287 | res@gsnLeftString = vNam(varn) |
---|
1288 | res@gsnRightString = unit(varn) |
---|
1289 | res@trXMinF = miniwpptp |
---|
1290 | res@trXMaxF = maxiwpptp |
---|
1291 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1292 | end if |
---|
1293 | end if |
---|
1294 | if (vNam(varn) .EQ. "wspts") then |
---|
1295 | miniwspts=min(data(varn,:,:)) |
---|
1296 | maxiwspts=max(data(varn,:,:)) |
---|
1297 | if (over .EQ. 1) then |
---|
1298 | res@xyDashPattern = 1 |
---|
1299 | plot_wspts = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1300 | else |
---|
1301 | res@gsnLeftString = vNam(varn) |
---|
1302 | res@gsnRightString = unit(varn) |
---|
1303 | res@trXMinF = miniwspts |
---|
1304 | res@trXMaxF = maxiwspts |
---|
1305 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1306 | end if |
---|
1307 | end if |
---|
1308 | if (vNam(varn) .EQ. "wpt") then |
---|
1309 | miniwpt=min(data(varn,:,:)) |
---|
1310 | maxiwpt=max(data(varn,:,:)) |
---|
1311 | if (over .EQ. 1) then |
---|
1312 | res@xyDashPattern = 2 |
---|
1313 | plot_wpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1314 | else |
---|
1315 | res@gsnLeftString = vNam(varn) |
---|
1316 | res@gsnRightString = unit(varn) |
---|
1317 | res@trXMinF = miniwpt |
---|
1318 | res@trXMaxF = maxiwpt |
---|
1319 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1320 | end if |
---|
1321 | end if |
---|
1322 | |
---|
1323 | if (vNam(varn) .EQ. "wsptsBC") then |
---|
1324 | miniwsptsBC=min(data(varn,:,:)) |
---|
1325 | maxiwsptsBC=max(data(varn,:,:)) |
---|
1326 | if (over .EQ. 1) then |
---|
1327 | res@xyDashPattern = 0 |
---|
1328 | plot_wsptsBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1329 | else |
---|
1330 | res@gsnLeftString = vNam(varn) |
---|
1331 | res@gsnRightString = unit(varn) |
---|
1332 | res@trXMinF = miniwsptsBC |
---|
1333 | res@trXMaxF = maxiwsptsBC |
---|
1334 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1335 | end if |
---|
1336 | end if |
---|
1337 | if (vNam(varn) .EQ. "wptBC") then |
---|
1338 | miniwptBC=min(data(varn,:,:)) |
---|
1339 | maxiwptBC=max(data(varn,:,:)) |
---|
1340 | if (over .EQ. 1) then |
---|
1341 | res@xyDashPattern = 1 |
---|
1342 | plot_wptBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1343 | else |
---|
1344 | res@gsnLeftString = vNam(varn) |
---|
1345 | res@gsnRightString = unit(varn) |
---|
1346 | res@trXMinF = miniwptBC |
---|
1347 | res@trXMaxF = maxiwptBC |
---|
1348 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1349 | end if |
---|
1350 | end if |
---|
1351 | |
---|
1352 | if (vNam(varn) .EQ. "wpvptp") then |
---|
1353 | miniwpvptp=min(data(varn,:,:)) |
---|
1354 | maxiwpvptp=max(data(varn,:,:)) |
---|
1355 | if (over .EQ. 1) then |
---|
1356 | res@xyDashPattern = 0 |
---|
1357 | plot_wpvptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1358 | else |
---|
1359 | res@gsnLeftString = vNam(varn) |
---|
1360 | res@gsnRightString = unit(varn) |
---|
1361 | res@trXMinF = miniwpvptp |
---|
1362 | res@trXMaxF = maxiwpvptp |
---|
1363 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1364 | end if |
---|
1365 | end if |
---|
1366 | if (vNam(varn) .EQ. "wsvpts") then |
---|
1367 | miniwsvpts=min(data(varn,:,:)) |
---|
1368 | maxiwsvpts=max(data(varn,:,:)) |
---|
1369 | if (over .EQ. 1) then |
---|
1370 | res@xyDashPattern = 1 |
---|
1371 | plot_wsvpts = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1372 | else |
---|
1373 | res@gsnLeftString = vNam(varn) |
---|
1374 | res@gsnRightString = unit(varn) |
---|
1375 | res@trXMinF = miniwsvpts |
---|
1376 | res@trXMaxF = maxiwsvpts |
---|
1377 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1378 | end if |
---|
1379 | end if |
---|
1380 | if (vNam(varn) .EQ. "wvpt") then |
---|
1381 | miniwvpt=min(data(varn,:,:)) |
---|
1382 | maxiwvpt=max(data(varn,:,:)) |
---|
1383 | if (over .EQ. 1) then |
---|
1384 | res@xyDashPattern = 2 |
---|
1385 | plot_wvpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1386 | else |
---|
1387 | res@gsnLeftString = vNam(varn) |
---|
1388 | res@gsnRightString = unit(varn) |
---|
1389 | res@trXMinF = miniwvpt |
---|
1390 | res@trXMaxF = maxiwvpt |
---|
1391 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1392 | end if |
---|
1393 | end if |
---|
1394 | |
---|
1395 | if (vNam(varn) .EQ. "wpqp") then |
---|
1396 | miniwpqp=min(data(varn,:,:)) |
---|
1397 | maxiwpqp=max(data(varn,:,:)) |
---|
1398 | if (over .EQ. 1) then |
---|
1399 | res@xyDashPattern = 0 |
---|
1400 | plot_wpqp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1401 | else |
---|
1402 | res@gsnLeftString = vNam(varn) |
---|
1403 | res@gsnRightString = unit(varn) |
---|
1404 | res@trXMinF = miniwpqp |
---|
1405 | res@trXMaxF = maxiwpqp |
---|
1406 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1407 | end if |
---|
1408 | end if |
---|
1409 | if (vNam(varn) .EQ. "wsqs") then |
---|
1410 | miniwsqs=min(data(varn,:,:)) |
---|
1411 | maxiwsqs=max(data(varn,:,:)) |
---|
1412 | if (over .EQ. 1) then |
---|
1413 | res@xyDashPattern = 1 |
---|
1414 | plot_wsqs = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1415 | else |
---|
1416 | res@gsnLeftString = vNam(varn) |
---|
1417 | res@gsnRightString = unit(varn) |
---|
1418 | res@trXMinF = miniwsqs |
---|
1419 | res@trXMaxF = maxiwsqs |
---|
1420 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1421 | end if |
---|
1422 | end if |
---|
1423 | if (vNam(varn) .EQ. "wq") then |
---|
1424 | miniwq=min(data(varn,:,:)) |
---|
1425 | maxiwq=max(data(varn,:,:)) |
---|
1426 | if (over .EQ. 1) then |
---|
1427 | res@xyDashPattern = 2 |
---|
1428 | plot_wq = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1429 | else |
---|
1430 | res@gsnLeftString = vNam(varn) |
---|
1431 | res@gsnRightString = unit(varn) |
---|
1432 | res@trXMinF = miniwq |
---|
1433 | res@trXMaxF = maxiwq |
---|
1434 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1435 | end if |
---|
1436 | end if |
---|
1437 | |
---|
1438 | if (vNam(varn) .EQ. "wpqvp") then |
---|
1439 | miniwpqvp=min(data(varn,:,:)) |
---|
1440 | maxiwpqvp=max(data(varn,:,:)) |
---|
1441 | if (over .EQ. 1) then |
---|
1442 | res@xyDashPattern = 0 |
---|
1443 | plot_wpqvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1444 | else |
---|
1445 | res@gsnLeftString = vNam(varn) |
---|
1446 | res@gsnRightString = unit(varn) |
---|
1447 | res@trXMinF = miniwpqvp |
---|
1448 | res@trXMaxF = maxiwpqvp |
---|
1449 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1450 | end if |
---|
1451 | end if |
---|
1452 | if (vNam(varn) .EQ. "wsqvs") then |
---|
1453 | miniwsqvs=min(data(varn,:,:)) |
---|
1454 | maxiwsqvs=max(data(varn,:,:)) |
---|
1455 | if (over .EQ. 1) then |
---|
1456 | res@xyDashPattern = 1 |
---|
1457 | plot_wsqvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1458 | else |
---|
1459 | res@gsnLeftString = vNam(varn) |
---|
1460 | res@gsnRightString = unit(varn) |
---|
1461 | res@trXMinF = miniwsqvs |
---|
1462 | res@trXMaxF = maxiwsqvs |
---|
1463 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1464 | end if |
---|
1465 | end if |
---|
1466 | if (vNam(varn) .EQ. "wqv") then |
---|
1467 | miniwqv=min(data(varn,:,:)) |
---|
1468 | maxiwqv=max(data(varn,:,:)) |
---|
1469 | if (over .EQ. 1) then |
---|
1470 | res@xyDashPattern = 2 |
---|
1471 | plot_wqv = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1472 | else |
---|
1473 | res@gsnLeftString = vNam(varn) |
---|
1474 | res@gsnRightString = unit(varn) |
---|
1475 | res@trXMinF = miniwqv |
---|
1476 | res@trXMaxF = maxiwqv |
---|
1477 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1478 | end if |
---|
1479 | end if |
---|
1480 | |
---|
1481 | if (vNam(varn) .EQ. "wpsp") then |
---|
1482 | miniwpsp=min(data(varn,:,:)) |
---|
1483 | maxiwpsp=max(data(varn,:,:)) |
---|
1484 | if (over .EQ. 1) then |
---|
1485 | res@xyDashPattern = 0 |
---|
1486 | plot_wpsp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1487 | else |
---|
1488 | res@gsnLeftString = vNam(varn) |
---|
1489 | res@gsnRightString = unit(varn) |
---|
1490 | res@trXMinF = miniwpsp |
---|
1491 | res@trXMaxF = maxiwpsp |
---|
1492 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1493 | end if |
---|
1494 | end if |
---|
1495 | if (vNam(varn) .EQ. "wsss") then |
---|
1496 | miniwsss=min(data(varn,:,:)) |
---|
1497 | maxiwsss=max(data(varn,:,:)) |
---|
1498 | if (over .EQ. 1) then |
---|
1499 | res@xyDashPattern = 1 |
---|
1500 | plot_wsss = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1501 | else |
---|
1502 | res@gsnLeftString = vNam(varn) |
---|
1503 | res@gsnRightString = unit(varn) |
---|
1504 | res@trXMinF = miniwsss |
---|
1505 | res@trXMaxF = maxiwsss |
---|
1506 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1507 | end if |
---|
1508 | end if |
---|
1509 | if (vNam(varn) .EQ. "ws") then |
---|
1510 | miniws=min(data(varn,:,:)) |
---|
1511 | maxiws=max(data(varn,:,:)) |
---|
1512 | if (over .EQ. 1) then |
---|
1513 | res@xyDashPattern = 2 |
---|
1514 | plot_ws = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1515 | else |
---|
1516 | res@gsnLeftString = vNam(varn) |
---|
1517 | res@gsnRightString = unit(varn) |
---|
1518 | res@trXMinF = miniws |
---|
1519 | res@trXMaxF = maxiws |
---|
1520 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1521 | end if |
---|
1522 | end if |
---|
1523 | |
---|
1524 | if (vNam(varn) .EQ. "wpsap") then |
---|
1525 | miniwpsap=min(data(varn,:,:)) |
---|
1526 | maxiwpsap=max(data(varn,:,:)) |
---|
1527 | if (over .EQ. 1) then |
---|
1528 | res@xyDashPattern = 0 |
---|
1529 | plot_wpsap = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1530 | else |
---|
1531 | res@gsnLeftString = vNam(varn) |
---|
1532 | res@gsnRightString = unit(varn) |
---|
1533 | res@trXMinF = miniwpsap |
---|
1534 | res@trXMaxF = maxiwpsap |
---|
1535 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1536 | end if |
---|
1537 | end if |
---|
1538 | if (vNam(varn) .EQ. "wssas") then |
---|
1539 | miniwssas=min(data(varn,:,:)) |
---|
1540 | maxiwssas=max(data(varn,:,:)) |
---|
1541 | if (over .EQ. 1) then |
---|
1542 | res@xyDashPattern = 1 |
---|
1543 | plot_wssas = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1544 | else |
---|
1545 | res@gsnLeftString = vNam(varn) |
---|
1546 | res@gsnRightString = unit(varn) |
---|
1547 | res@trXMinF = miniwssas |
---|
1548 | res@trXMaxF = maxiwssas |
---|
1549 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1550 | end if |
---|
1551 | end if |
---|
1552 | if (vNam(varn) .EQ. "wsa") then |
---|
1553 | miniwsa=min(data(varn,:,:)) |
---|
1554 | maxiwsa=max(data(varn,:,:)) |
---|
1555 | if (over .EQ. 1) then |
---|
1556 | res@xyDashPattern = 2 |
---|
1557 | plot_wsa = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1558 | else |
---|
1559 | res@gsnLeftString = vNam(varn) |
---|
1560 | res@gsnRightString = unit(varn) |
---|
1561 | res@trXMinF = miniwsa |
---|
1562 | res@trXMaxF = maxiwsa |
---|
1563 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1564 | end if |
---|
1565 | end if |
---|
1566 | |
---|
1567 | if (vNam(varn) .EQ. "wses") then |
---|
1568 | if (over .EQ. 1) then |
---|
1569 | res@gsnLeftString = vNam(varn) |
---|
1570 | res@gsnRightString = unit(varn) |
---|
1571 | res@trXMinF = min(data(varn,:,:)) |
---|
1572 | res@trXMaxF = max(data(varn,:,:)) |
---|
1573 | plot_wses = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1574 | else |
---|
1575 | res@gsnLeftString = vNam(varn) |
---|
1576 | res@gsnRightString = unit(varn) |
---|
1577 | res@trXMinF = min(data(varn,:,:)) |
---|
1578 | res@trXMaxF = max(data(varn,:,:)) |
---|
1579 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1580 | end if |
---|
1581 | end if |
---|
1582 | |
---|
1583 | if (vNam(varn) .EQ. "us2") then |
---|
1584 | minius2=min(data(varn,:,:)) |
---|
1585 | maxius2=max(data(varn,:,:)) |
---|
1586 | if (over .EQ. 1) then |
---|
1587 | res@xyDashPattern = 0 |
---|
1588 | plot_us2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1589 | else |
---|
1590 | res@gsnLeftString = vNam(varn) |
---|
1591 | res@gsnRightString = unit(varn) |
---|
1592 | res@trXMinF = minius2 |
---|
1593 | res@trXMaxF = maxius2 |
---|
1594 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1595 | end if |
---|
1596 | end if |
---|
1597 | if (vNam(varn) .EQ. "vs2") then |
---|
1598 | minivs2=min(data(varn,:,:)) |
---|
1599 | maxivs2=max(data(varn,:,:)) |
---|
1600 | if (over .EQ. 1) then |
---|
1601 | res@xyDashPattern = 1 |
---|
1602 | plot_vs2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1603 | else |
---|
1604 | res@gsnLeftString = vNam(varn) |
---|
1605 | res@gsnRightString = unit(varn) |
---|
1606 | res@trXMinF = minivs2 |
---|
1607 | res@trXMaxF = maxivs2 |
---|
1608 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1609 | end if |
---|
1610 | end if |
---|
1611 | if (vNam(varn) .EQ. "ws2") then |
---|
1612 | miniws2=min(data(varn,:,:)) |
---|
1613 | maxiws2=max(data(varn,:,:)) |
---|
1614 | if (over .EQ. 1) then |
---|
1615 | res@xyDashPattern = 2 |
---|
1616 | plot_ws2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1617 | else |
---|
1618 | res@gsnLeftString = vNam(varn) |
---|
1619 | res@gsnRightString = unit(varn) |
---|
1620 | res@trXMinF = miniws2 |
---|
1621 | res@trXMaxF = maxiws2 |
---|
1622 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1623 | end if |
---|
1624 | end if |
---|
1625 | |
---|
1626 | if (vNam(varn) .EQ. "pts2") then |
---|
1627 | minipts2=min(data(varn,:,:)) |
---|
1628 | maxipts2=max(data(varn,:,:)) |
---|
1629 | if (over .EQ. 1) then |
---|
1630 | res@gsnLeftString = vNam(varn) |
---|
1631 | res@gsnRightString = unit(varn) |
---|
1632 | res@trXMinF = minipts2 |
---|
1633 | res@trXMaxF = maxipts2 |
---|
1634 | plot_pts2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1635 | else |
---|
1636 | res@gsnLeftString = vNam(varn) |
---|
1637 | res@gsnRightString = unit(varn) |
---|
1638 | res@trXMinF = minipts2 |
---|
1639 | res@trXMaxF = maxipts2 |
---|
1640 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1641 | end if |
---|
1642 | end if |
---|
1643 | |
---|
1644 | if (vNam(varn) .EQ. "ws3") then |
---|
1645 | miniws3=min(data(varn,:,:)) |
---|
1646 | maxiws3=max(data(varn,:,:)) |
---|
1647 | if (over .EQ. 1) then |
---|
1648 | res@gsnLeftString = vNam(varn) |
---|
1649 | res@gsnRightString = unit(varn) |
---|
1650 | res@trXMinF = miniws3 |
---|
1651 | res@trXMaxF = maxiws3 |
---|
1652 | plot_ws3 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1653 | else |
---|
1654 | res@gsnLeftString = vNam(varn) |
---|
1655 | res@gsnRightString = unit(varn) |
---|
1656 | res@trXMinF = miniws3 |
---|
1657 | res@trXMaxF = maxiws3 |
---|
1658 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1659 | end if |
---|
1660 | end if |
---|
1661 | |
---|
1662 | if (vNam(varn) .EQ. "Sw") then |
---|
1663 | miniSw=min(data(varn,:,:)) |
---|
1664 | maxiSw=max(data(varn,:,:)) |
---|
1665 | if (over .EQ. 1) then |
---|
1666 | res@gsnLeftString = vNam(varn) |
---|
1667 | res@gsnRightString = unit(varn) |
---|
1668 | res@trXMinF = miniSw |
---|
1669 | res@trXMaxF = maxiSw |
---|
1670 | plot_Sw = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1671 | else |
---|
1672 | res@gsnLeftString = vNam(varn) |
---|
1673 | res@gsnRightString = unit(varn) |
---|
1674 | res@trXMinF = miniSw |
---|
1675 | res@trXMaxF = maxiSw |
---|
1676 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1677 | end if |
---|
1678 | end if |
---|
1679 | |
---|
1680 | if (vNam(varn) .EQ. "ws2pts") then |
---|
1681 | miniws2pts=min(data(varn,:,:)) |
---|
1682 | maxiws2pts=max(data(varn,:,:)) |
---|
1683 | if (over .EQ. 1) then |
---|
1684 | res@gsnLeftString = vNam(varn) |
---|
1685 | res@gsnRightString = unit(varn) |
---|
1686 | res@trXMinF = miniws2pts |
---|
1687 | res@trXMaxF = maxiws2pts |
---|
1688 | plot_ws2pts = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1689 | else |
---|
1690 | res@gsnLeftString = vNam(varn) |
---|
1691 | res@gsnRightString = unit(varn) |
---|
1692 | res@trXMinF = miniws2pts |
---|
1693 | res@trXMaxF = maxiws2pts |
---|
1694 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1695 | end if |
---|
1696 | end if |
---|
1697 | |
---|
1698 | if (vNam(varn) .EQ. "wspts2") then |
---|
1699 | miniwspts2=min(data(varn,:,:)) |
---|
1700 | maxiwspts2=max(data(varn,:,:)) |
---|
1701 | if (over .EQ. 1) then |
---|
1702 | res@gsnLeftString = vNam(varn) |
---|
1703 | res@gsnRightString = unit(varn) |
---|
1704 | res@trXMinF = miniwspts2 |
---|
1705 | res@trXMaxF = maxiwspts2 |
---|
1706 | plot_wspts2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1707 | else |
---|
1708 | res@gsnLeftString = vNam(varn) |
---|
1709 | res@gsnRightString = unit(varn) |
---|
1710 | res@trXMinF = miniwspts2 |
---|
1711 | res@trXMaxF = maxiwspts2 |
---|
1712 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1713 | end if |
---|
1714 | end if |
---|
1715 | |
---|
1716 | if (vNam(varn) .EQ. "wsususodz") then |
---|
1717 | miniwsususodz=min(data(varn,:,:)) |
---|
1718 | maxiwsususodz=max(data(varn,:,:)) |
---|
1719 | if (over .EQ. 1) then |
---|
1720 | res@xyDashPattern = 0 |
---|
1721 | plot_wsususodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1722 | else |
---|
1723 | res@gsnLeftString = vNam(varn) |
---|
1724 | res@gsnRightString = unit(varn) |
---|
1725 | res@trXMinF = miniwsususodz |
---|
1726 | res@trXMaxF = maxiwsususodz |
---|
1727 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1728 | end if |
---|
1729 | end if |
---|
1730 | if (vNam(varn) .EQ. "wspsodz") then |
---|
1731 | miniwspsodz=min(data(varn,:,:)) |
---|
1732 | maxiwspsodz=max(data(varn,:,:)) |
---|
1733 | if (over .EQ. 1) then |
---|
1734 | res@xyDashPattern = 1 |
---|
1735 | plot_wspsodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1736 | else |
---|
1737 | res@gsnLeftString = vNam(varn) |
---|
1738 | res@gsnRightString = unit(varn) |
---|
1739 | res@trXMinF = miniwspsodz |
---|
1740 | res@trXMaxF = maxiwspsodz |
---|
1741 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1742 | end if |
---|
1743 | end if |
---|
1744 | if (vNam(varn) .EQ. "wpeodz") then |
---|
1745 | miniwpeodz=min(data(varn,:,:)) |
---|
1746 | maxiwpeodz=max(data(varn,:,:)) |
---|
1747 | if (over .EQ. 1) then |
---|
1748 | res@xyDashPattern = 2 |
---|
1749 | plot_wpeodz = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1750 | else |
---|
1751 | res@gsnLeftString = vNam(varn) |
---|
1752 | res@gsnRightString = unit(varn) |
---|
1753 | res@trXMinF = miniwpeodz |
---|
1754 | res@trXMaxF = maxiwpeodz |
---|
1755 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1756 | end if |
---|
1757 | end if |
---|
1758 | if (over .EQ. 0) then |
---|
1759 | n=n+1 |
---|
1760 | end if |
---|
1761 | varn=varn+1 |
---|
1762 | end if |
---|
1763 | delete(temp) |
---|
1764 | end do |
---|
1765 | |
---|
1766 | if (over .EQ. 1 ) then |
---|
1767 | |
---|
1768 | overlay(plot_u,plot_v) |
---|
1769 | overlay(plot_u,plot_w) |
---|
1770 | u=0 |
---|
1771 | overlay(plot_pt,plot_vpt) |
---|
1772 | overlay(plot_pt,plot_lpt) |
---|
1773 | pt=0 |
---|
1774 | overlay(plot_q,plot_qv) |
---|
1775 | overlay(plot_q,plot_ql) |
---|
1776 | q=0 |
---|
1777 | overlay(plot_e,plot_es) |
---|
1778 | e=0 |
---|
1779 | overlay(plot_km,plot_kh) |
---|
1780 | km=0 |
---|
1781 | overlay(plot_wpup,plot_wsus) |
---|
1782 | overlay(plot_wpup,plot_wu) |
---|
1783 | wpup=0 |
---|
1784 | overlay(plot_wpvp,plot_wsvs) |
---|
1785 | overlay(plot_wpvp,plot_wv) |
---|
1786 | wpvp=0 |
---|
1787 | overlay(plot_wpptp,plot_wspts) |
---|
1788 | overlay(plot_wpptp,plot_wpt) |
---|
1789 | wpptp=0 |
---|
1790 | overlay(plot_wsptsBC,plot_wptBC) |
---|
1791 | wsptsBC=0 |
---|
1792 | overlay(plot_wpvptp,plot_wsvpts) |
---|
1793 | overlay(plot_wpvptp,plot_wvpt) |
---|
1794 | wpvptp=0 |
---|
1795 | overlay(plot_wpqp,plot_wsqs) |
---|
1796 | overlay(plot_wpqp,plot_wq) |
---|
1797 | wpqp=0 |
---|
1798 | overlay(plot_wpqvp,plot_wsqvs) |
---|
1799 | overlay(plot_wpqvp,plot_wqv) |
---|
1800 | wpqvp=0 |
---|
1801 | overlay(plot_wpsp,plot_wsss) |
---|
1802 | overlay(plot_wpsp,plot_ws) |
---|
1803 | wpsp=0 |
---|
1804 | overlay(plot_wpsap,plot_wssas) |
---|
1805 | overlay(plot_wpsap,plot_wsa) |
---|
1806 | wpsap=0 |
---|
1807 | overlay(plot_us2,plot_vs2) |
---|
1808 | overlay(plot_us2,plot_ws2) |
---|
1809 | us2=0 |
---|
1810 | overlay(plot_wsususodz,plot_wspsodz) |
---|
1811 | overlay(plot_wsususodz,plot_wpeodz) |
---|
1812 | wsususodz=0 |
---|
1813 | |
---|
1814 | end if |
---|
1815 | |
---|
1816 | if (over .EQ. 1) then |
---|
1817 | |
---|
1818 | do varn = 0,dim-1 |
---|
1819 | |
---|
1820 | if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then |
---|
1821 | check = False |
---|
1822 | else |
---|
1823 | if (.not. isvar("var")) then |
---|
1824 | check = True |
---|
1825 | if (parameter(21) .NE. "variables") then |
---|
1826 | var=parameter(21) |
---|
1827 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
1828 | end if |
---|
1829 | else |
---|
1830 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
1831 | end if |
---|
1832 | end if |
---|
1833 | |
---|
1834 | if (check) |
---|
1835 | |
---|
1836 | z = f->$vNam(varn+1)$ |
---|
1837 | |
---|
1838 | ;res@gsnLeftString = vNam(varn) |
---|
1839 | res@gsnRightString = unit(varn) |
---|
1840 | res@trYMinF = min_z |
---|
1841 | res@trYMaxF = max_z |
---|
1842 | res@trXMinF = min(data(varn,:,:)) |
---|
1843 | res@trXMaxF = max(data(varn,:,:)) |
---|
1844 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) |
---|
1845 | |
---|
1846 | if (vNam(varn) .EQ. "u" .OR. vNam(varn) .EQ. "v" .OR. vNam(varn) .EQ. "w") then |
---|
1847 | if (u .EQ. 0) then |
---|
1848 | res@gsnLeftString = "u, v and w" |
---|
1849 | res@gsnRightString = unit(varn) |
---|
1850 | res@trXMinF = min((/miniu,miniv,miniw/)) |
---|
1851 | res@trXMaxF = max((/maxiu,maxiv,maxiw/)) |
---|
1852 | if (vNam(varn) .EQ. "w") then |
---|
1853 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
1854 | else |
---|
1855 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1856 | end if |
---|
1857 | |
---|
1858 | ; *************************************************** |
---|
1859 | ; legend for combined plot |
---|
1860 | ; *************************************************** |
---|
1861 | |
---|
1862 | lgres = True |
---|
1863 | lgMonoDashIndex = False |
---|
1864 | lgres@lgLabelFont = "helvetica" |
---|
1865 | lgres@lgLabelFontHeightF = .1 |
---|
1866 | lgres@vpWidthF = 0.12 |
---|
1867 | lgres@vpHeightF = 0.1 |
---|
1868 | lgres@lgDashIndexes = (/0,1,2/) |
---|
1869 | lbid = gsn_create_legend(wks,3,(/"u","v","w"/),lgres) |
---|
1870 | |
---|
1871 | amres = True |
---|
1872 | amres@amParallelPosF = 0.65 |
---|
1873 | amres@amOrthogonalPosF = -0.2 |
---|
1874 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
1875 | overlay(plot(n),plot_u) |
---|
1876 | u=1 |
---|
1877 | else |
---|
1878 | varn=varn+1 |
---|
1879 | continue |
---|
1880 | end if |
---|
1881 | end if |
---|
1882 | |
---|
1883 | if (vNam(varn) .EQ. "pt" .OR. vNam(varn) .EQ. "vpt" .OR. vNam(varn) .EQ. "lpt") then |
---|
1884 | if (pt .EQ. 0) then |
---|
1885 | res@gsnLeftString = "pt, vpt and lpt" |
---|
1886 | res@gsnRightString = unit(varn) |
---|
1887 | res@trXMinF = min((/minipt,minivpt,minilpt/)) |
---|
1888 | res@trXMaxF = max((/maxipt,maxivpt,maxilpt/)) |
---|
1889 | |
---|
1890 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1891 | |
---|
1892 | ; *************************************************** |
---|
1893 | ; legend for combined plot |
---|
1894 | ; *************************************************** |
---|
1895 | |
---|
1896 | lgres = True |
---|
1897 | lgMonoDashIndex = False |
---|
1898 | lgres@lgLabelFont = "helvetica" |
---|
1899 | lgres@lgLabelFontHeightF = .1 |
---|
1900 | lgres@vpWidthF = 0.12 |
---|
1901 | lgres@vpHeightF = 0.1 |
---|
1902 | lgres@lgDashIndexes = (/0,1,2/) |
---|
1903 | lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres) |
---|
1904 | |
---|
1905 | amres = True |
---|
1906 | amres@amParallelPosF = 0.65 |
---|
1907 | amres@amOrthogonalPosF = -0.2 |
---|
1908 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
1909 | overlay(plot(n),plot_pt) |
---|
1910 | pt=1 |
---|
1911 | else |
---|
1912 | varn=varn+1 |
---|
1913 | continue |
---|
1914 | end if |
---|
1915 | end if |
---|
1916 | if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then |
---|
1917 | if (q .EQ. 0) then |
---|
1918 | res@gsnLeftString = "q, qv and ql" |
---|
1919 | res@gsnRightString = unit(varn) |
---|
1920 | res@trXMinF = min((/miniq,miniqv,miniql/)) |
---|
1921 | res@trXMaxF = max((/maxiq,maxiqv,maxiql/)) |
---|
1922 | |
---|
1923 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1924 | |
---|
1925 | ; *************************************************** |
---|
1926 | ; legend for combined plot |
---|
1927 | ; *************************************************** |
---|
1928 | |
---|
1929 | lgres = True |
---|
1930 | lgMonoDashIndex = False |
---|
1931 | lgres@lgLabelFont = "helvetica" |
---|
1932 | lgres@lgLabelFontHeightF = .1 |
---|
1933 | lgres@vpWidthF = 0.12 |
---|
1934 | lgres@vpHeightF = 0.1 |
---|
1935 | lgres@lgDashIndexes = (/0,1,2/) |
---|
1936 | lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres) |
---|
1937 | |
---|
1938 | amres = True |
---|
1939 | amres@amParallelPosF = 0.65 |
---|
1940 | amres@amOrthogonalPosF = -0.2 |
---|
1941 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
1942 | overlay(plot(n),plot_q) |
---|
1943 | q=1 |
---|
1944 | else |
---|
1945 | varn=varn+1 |
---|
1946 | continue |
---|
1947 | end if |
---|
1948 | end if |
---|
1949 | if (vNam(varn) .EQ. "rho") then |
---|
1950 | plot(n) = plot_rho |
---|
1951 | end if |
---|
1952 | if (vNam(varn) .EQ. "s") then |
---|
1953 | plot(n) = plot_s |
---|
1954 | end if |
---|
1955 | if (vNam(varn) .EQ. "sa") then |
---|
1956 | plot(n) = plot_sa |
---|
1957 | end if |
---|
1958 | if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then |
---|
1959 | if (e .EQ. 0) then |
---|
1960 | res@gsnLeftString = "e and es" |
---|
1961 | res@gsnRightString = unit(varn) |
---|
1962 | res@trXMinF = min((/minie,minies/)) |
---|
1963 | res@trXMaxF = max((/maxie,maxies/)) |
---|
1964 | |
---|
1965 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1966 | |
---|
1967 | ; *************************************************** |
---|
1968 | ; legend for combined plot |
---|
1969 | ; *************************************************** |
---|
1970 | |
---|
1971 | lgres = True |
---|
1972 | lgMonoDashIndex = False |
---|
1973 | lgres@lgLabelFont = "helvetica" |
---|
1974 | lgres@lgLabelFontHeightF = .1 |
---|
1975 | lgres@vpWidthF = 0.12 |
---|
1976 | lgres@vpHeightF = 0.1 |
---|
1977 | lgres@lgDashIndexes = (/0,1,2/) |
---|
1978 | lbid = gsn_create_legend(wks,3,(/"e","es"/),lgres) |
---|
1979 | |
---|
1980 | amres = True |
---|
1981 | amres@amParallelPosF = 0.65 |
---|
1982 | amres@amOrthogonalPosF = -0.2 |
---|
1983 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
1984 | overlay(plot(n),plot_e) |
---|
1985 | e=1 |
---|
1986 | else |
---|
1987 | varn=varn+1 |
---|
1988 | continue |
---|
1989 | end if |
---|
1990 | end if |
---|
1991 | if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then |
---|
1992 | if (km .EQ. 0) then |
---|
1993 | res@gsnLeftString = "km and kh" |
---|
1994 | res@gsnRightString = unit(varn) |
---|
1995 | res@trXMinF = min((/minikm,minikh/)) |
---|
1996 | res@trXMaxF = max((/maxikm,maxikh/)) |
---|
1997 | |
---|
1998 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) |
---|
1999 | |
---|
2000 | ; *************************************************** |
---|
2001 | ; legend for combined plot |
---|
2002 | ; *************************************************** |
---|
2003 | |
---|
2004 | lgres = True |
---|
2005 | lgMonoDashIndex = False |
---|
2006 | lgres@lgLabelFont = "helvetica" |
---|
2007 | lgres@lgLabelFontHeightF = .1 |
---|
2008 | lgres@vpWidthF = 0.12 |
---|
2009 | lgres@vpHeightF = 0.1 |
---|
2010 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2011 | lbid = gsn_create_legend(wks,3,(/"km","kh"/),lgres) |
---|
2012 | |
---|
2013 | amres = True |
---|
2014 | amres@amParallelPosF = 0.65 |
---|
2015 | amres@amOrthogonalPosF = -0.2 |
---|
2016 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2017 | overlay(plot(n),plot_km) |
---|
2018 | km=1 |
---|
2019 | else |
---|
2020 | varn=varn+1 |
---|
2021 | continue |
---|
2022 | end if |
---|
2023 | end if |
---|
2024 | if (vNam(varn) .EQ. "l") then |
---|
2025 | plot(n) = plot_l |
---|
2026 | end if |
---|
2027 | if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then |
---|
2028 | if (wpup .EQ. 0) then |
---|
2029 | res@gsnLeftString = "wpup, wsus and wu" |
---|
2030 | res@gsnRightString = unit(varn) |
---|
2031 | res@trXMinF = min((/miniwpup,miniwsus,miniwu/)) |
---|
2032 | res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/)) |
---|
2033 | |
---|
2034 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2035 | |
---|
2036 | ; *************************************************** |
---|
2037 | ; legend for combined plot |
---|
2038 | ; *************************************************** |
---|
2039 | |
---|
2040 | lgres = True |
---|
2041 | lgMonoDashIndex = False |
---|
2042 | lgres@lgLabelFont = "helvetica" |
---|
2043 | lgres@lgLabelFontHeightF = .1 |
---|
2044 | lgres@vpWidthF = 0.12 |
---|
2045 | lgres@vpHeightF = 0.1 |
---|
2046 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2047 | lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres) |
---|
2048 | |
---|
2049 | amres = True |
---|
2050 | amres@amParallelPosF = 0.65 |
---|
2051 | amres@amOrthogonalPosF = -0.2 |
---|
2052 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2053 | overlay(plot(n),plot_wpup) |
---|
2054 | wpup=1 |
---|
2055 | else |
---|
2056 | varn=varn+1 |
---|
2057 | continue |
---|
2058 | end if |
---|
2059 | end if |
---|
2060 | if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then |
---|
2061 | if (wpvp .EQ. 0) then |
---|
2062 | res@gsnLeftString = "wpvp, wsus and wv" |
---|
2063 | res@gsnRightString = unit(varn) |
---|
2064 | res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/)) |
---|
2065 | res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/)) |
---|
2066 | |
---|
2067 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2068 | |
---|
2069 | ; *************************************************** |
---|
2070 | ; legend for combined plot |
---|
2071 | ; *************************************************** |
---|
2072 | |
---|
2073 | lgres = True |
---|
2074 | lgMonoDashIndex = False |
---|
2075 | lgres@lgLabelFont = "helvetica" |
---|
2076 | lgres@lgLabelFontHeightF = .1 |
---|
2077 | lgres@vpWidthF = 0.12 |
---|
2078 | lgres@vpHeightF = 0.1 |
---|
2079 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2080 | lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres) |
---|
2081 | |
---|
2082 | amres = True |
---|
2083 | amres@amParallelPosF = 0.65 |
---|
2084 | amres@amOrthogonalPosF = -0.2 |
---|
2085 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2086 | overlay(plot(n),plot_wpvp) |
---|
2087 | wpup=0 |
---|
2088 | else |
---|
2089 | varn=varn+1 |
---|
2090 | continue |
---|
2091 | end if |
---|
2092 | end if |
---|
2093 | if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then |
---|
2094 | if (wpptp .EQ. 0) then |
---|
2095 | res@gsnLeftString = "wpptp, wspts and wv" |
---|
2096 | res@gsnRightString = unit(varn) |
---|
2097 | res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/)) |
---|
2098 | res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/)) |
---|
2099 | |
---|
2100 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2101 | |
---|
2102 | ; *************************************************** |
---|
2103 | ; legend for combined plot |
---|
2104 | ; *************************************************** |
---|
2105 | |
---|
2106 | lgres = True |
---|
2107 | lgMonoDashIndex = False |
---|
2108 | lgres@lgLabelFont = "helvetica" |
---|
2109 | lgres@lgLabelFontHeightF = .1 |
---|
2110 | lgres@vpWidthF = 0.12 |
---|
2111 | lgres@vpHeightF = 0.1 |
---|
2112 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2113 | lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres) |
---|
2114 | |
---|
2115 | amres = True |
---|
2116 | amres@amParallelPosF = 0.65 |
---|
2117 | amres@amOrthogonalPosF = -0.2 |
---|
2118 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2119 | overlay(plot(n),plot_wpptp) |
---|
2120 | wpptp=1 |
---|
2121 | else |
---|
2122 | varn=varn+1 |
---|
2123 | continue |
---|
2124 | end if |
---|
2125 | end if |
---|
2126 | if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then |
---|
2127 | if (wsptsBC .EQ. 0) then |
---|
2128 | res@gsnLeftString = "wsptsBC and wptBC" |
---|
2129 | res@gsnRightString = unit(varn) |
---|
2130 | res@trXMinF = min((/miniwsptsBC,miniwptBC/)) |
---|
2131 | res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/)) |
---|
2132 | |
---|
2133 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2134 | |
---|
2135 | ; *************************************************** |
---|
2136 | ; legend for combined plot |
---|
2137 | ; *************************************************** |
---|
2138 | |
---|
2139 | lgres = True |
---|
2140 | lgMonoDashIndex = False |
---|
2141 | lgres@lgLabelFont = "helvetica" |
---|
2142 | lgres@lgLabelFontHeightF = .1 |
---|
2143 | lgres@vpWidthF = 0.12 |
---|
2144 | lgres@vpHeightF = 0.1 |
---|
2145 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2146 | lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres) |
---|
2147 | |
---|
2148 | amres = True |
---|
2149 | amres@amParallelPosF = 0.65 |
---|
2150 | amres@amOrthogonalPosF = -0.2 |
---|
2151 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2152 | overlay(plot(n),plot_wsptsBC) |
---|
2153 | wsptsBC=1 |
---|
2154 | else |
---|
2155 | varn=varn+1 |
---|
2156 | continue |
---|
2157 | end if |
---|
2158 | end if |
---|
2159 | if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then |
---|
2160 | if (wpvptp .EQ. 0) then |
---|
2161 | res@gsnLeftString = "wpvptp, wsvpts and wv" |
---|
2162 | res@gsnRightString = unit(varn) |
---|
2163 | res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/)) |
---|
2164 | res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/)) |
---|
2165 | |
---|
2166 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2167 | |
---|
2168 | ; *************************************************** |
---|
2169 | ; legend for combined plot |
---|
2170 | ; *************************************************** |
---|
2171 | |
---|
2172 | lgres = True |
---|
2173 | lgMonoDashIndex = False |
---|
2174 | lgres@lgLabelFont = "helvetica" |
---|
2175 | lgres@lgLabelFontHeightF = .1 |
---|
2176 | lgres@vpWidthF = 0.12 |
---|
2177 | lgres@vpHeightF = 0.1 |
---|
2178 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2179 | lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres) |
---|
2180 | |
---|
2181 | amres = True |
---|
2182 | amres@amParallelPosF = 0.65 |
---|
2183 | amres@amOrthogonalPosF = -0.2 |
---|
2184 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2185 | overlay(plot(n),plot_wpvptp) |
---|
2186 | wpvptp=1 |
---|
2187 | else |
---|
2188 | varn=varn+1 |
---|
2189 | continue |
---|
2190 | end if |
---|
2191 | end if |
---|
2192 | if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then |
---|
2193 | if (wpqp .EQ. 0) then |
---|
2194 | res@gsnLeftString = "wpqp, wsqs and wq" |
---|
2195 | res@gsnRightString = unit(varn) |
---|
2196 | res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/)) |
---|
2197 | res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/)) |
---|
2198 | |
---|
2199 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2200 | |
---|
2201 | ; *************************************************** |
---|
2202 | ; legend for combined plot |
---|
2203 | ; *************************************************** |
---|
2204 | |
---|
2205 | lgres = True |
---|
2206 | lgMonoDashIndex = False |
---|
2207 | lgres@lgLabelFont = "helvetica" |
---|
2208 | lgres@lgLabelFontHeightF = .1 |
---|
2209 | lgres@vpWidthF = 0.12 |
---|
2210 | lgres@vpHeightF = 0.1 |
---|
2211 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2212 | lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres) |
---|
2213 | |
---|
2214 | amres = True |
---|
2215 | amres@amParallelPosF = 0.65 |
---|
2216 | amres@amOrthogonalPosF = -0.2 |
---|
2217 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2218 | overlay(plot(n),plot_wpqp) |
---|
2219 | wpqp=1 |
---|
2220 | else |
---|
2221 | varn=varn+1 |
---|
2222 | continue |
---|
2223 | end if |
---|
2224 | end if |
---|
2225 | if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then |
---|
2226 | if (wpqvp .EQ. 0) then |
---|
2227 | res@gsnLeftString = "wpqvp, wsqvs and wqv" |
---|
2228 | res@gsnRightString = unit(varn) |
---|
2229 | res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/)) |
---|
2230 | res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/)) |
---|
2231 | |
---|
2232 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2233 | |
---|
2234 | ; *************************************************** |
---|
2235 | ; legend for combined plot |
---|
2236 | ; *************************************************** |
---|
2237 | |
---|
2238 | lgres = True |
---|
2239 | lgMonoDashIndex = False |
---|
2240 | lgres@lgLabelFont = "helvetica" |
---|
2241 | lgres@lgLabelFontHeightF = .1 |
---|
2242 | lgres@vpWidthF = 0.12 |
---|
2243 | lgres@vpHeightF = 0.1 |
---|
2244 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2245 | lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres) |
---|
2246 | |
---|
2247 | amres = True |
---|
2248 | amres@amParallelPosF = 0.65 |
---|
2249 | amres@amOrthogonalPosF = -0.2 |
---|
2250 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2251 | overlay(plot(n),plot_wpqvp) |
---|
2252 | wpqvp=1 |
---|
2253 | else |
---|
2254 | varn=varn+1 |
---|
2255 | continue |
---|
2256 | end if |
---|
2257 | end if |
---|
2258 | if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then |
---|
2259 | if (wpsp .EQ. 0) then |
---|
2260 | res@gsnLeftString = "wpsp, wsss and ws" |
---|
2261 | res@gsnRightString = unit(varn) |
---|
2262 | res@trXMinF = min((/miniwpsp,miniwsss,miniws/)) |
---|
2263 | res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/)) |
---|
2264 | |
---|
2265 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2266 | |
---|
2267 | ; *************************************************** |
---|
2268 | ; legend for combined plot |
---|
2269 | ; *************************************************** |
---|
2270 | |
---|
2271 | lgres = True |
---|
2272 | lgMonoDashIndex = False |
---|
2273 | lgres@lgLabelFont = "helvetica" |
---|
2274 | lgres@lgLabelFontHeightF = .1 |
---|
2275 | lgres@vpWidthF = 0.12 |
---|
2276 | lgres@vpHeightF = 0.1 |
---|
2277 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2278 | lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres) |
---|
2279 | |
---|
2280 | amres = True |
---|
2281 | amres@amParallelPosF = 0.65 |
---|
2282 | amres@amOrthogonalPosF = -0.2 |
---|
2283 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2284 | overlay(plot(n),plot_wpsp) |
---|
2285 | wpsp=1 |
---|
2286 | else |
---|
2287 | varn=varn+1 |
---|
2288 | continue |
---|
2289 | end if |
---|
2290 | end if |
---|
2291 | if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then |
---|
2292 | if (wpsap .EQ. 0) then |
---|
2293 | res@gsnLeftString = "wpsap, wssas and wsa" |
---|
2294 | res@gsnRightString = unit(varn) |
---|
2295 | res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/)) |
---|
2296 | res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/)) |
---|
2297 | |
---|
2298 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2299 | |
---|
2300 | ; *************************************************** |
---|
2301 | ; legend for combined plot |
---|
2302 | ; *************************************************** |
---|
2303 | |
---|
2304 | lgres = True |
---|
2305 | lgMonoDashIndex = False |
---|
2306 | lgres@lgLabelFont = "helvetica" |
---|
2307 | lgres@lgLabelFontHeightF = .1 |
---|
2308 | lgres@vpWidthF = 0.12 |
---|
2309 | lgres@vpHeightF = 0.1 |
---|
2310 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2311 | lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres) |
---|
2312 | |
---|
2313 | amres = True |
---|
2314 | amres@amParallelPosF = 0.65 |
---|
2315 | amres@amOrthogonalPosF = -0.2 |
---|
2316 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2317 | overlay(plot(n),plot_wpsap) |
---|
2318 | wpsap=1 |
---|
2319 | else |
---|
2320 | varn=varn+1 |
---|
2321 | continue |
---|
2322 | end if |
---|
2323 | end if |
---|
2324 | if (vNam(varn) .EQ. "wses") then |
---|
2325 | plot(n) = plot_wses |
---|
2326 | end if |
---|
2327 | if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then |
---|
2328 | if (us2 .EQ. 0) then |
---|
2329 | res@gsnLeftString = "us2, vs2 and ws2" |
---|
2330 | res@gsnRightString = unit(varn) |
---|
2331 | res@trXMinF = min((/minius2,minivs2,miniws2/)) |
---|
2332 | res@trXMaxF = max((/maxius2,maxivs2,maxiws2/)) |
---|
2333 | |
---|
2334 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2335 | |
---|
2336 | ; *************************************************** |
---|
2337 | ; legend for combined plot |
---|
2338 | ; *************************************************** |
---|
2339 | |
---|
2340 | lgres = True |
---|
2341 | lgMonoDashIndex = False |
---|
2342 | lgres@lgLabelFont = "helvetica" |
---|
2343 | lgres@lgLabelFontHeightF = .1 |
---|
2344 | lgres@vpWidthF = 0.12 |
---|
2345 | lgres@vpHeightF = 0.1 |
---|
2346 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2347 | lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres) |
---|
2348 | |
---|
2349 | amres = True |
---|
2350 | amres@amParallelPosF = 0.65 |
---|
2351 | amres@amOrthogonalPosF = -0.2 |
---|
2352 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2353 | overlay(plot(n),plot_us2) |
---|
2354 | us2=1 |
---|
2355 | else |
---|
2356 | varn=varn+1 |
---|
2357 | continue |
---|
2358 | end if |
---|
2359 | end if |
---|
2360 | if (vNam(varn) .EQ. "pts2") then |
---|
2361 | plot(n) = plot_pts2 |
---|
2362 | end if |
---|
2363 | if (vNam(varn) .EQ. "ws3") then |
---|
2364 | plot(n) = plot_ws3 |
---|
2365 | end if |
---|
2366 | if (vNam(varn) .EQ. "Sw") then |
---|
2367 | plot(n) = plot_Sw |
---|
2368 | end if |
---|
2369 | if (vNam(varn) .EQ. "ws2pts") then |
---|
2370 | plot(n) = plot_ws2pts |
---|
2371 | end if |
---|
2372 | if (vNam(varn) .EQ. "wspts2") then |
---|
2373 | plot(n) = plot_wspts2 |
---|
2374 | end if |
---|
2375 | if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then |
---|
2376 | if (wsususodz .EQ. 0) then |
---|
2377 | res@gsnLeftString = "wsususodz, wspsodz and ws2" |
---|
2378 | res@gsnRightString = unit(varn) |
---|
2379 | res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/)) |
---|
2380 | res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/)) |
---|
2381 | |
---|
2382 | plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) |
---|
2383 | |
---|
2384 | ; *************************************************** |
---|
2385 | ; legend for combined plot |
---|
2386 | ; *************************************************** |
---|
2387 | |
---|
2388 | lgres = True |
---|
2389 | lgMonoDashIndex = False |
---|
2390 | lgres@lgLabelFont = "helvetica" |
---|
2391 | lgres@lgLabelFontHeightF = .1 |
---|
2392 | lgres@vpWidthF = 0.12 |
---|
2393 | lgres@vpHeightF = 0.1 |
---|
2394 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2395 | lbid = gsn_create_legend(wks,3,(/"wsususodz","wspsodz","wpeodz"/),lgres) |
---|
2396 | |
---|
2397 | amres = True |
---|
2398 | amres@amParallelPosF = 0.65 |
---|
2399 | amres@amOrthogonalPosF = -0.2 |
---|
2400 | annoid1 = gsn_add_annotation(plot(n),lbid,amres) |
---|
2401 | overlay(plot(n),plot_wsususodz) |
---|
2402 | wsususodz=1 |
---|
2403 | else |
---|
2404 | varn=varn+1 |
---|
2405 | continue |
---|
2406 | end if |
---|
2407 | end if |
---|
2408 | n=n+1 |
---|
2409 | varn=varn+1 |
---|
2410 | end if |
---|
2411 | end do |
---|
2412 | end if |
---|
2413 | |
---|
2414 | if (combine .EQ. 1) then |
---|
2415 | |
---|
2416 | check=True |
---|
2417 | n_o=0 |
---|
2418 | do varn = 0,dim-1 |
---|
2419 | |
---|
2420 | if ( isStrSubset( vNam(varn), "NORM") .or. isStrSubset( vNam(varn), "time") ) |
---|
2421 | check = False |
---|
2422 | else |
---|
2423 | if (.not. isvar("var")) then |
---|
2424 | check = True |
---|
2425 | if (parameter(21) .NE. "variables") then |
---|
2426 | var=parameter(21) |
---|
2427 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
2428 | end if |
---|
2429 | else |
---|
2430 | check = isStrSubset( var,","+vNam(varn)+"," ) |
---|
2431 | end if |
---|
2432 | end if |
---|
2433 | |
---|
2434 | if(check) then |
---|
2435 | |
---|
2436 | z = f->$vNam(varn+1)$ |
---|
2437 | com=isStrSubset( c_var,","+vNam(varn)+"," ) |
---|
2438 | |
---|
2439 | if (com) |
---|
2440 | |
---|
2441 | if (n_o .EQ. 1) then |
---|
2442 | res@xyDashPattern = 1 |
---|
2443 | |
---|
2444 | else |
---|
2445 | if (n_o .EQ. 2) then |
---|
2446 | res@xyDashPattern = 2 |
---|
2447 | else |
---|
2448 | res@xyDashPattern = 0 |
---|
2449 | res@gsnLeftString = c_var |
---|
2450 | res@gsnRightString = unit(varn) |
---|
2451 | res@trXMinF = min(mini) |
---|
2452 | res@trXMaxF = max(maxi) |
---|
2453 | end if |
---|
2454 | end if |
---|
2455 | label(n_o)=vNam(varn) |
---|
2456 | color_o(n_o)=237 |
---|
2457 | plot_o(n_o)=gsn_csm_xy(wks,data(varn,:,:),z,res) |
---|
2458 | n_o=n_o+1 |
---|
2459 | end if |
---|
2460 | varn=varn+1 |
---|
2461 | end if |
---|
2462 | end do |
---|
2463 | |
---|
2464 | if(number_comb .EQ. 2)then |
---|
2465 | overlay(plot_o(0),plot_o(1)) |
---|
2466 | end if |
---|
2467 | if(number_comb .EQ. 3)then |
---|
2468 | overlay(plot_o(0),plot_o(1)) |
---|
2469 | overlay(plot_o(0),plot_o(2)) |
---|
2470 | end if |
---|
2471 | |
---|
2472 | ; *************************************************** |
---|
2473 | ; legend for combined plot |
---|
2474 | ; *************************************************** |
---|
2475 | |
---|
2476 | lgres = True |
---|
2477 | lgMonoDashIndex = False |
---|
2478 | lgres@lgDashIndexes = (/0,1,2/) |
---|
2479 | lgres@lgLabelFont = "helvetica" |
---|
2480 | lgres@lgLabelFontHeightF = .1 |
---|
2481 | lgres@vpWidthF = 0.12 |
---|
2482 | lgres@vpHeightF = 0.1 |
---|
2483 | |
---|
2484 | lbid = gsn_create_legend(wks,number_comb,label,lgres) |
---|
2485 | |
---|
2486 | amres = True |
---|
2487 | amres@amParallelPosF = 0.65 |
---|
2488 | amres@amOrthogonalPosF = -0.2 |
---|
2489 | annoid1 = gsn_add_annotation(plot_o(0),lbid,amres) |
---|
2490 | |
---|
2491 | plot(0) = plot_o(0) |
---|
2492 | |
---|
2493 | end if |
---|
2494 | |
---|
2495 | ; *************************************************** |
---|
2496 | ; merge plots onto one page |
---|
2497 | ; *************************************************** |
---|
2498 | print(n) |
---|
2499 | do m=0,n-1 |
---|
2500 | plot_(m)=plot(n-1-m) |
---|
2501 | end do |
---|
2502 | |
---|
2503 | if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then |
---|
2504 | gsn_panel(wks,plot_(0:n-1),(/n,1/),resP) |
---|
2505 | else |
---|
2506 | do i = 0,n-1, no_lines*no_columns |
---|
2507 | if( (i+no_lines*no_columns) .gt. (n-1)) then |
---|
2508 | gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP) |
---|
2509 | else |
---|
2510 | gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP) |
---|
2511 | end if |
---|
2512 | end do |
---|
2513 | end if |
---|
2514 | |
---|
2515 | print(" ") |
---|
2516 | print("Output to: " + file_out +"."+ format_out) |
---|
2517 | print(" ") |
---|
2518 | |
---|
2519 | end |
---|