Home

Context Navigation

profiles.ncl @ 174

Last change on this file since 174 was 174, checked in by letzel, 17 years ago
NCL scripts in trunk/SCRIPTS/NCL updated
File size: 144.9 KB

Line
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",77,"string")
13	delete(parameter@_FillValue)
14	else
15	print(" ")
16	print("'.ncl_preferences' is not existent")
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(37)) .NE. 0) then
110	over = stringtointeger(parameter(37))
111	if (stringtointeger(parameter(37)) .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	if (.not. isvar("prof3d"))then
121	prof3d = 0
122	if (stringtointeger(parameter(75)) .NE. 0) then
123	prof3d = stringtointeger(parameter(75))
124	if (stringtointeger(parameter(75)) .NE. 1) then
125	print(" ")
126	print("Please set 'prof3d' to 0 or 1")
127	print(" ")
128	exit
129	end if
130	end if
131	end if
132
133	; ***************************************************
134	; open input file
135	; ***************************************************
136
137	f=addfile( file_in,"r")
138
139	vNam = getfilevarnames(f)
140	print(" ")
141	print("Variable on netCDF file: " + vNam)
142	print(" ")
143	dim = dimsizes(vNam)
144
145	if (dim .EQ. 0) then
146	print(" ")
147	print("There are no data on file")
148	print(" ")
149	end if
150
151	if (prof3d .EQ. 0) then
152
153	do varn = dim-1,0,1
154	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
155	varn=varn-1
156	continue
157	end if
158	if (vNam(varn) .EQ. "u" .OR. isStrSubset(vNam(varn), "u_"))
159	z_u = f->$vNam(varn+1)$
160	break
161	else
162	if (vNam(varn) .EQ. "v" .OR. isStrSubset(vNam(varn), "v_"))then
163
164	z_u = f->$vNam(varn+1)$
165	break
166	else
167	if(vNam(varn) .EQ. "pt" .OR. isStrSubset(vNam(varn), "pt_"))then
168	z_u = f->$vNam(varn+1)$
169	break
170	else
171	if(vNam(varn) .EQ. "vpt" .OR. isStrSubset(vNam(varn), "vpt_"))then
172	z_u = f->$vNam(varn+1)$
173	break
174	else
175	if(vNam(varn) .EQ. "lpt" .OR. isStrSubset(vNam(varn), "lpt_"))then
176	z_u = f->$vNam(varn+1)$
177	break
178	else
179	if(vNam(varn) .EQ. "q" .OR. isStrSubset(vNam(varn), "q_") )then
180	z_u = f->$vNam(varn+1)$
181	break
182	else
183	if(vNam(varn) .EQ. "qv" .OR. isStrSubset(vNam(varn), "qv_"))then
184	z_u = f->$vNam(varn+1)$
185	break
186	else
187	if(vNam(varn) .EQ. "ql" .OR. isStrSubset(vNam(varn), "ql_"))then
188	z_u = f->$vNam(varn+1)$
189	break
190	else
191	if(vNam(varn) .EQ. "rho" .OR. isStrSubset(vNam(varn), "rho_"))then
192	z_u = f->$vNam(varn+1)$
193	break
194	else
195	if(vNam(varn) .EQ. "s" .OR. isStrSubset(vNam(varn), "s_") )then
196	z_u = f->$vNam(varn+1)$
197	break
198	else
199	if(vNam(varn) .EQ. "sa" .OR. isStrSubset(vNam(varn), "sa_"))then
200	z_u = f->$vNam(varn+1)$
201	break
202	else
203	if(vNam(varn) .EQ. "e" .OR. isStrSubset(vNam(varn), "e_"))then
204	z_u = f->$vNam(varn+1)$
205	break
206	else
207	if(vNam(varn) .EQ. "es" .OR. isStrSubset(vNam(varn), "es_"))then
208	z_u = f->$vNam(varn+1)$
209	break
210	else
211	if(vNam(varn) .EQ. "km" .OR. isStrSubset(vNam(varn), "km_"))then
212	z_u = f->$vNam(varn+1)$
213	break
214	else
215	if(vNam(varn) .EQ. "kh" .OR. isStrSubset(vNam(varn), "kh_"))then
216	z_u = f->$vNam(varn+1)$
217	break
218	else
219	if(vNam(varn) .EQ. "l" .OR. isStrSubset(vNam(varn), "l_"))then
220	z_u = f->$vNam(varn+1)$
221	break
222	else
223	if(vNam(varn) .EQ. "us2" .OR. isStrSubset(vNam(varn), "us2_"))then
224	z_u = f->$vNam(varn+1)$
225	break
226	else
227	if(vNam(varn) .EQ. "vs2" .OR. isStrSubset(vNam(varn), "vs2_") )then
228	z_u = f->$vNam(varn+1)$
229	break
230	else
231	if(vNam(varn) .EQ. "pts2" .OR. isStrSubset(vNam(varn), "pts2_"))then
232	z_u = f->$vNam(varn+1)$
233	break
234	else
235	if(vNam(varn) .EQ. "wsususodz" .OR. isStrSubset(vNam(varn), "wsususodz_"))then
236	z_u = f->$vNam(varn+1)$
237	break
238	else
239	if(vNam(varn) .EQ. "wspsodz" .OR. isStrSubset(vNam(varn), "wspsodz_"))then
240	z_u = f->$vNam(varn+1)$
241	break
242	else
243	if(vNam(varn) .EQ. "wpeodz" .OR. isStrSubset(vNam(varn), "wpeodz_"))then
244	z_u = f->$vNam(varn+1)$
245	break
246	end if
247	end if
248	end if
249	end if
250	end if
251	end if
252	end if
253	end if
254	end if
255	end if
256	end if
257	end if
258	end if
259	end if
260	end if
261	end if
262	end if
263	end if
264	end if
265	end if
266	end if
267	end if
268	varn=varn-1
269	end do
270
271	do varn=dim-1,0,1
272	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
273	varn=varn-1
274	continue
275	end if
276	if (vNam(varn) .EQ. "w" .OR. isStrSubset(vNam(varn), "w_"))
277	z_w = f->$vNam(varn+1)$
278	break
279	else
280	if (vNam(varn) .EQ. "wpup" .OR. isStrSubset(vNam(varn), "wpup_"))then
281	z_w = f->$vNam(varn+1)$
282	break
283	else
284	if(vNam(varn) .EQ. "wsus" .OR. isStrSubset(vNam(varn), "wsus_"))then
285	z_w = f->$vNam(varn+1)$
286	break
287	else
288	if(vNam(varn) .EQ. "wu" .OR. isStrSubset(vNam(varn), "wu_"))then
289	z_w = f->$vNam(varn+1)$
290	break
291	else
292	if(vNam(varn) .EQ. "wpvp" .OR. isStrSubset(vNam(varn), "wpvp_"))then
293	z_w = f->$vNam(varn+1)$
294	break
295	else
296	if(vNam(varn) .EQ. "wsvs" .OR. isStrSubset(vNam(varn), "wsvs_"))then
297	z_w = f->$vNam(varn+1)$
298	break
299	else
300	if(vNam(varn) .EQ. "wv" .OR. isStrSubset(vNam(varn), "wv_"))then
301	z_w = f->$vNam(varn+1)$
302	break
303	else
304	if(vNam(varn) .EQ. "wptpp" .OR. isStrSubset(vNam(varn), "wptpp_"))then
305	z_w = f->$vNam(varn+1)$
306	break
307	else
308	if(vNam(varn) .EQ. "wspts" .OR. isStrSubset(vNam(varn), "wspts_"))then
309	z_w = f->$vNam(varn+1)$
310	break
311	else
312	if(vNam(varn) .EQ. "wpt" .OR. isStrSubset(vNam(varn), "wpz_"))then
313	z_w = f->$vNam(varn+1)$
314	break
315	else
316	if(vNam(varn) .EQ. "wsptsBC" .OR. isStrSubset(vNam(varn), "wsptsBC_"))then
317	z_w = f->$vNam(varn+1)$
318	break
319	else
320	if(vNam(varn) .EQ. "wptBC" .OR. isStrSubset(vNam(varn), "wptBC_"))then
321	z_w = f->$vNam(varn+1)$
322	break
323	else
324	if(vNam(varn) .EQ. "wpvptp" .OR. isStrSubset(vNam(varn), "wpvptp_"))then
325	z_w = f->$vNam(varn+1)$
326	break
327	else
328	if(vNam(varn) .EQ. "wsvpts" .OR. isStrSubset(vNam(varn), "wsvpts_"))then
329	z_w = f->$vNam(varn+1)$
330	break
331	else
332	if(vNam(varn) .EQ. "wvpt" .OR. isStrSubset(vNam(varn), "wvpt_"))then
333	z_w = f->$vNam(varn+1)$
334	break
335	else
336	if(vNam(varn) .EQ. "wpqp" .OR. isStrSubset(vNam(varn), "wpqp_"))then
337	z_w = f->$vNam(varn+1)$
338	break
339	else
340	if(vNam(varn) .EQ. "wsqs" .OR. isStrSubset(vNam(varn), "wsqs_"))then
341	z_w = f->$vNam(varn+1)$
342	break
343	else
344	if(vNam(varn) .EQ. "wq" .OR. isStrSubset(vNam(varn), "wq_"))then
345	z_w = f->$vNam(varn+1)$
346	break
347	else
348	if(vNam(varn) .EQ. "wpqvp" .OR. isStrSubset(vNam(varn), "wpqvp_"))then
349	z_w = f->$vNam(varn+1)$
350	break
351	else
352	if(vNam(varn) .EQ. "wsqvs" .OR. isStrSubset(vNam(varn), "wsqvs_"))then
353	z_w = f->$vNam(varn+1)$
354	break
355	else
356	if(vNam(varn) .EQ. "wqv" .OR. isStrSubset(vNam(varn), "wqv_"))then
357	z_w = f->$vNam(varn+1)$
358	break
359	else
360	if(vNam(varn) .EQ. "wpsp" .OR. isStrSubset(vNam(varn), "wpsp_"))then
361	z_w = f->$vNam(varn+1)$
362	break
363	else
364	if(vNam(varn) .EQ. "wsss" .OR. isStrSubset(vNam(varn), "wsss_"))then
365	z_w = f->$vNam(varn+1)$
366	break
367	else
368	if(vNam(varn) .EQ. "ws" .OR. isStrSubset(vNam(varn), "ws_"))then
369	z_w = f->$vNam(varn+1)$
370	break
371	else
372	if(vNam(varn) .EQ. "wpsap" .OR. isStrSubset(vNam(varn), "wpsap_"))then
373	z_w = f->$vNam(varn+1)$
374	break
375	else
376	if(vNam(varn) .EQ. "wssas" .OR. isStrSubset(vNam(varn), "wssas_") )then
377	z_w = f->$vNam(varn+1)$
378	break
379	else
380	if(vNam(varn) .EQ. "wsa" .OR. isStrSubset(vNam(varn), "wsa_"))then
381	z_w = f->$vNam(varn+1)$
382	break
383	else
384	if(vNam(varn) .EQ. "wses" .OR. isStrSubset(vNam(varn), "wses_"))then
385	z_w = f->$vNam(varn+1)$
386	break
387	else
388	if(vNam(varn) .EQ. "ws2" .OR. isStrSubset(vNam(varn), "ws2_"))then
389	z_w = f->$vNam(varn+1)$
390	break
391	else
392	if(vNam(varn) .EQ. "ws3" .OR. isStrSubset(vNam(varn), "ws3_"))then
393	z_w = f->$vNam(varn+1)$
394	break
395	else
396	if(vNam(varn) .EQ. "Sw" .OR. isStrSubset(vNam(varn), "Sw_"))then
397	z_w = f->$vNam(varn+1)$
398	break
399	else
400	if(vNam(varn) .EQ. "ws2pts".OR. isStrSubset(vNam(varn), "ws2pts_") )then
401	z_w = f->$vNam(varn+1)$
402	break
403	else
404	if(vNam(varn) .EQ. "wspts2" .OR. isStrSubset(vNam(varn), "wspts2_"))then
405	z_w = f->$vNam(varn+1)$
406	break
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	end if
424	end if
425	end if
426	end if
427	end if
428	end if
429	end if
430	end if
431	end if
432	end if
433	end if
434	end if
435	end if
436	end if
437	end if
438	end if
439	end if
440	varn=varn-1
441	end do
442	if ( .not. isvar("z_u") ) then
443	zu = 0
444	else
445	zu = 1
446	end if
447	if ( .not. isvar("z_w") ) then
448	zw = 0
449	else
450	zw = 1
451	end if
452
453	if (zu .EQ. 0 .AND. zw .EQ. 0) then
454	do varn=0,dim-1
455	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
456	check = False
457	else
458	if (.not. isvar("var")) then
459	check = True
460	if (parameter(21) .NE. "variables") then
461	var=parameter(21)
462	check = isStrSubset( var,","+vNam(varn)+"," )
463	end if
464	else
465	check = isStrSubset( var,","+vNam(varn)+"," )
466	end if
467	end if
468	if (check)
469	z = f->$vNam(varn+1)$
470	dimz = dimsizes(z)
471	break
472	end if
473	end do
474	end if
475	if (isvar("z_u") ) then
476	dimz = dimsizes(z_u)
477	else
478	if (isvar("z_w"))then
479	dimz = dimsizes(z_w)
480	end if
481	end if
482
483	else
484
485	do varn = dim-1,0,1
486	if (vNam(varn) .EQ. "zu_3d")
487	z_u = f->zu_3d
488	dimz = dimsizes(z_u)
489	else
490	if (vNam(varn) .EQ. "zw_3d")
491	z_w = f->zw_3d
492	dimz = dimsizes(z_w)
493	end if
494	end if
495	end do
496
497	end if
498
499	t_all = f->time
500	nt = dimsizes(t_all)
501	delta_t=t_all(nt-1)/nt
502
503	; ****************************************************
504	; start of time step and different types of mistakes that could be done
505	; ****************************************************
506
507	if ( .not. isvar("start_time_step") ) then
508	start_time_step=t_all(0)/3600
509	if (parameter(13) .NE. "t(0)") then
510	if (stringtodouble(parameter(13)) .GT. t_all(nt-1)/3600)
511	print(" ")
512	print("'start_time_step' = "+ parameter(13) +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
513	print(" ")
514	print("Please select another 'start_time_step'")
515	print(" ")
516	exit
517	end if
518	if (stringtofloat(parameter(13)) .LT. t_all(0)/3600)
519	print(" ")
520	print("'start_time_step' = "+ parameter(13) +"h is lower than first time step = " + t_all(0)+"s = "+t_all(0)/3600+"h")
521	exit
522	end if
523	start_time_step=stringtodouble(parameter(13))
524	end if
525	else
526	if (start_time_step .GT. t_all(nt-1)/3600)
527	print(" ")
528	print("'start_time_step' = "+ start_time_step +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
529	print(" ")
530	print("Please select another 'start_time_step'")
531	print(" ")
532	exit
533	end if
534	if (start_time_step .LT. t_all(0)/3600)
535	print(" ")
536	print("'start_time_step' = "+ start_time_step +"h is lower than first time step = " + t_all(0)+"s = "+t_all(0)/3600+"h")
537
538	exit
539	end if
540	end if
541	start_time_step = start_time_step*3600
542	;if (start_time_step .GE. t_all(0) .AND. start_time_step .LT. t_all(1)+delta_t/2)then
543	; st=0
544	;end if
545	do i=0,nt-1
546	if (start_time_step .GE. t_all(i)-delta_t/2 .AND. start_time_step .LT. t_all(i)+delta_t/2)then
547	st=i
548	break
549	end if
550	end do
551
552	; ****************************************************
553	; end of time step and different types of mistakes that could be done
554	; ****************************************************
555
556	if ( .not. isvar("end_time_step") ) then
557	end_time_step = t_all(nt-1)/3600
558	if (parameter(15) .NE. "t(end)") then
559	if (stringtodouble(parameter(15)) .GT. t_all(nt-1)/3600)
560	print(" ")
561	print("'end_time_step' = "+ parameter(15) +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
562	print(" ")
563	print("Please select another 'end_time_step'")
564	print(" ")
565	exit
566	end if
567	if (stringtodouble(parameter(15)) .LT. start_time_step/3600)
568	print(" ")
569	print("'end_time_step' = "+ parameter(15) +"h is lower than 'start_time_step' = "+start_time_step/3600+"h")
570	print(" ")
571	print("Please select another 'start_time_step' or 'end_time_step'")
572	print(" ")
573	exit
574	end if
575	end_time_step = stringtodouble(parameter(15))
576	end if
577	else
578	if (end_time_step .GT. t_all(nt-1)/3600)
579	print(" ")
580	print("'end_time_step' = "+ end_time_step +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
581	print(" ")
582	print("Please select another 'end_time_step'")
583	print(" ")
584	exit
585	end if
586	if (end_time_step .LT. start_time_step/3600)
587	print(" ")
588	print("'end_time_step' = "+ end_time_step +"h is lower than 'start_time_step' = "+start_time_step+"h")
589	print(" ")
590	print("Please select another 'start_time_step' or 'end_time_step'")
591	print(" ")
592	exit
593	end if
594	end if
595	end_time_step = end_time_step*3600
596	;if (end_time_step .GE. t_all(0) .AND. end_time_step .LT. t_all(1)+delta_t/2)then
597	; et=1
598	;end if
599	do i=0,nt-1
600	if (end_time_step .GE. t_all(i)-delta_t/2 .AND. end_time_step .LT. t_all(i)+delta_t/2)then
601	et=i
602	break
603	end if
604	end do
605
606	delete(start_time_step)
607	start_time_step=round(st,3)
608	delete(end_time_step)
609	end_time_step=round(et,3)
610
611	; ****************************************************
612	; set up legend and colors
613	; ****************************************************
614
615	legend_label=new(nt,double)
616	do p=start_time_step,end_time_step
617	legend_label(p-start_time_step)=t_all(p)
618	end do
619
620	np = end_time_step-start_time_step+1
621	if ( black .eq. 0 ) then
622	res=True
623	res@xyLineColors = ispan(2,237,235/np)
624	end if
625
626	; ***************************************************
627	; set up recourses
628	; ***************************************************
629
630	res@gsnDraw = False
631	res@gsnFrame = False
632	res@gsnPaperOrientation = "portrait"
633	res@gsnPaperWidth = 8.27
634	res@gsnPaperHeight = 11.69
635	res@gsnPaperMargin = 0.79
636	res@txFont = "helvetica"
637	res@tiMainFont = "helvetica"
638	res@tiXAxisFont = "helvetica"
639	res@tiYAxisFont = "helvetica"
640	res@tmXBLabelFont = "helvetica"
641	res@tmYLLabelFont = "helvetica"
642	res@lgLabelFont = "helvetica"
643	res@tmLabelAutoStride = True
644	res@pmLegendDisplayMode = "Always"
645	res@pmLegendSide = "Top"
646	res@xyExplicitLegendLabels = legend_label
647	res@pmLegendParallelPosF = 1.15
648	res@pmLegendOrthogonalPosF = -1.0
649	res@pmLegendWidthF = 0.12
650	res@pmLegendHeightF = 0.3
651	res@lgLabelFontHeightF = .02
652	res@txFontHeightF = 0.02
653	res@tiXAxisFontHeightF = 0.02
654	res@tiYAxisFontHeightF = 0.02
655	res@tiXAxisString = " "
656	res@tiYAxisString = "Height [z]"
657
658	if ( dash .eq. 0 ) then
659	res@xyMonoDashPattern = True
660	end if
661
662	resP = True
663	resP@txFont = "helvetica"
664	resP@txString = f@title
665	resP@txFuncCode = "~"
666	resP@txFontHeightF = 0.014
667
668	; ***************************************************
669	; set up graphics for plot
670	; ***************************************************
671
672	plot = new(dim,graphic)
673	plot_ = new(dim,graphic)
674
675	if (combine .EQ. 1) then
676	if ( .not. isvar("number_comb") ) then
677	if (parameter(25) .EQ. "0") then
678	print(" ")
679	print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot")
680	print(" ")
681	exit
682	else
683	number_comb=stringtointeger(parameter(25))
684	plot_o = new(number_comb,graphic)
685	end if
686	else
687	if(number_comb .EQ. 2 .OR. number_comb .EQ. 3) then
688	plot_o = new(number_comb,graphic)
689	else
690	print(" ")
691	print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot")
692	print(" ")
693	exit
694	end if
695	end if
696	label=new(number_comb,string)
697	color_o=new(number_comb,integer)
698	mini=new(number_comb,float)
699	maxi=new(number_comb,float)
700	end if
701
702	wks=gsn_open_wks(format_out,file_out)
703	gsn_define_colormap(wks,"rainbow+white")
704
705	; ***************************************************
706	; indicate plot number
707	; ***************************************************
708
709	if (combine .EQ. 1) then
710	n = 1
711	else
712	n = 0
713	end if
714
715	; ***************************************************
716	; set up minimum and maximum height
717	; ***************************************************
718
719	if (.not. isvar("min_z"))
720	min_z=0
721	if (stringtointeger(parameter(33)) .NE. 0) then
722	if (stringtointeger(parameter(33)) .GE. max(z_u) ) then
723	print(" ")
724	print("Minimum of height ('min_z'="+stringtointeger(parameter(33))+") is greater than available heights (="+max(z_u)+")")
725	print(" ")
726	exit
727	end if
728	if (stringtointeger(parameter(33)) .LT. 0 ) then
729	print(" ")
730	print("Begin minimum of height 'min_z' with 0")
731	print(" ")
732	exit
733	end if
734	min_z=stringtointeger(parameter(33))
735	end if
736	else
737	if (min_z .GE. max(z_u) ) then
738	print(" ")
739	print("Minimum of height ('min_z'="+min_z+") is greater than available heights (="+max(z_u)+")")
740	print(" ")
741	exit
742	end if
743	if (min_z .LT. 0 ) then
744	print(" ")
745	print("Begin minimum of height 'min_z' with 0")
746	print(" ")
747	exit
748	end if
749	end if
750
751	if (.not. isvar("max_z"))
752	max_z=max(z_u)
753	if ((parameter(35)) .NE. "max(z_u)") then
754	if (stringtofloat(parameter(35)) .GE. max(z_u) ) then
755	print(" ")
756	print("Maximum of height ('max_z'="+parameter(35)+") is greater than available heights (="+max(z_u)+")")
757	print(" ")
758	exit
759	end if
760	if (stringtointeger(parameter(35)) .LE. 0 ) then
761	print(" ")
762	print("Maximum of height 'max_z' should be at least 1")
763	print(" ")
764	exit
765	end if
766	max_z=stringtointeger(parameter(35))
767	end if
768	else
769	if (max_z .GE. max(z_pr) ) then
770	print(" ")
771	print("Maximum of height ('max_z'="+max_z+") is greater than available heights (="+max(z_u)+")")
772	print(" ")
773	exit
774	end if
775	if (max_z .LE. 0 ) then
776	print(" ")
777	print("Maximum of height 'max_z' should be at least 1")
778	print(" ")
779	exit
780	end if
781	end if
782
783	; ***************************************************
784	; read data and create plots
785	; ***************************************************
786
787	do ti = start_time_step, end_time_step
788	if( t_all(ti) .lt. 10^36) then
789	start_time_step = ti
790	break
791	end if
792	end do
793
794	data = new((/dim,(end_time_step-start_time_step)+1,dimz/),float)
795	data_0 = new((/(end_time_step-start_time_step)+1,dimz/),float)
796	data_0 = 0.0
797	t = new((/(end_time_step-start_time_step)+1,dimz/),float)
798	t = 0.0
799	unit = new(dim,string)
800	z_v = new((/dim,dimz/),double)
801
802	if (over .EQ. 1) then
803	plot_u = gsn_csm_xy(wks,t,data_0(:,:),res)
804	miniu = 1.E27
805	maxiu =-1.E27
806	plot_v = gsn_csm_xy(wks,t,data_0(:,:),res)
807	miniv = 1.E27
808	maxiv =-1.E27
809	plot_w = gsn_csm_xy(wks,t,data_0(:,:),res)
810	miniw = 1.E27
811	maxiw =-1.E27
812	plot_pt = gsn_csm_xy(wks,t,data_0(:,:),res)
813	minipt = 1.E27
814	maxipt =-1.E27
815	plot_vpt = gsn_csm_xy(wks,t,data_0(:,:),res)
816	minivpt = 1.E27
817	maxivpt =-1.E27
818	plot_lpt = gsn_csm_xy(wks,t,data_0(:,:),res)
819	minilpt = 1.E27
820	maxilpt =-1.E27
821	plot_q = gsn_csm_xy(wks,t,data_0(:,:),res)
822	miniq = 1.E27
823	maxiq =-1.E27
824	plot_qv = gsn_csm_xy(wks,t,data_0(:,:),res)
825	miniqv = 1.E27
826	maxiqv =-1.E27
827	plot_ql = gsn_csm_xy(wks,t,data_0(:,:),res)
828	miniql = 1.E27
829	maxiql =-1.E27
830	plot_rho = gsn_csm_xy(wks,t,data_0(:,:),res)
831	plot_s = gsn_csm_xy(wks,t,data_0(:,:),res)
832	plot_sa = gsn_csm_xy(wks,t,data_0(:,:),res)
833	plot_e = gsn_csm_xy(wks,t,data_0(:,:),res)
834	minie = 1.E27
835	maxie =-1.E27
836	plot_es = gsn_csm_xy(wks,t,data_0(:,:),res)
837	minie = 1.E27
838	maxie =-1.E27
839	plot_km = gsn_csm_xy(wks,t,data_0(:,:),res)
840	minie = 1.E27
841	maxie =-1.E27
842	plot_kh = gsn_csm_xy(wks,t,data_0(:,:),res)
843	minie = 1.E27
844	maxie =-1.E27
845	plot_l = gsn_csm_xy(wks,t,data_0(:,:),res)
846	plot_wpup = gsn_csm_xy(wks,t,data_0(:,:),res)
847	miniwpup = 1.E27
848	maxiwpup =-1.E27
849	plot_wsus = gsn_csm_xy(wks,t,data_0(:,:),res)
850	miniwsus = 1.E27
851	maxiwsus =-1.E27
852	plot_wu = gsn_csm_xy(wks,t,data_0(:,:),res)
853	miniwu = 1.E27
854	maxiwu =-1.E27
855	plot_wpvp = gsn_csm_xy(wks,t,data_0(:,:),res)
856	miniwpvp = 1.E27
857	maxiwpvp =-1.E27
858	plot_wsvs = gsn_csm_xy(wks,t,data_0(:,:),res)
859	miniwsvs = 1.E27
860	maxiwsvs =-1.E27
861	plot_wv = gsn_csm_xy(wks,t,data_0(:,:),res)
862	miniwv = 1.E27
863	maxiwv =-1.E27
864	plot_wpptp = gsn_csm_xy(wks,t,data_0(:,:),res)
865	miniwpptp = 1.E27
866	maxiwpptp =-1.E27
867	plot_wspts = gsn_csm_xy(wks,t,data_0(:,:),res)
868	miniwspts = 1.E27
869	maxiwspts =-1.E27
870	plot_wpt = gsn_csm_xy(wks,t,data_0(:,:),res)
871	miniwpt = 1.E27
872	maxiwpt =-1.E27
873	plot_wsptsBC = gsn_csm_xy(wks,t,data_0(:,:),res)
874	miniwsptsBC = 1.E27
875	maxiwsptsBC =-1.E27
876	plot_wptBC = gsn_csm_xy(wks,t,data_0(:,:),res)
877	miniwptBC = 1.E27
878	maxiwptBC =-1.E27
879	plot_wpvptp = gsn_csm_xy(wks,t,data_0(:,:),res)
880	miniwpvptp = 1.E27
881	maxiwpvptp =-1.E27
882	plot_wsvpts = gsn_csm_xy(wks,t,data_0(:,:),res)
883	miniwsvpts = 1.E27
884	maxiewsvpts=-1.E27
885	plot_wvpt = gsn_csm_xy(wks,t,data_0(:,:),res)
886	miniwvpt = 1.E27
887	maxiwvpt =-1.E27
888	plot_wpqp = gsn_csm_xy(wks,t,data_0(:,:),res)
889	miniwpqp = 1.E27
890	maxiwpqp =-1.E27
891	plot_wsqs = gsn_csm_xy(wks,t,data_0(:,:),res)
892	miniwsqs = 1.E27
893	maxiwsqs =-1.E27
894	plot_wq = gsn_csm_xy(wks,t,data_0(:,:),res)
895	miniwq = 1.E27
896	maxiwq =-1.E27
897	plot_wpqvp = gsn_csm_xy(wks,t,data_0(:,:),res)
898	miniwpqvp = 1.E27
899	maxiwpqvp =-1.E27
900	plot_wsqvs = gsn_csm_xy(wks,t,data_0(:,:),res)
901	miniwsqvs = 1.E27
902	maxiwsqvs =-1.E27
903	plot_wqv = gsn_csm_xy(wks,t,data_0(:,:),res)
904	miniwqv = 1.E27
905	maxiwqv =-1.E27
906	plot_wpsp = gsn_csm_xy(wks,t,data_0(:,:),res)
907	miniwpsp = 1.E27
908	maxiwpsp =-1.E27
909	plot_wsss = gsn_csm_xy(wks,t,data_0(:,:),res)
910	miniwsss = 1.E27
911	maxiwsss =-1.E27
912	plot_ws = gsn_csm_xy(wks,t,data_0(:,:),res)
913	miniws = 1.E27
914	maxiws =-1.E27
915	plot_wpsap = gsn_csm_xy(wks,t,data_0(:,:),res)
916	miniwpsap = 1.E27
917	maxiwpsap =-1.E27
918	plot_wssas = gsn_csm_xy(wks,t,data_0(:,:),res)
919	miniwssas = 1.E27
920	maxiwssas =-1.E27
921	plot_wsa = gsn_csm_xy(wks,t,data_0(:,:),res)
922	miniwsa = 1.E27
923	maxiwsa =-1.E27
924	plot_wses = gsn_csm_xy(wks,t,data_0(:,:),res)
925	plot_us2 = gsn_csm_xy(wks,t,data_0(:,:),res)
926	minius2 = 1.E27
927	maxius2 =-1.E27
928	plot_vs2 = gsn_csm_xy(wks,t,data_0(:,:),res)
929	minivs2 = 1.E27
930	maxivs2 =-1.E27
931	plot_ws2 = gsn_csm_xy(wks,t,data_0(:,:),res)
932	miniws2 = 1.E27
933	maxiws2 =-1.E27
934	plot_pts2 = gsn_csm_xy(wks,t,data_0(:,:),res)
935	plot_ws3 = gsn_csm_xy(wks,t,data_0(:,:),res)
936	plot_Sw = gsn_csm_xy(wks,t,data_0(:,:),res)
937	plot_ws2pts = gsn_csm_xy(wks,t,data_0(:,:),res)
938	plot_wspts2 = gsn_csm_xy(wks,t,data_0(:,:),res)
939	plot_wsususodz = gsn_csm_xy(wks,t,data_0(:,:),res)
940	miniwsususodz = 1.E27
941	maxiwsususodz =-1.E27
942	plot_wspsodz = gsn_csm_xy(wks,t,data_0(:,:),res)
943	miniwspsodz = 1.E27
944	maxiwspsodz =-1.E27
945	plot_wpeodz = gsn_csm_xy(wks,t,data_0(:,:),res)
946	miniwpeodz = 1.E27
947	maxiwpeodz =-1.E27
948	end if
949
950	n_o=0
951	count_var=0
952
953	do varn = 0,dim-1
954
955	check = True
956
957	if (prof3d .EQ. 0) then
958	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
959	check = False
960	end if
961	else
962	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "zusi") .OR. isStrSubset( vNam(varn), "zwwi") .OR. isStrSubset( vNam(varn), "x") .OR. isStrSubset( vNam(varn), "xu") .OR. isStrSubset( vNam(varn), "y") .OR. isStrSubset( vNam(varn), "yv") .OR. isStrSubset( vNam(varn), "zu_3d") .OR. isStrSubset( vNam(varn), "zw_3d")) then
963	check = False
964	end if
965	end if
966
967	if (.not. isvar("var")) then
968	if (parameter(21) .NE. "variables") then
969	var=parameter(21)
970	check = isStrSubset( var,","+vNam(varn)+"," )
971	end if
972	else
973	check = isStrSubset( var,","+vNam(varn)+"," )
974	end if
975
976	if (combine .EQ. 1) then
977	com=isStrSubset(c_var,","+vNam(varn)+"," )
978	if (com)
979	if (prof3d .EQ. 0) then
980	temp = f->$vNam(varn)$
981	data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1)
982	else
983	do i=0,dimz-1
984	do j=start_time_step,end_time_step
985	temp= f->$vNam(varn)$(j,i,:,:)
986	data(varn,j-start_time_step,i) = dim_avg_Wrap(dim_avg_Wrap(temp(:,:)))
987	end do
988	end do
989	print(" ")
990	print("Variable for combine '"+vNam(varn)+"' is read")
991	print(" ")
992	end if
993	unit(varn) = temp@units
994	if (n_o .GT. number_comb-1) then
995	print(" ")
996	print("Please set 'number_comb' to the right number of overlaying variables ('c_var')")
997	print(" ")
998	exit
999	end if
1000	mini(n_o)=min(data(varn,:,:))
1001	maxi(n_o)=max(data(varn,:,:))
1002	n_o=n_o+1
1003	end if
1004	end if
1005
1006	if(check) then
1007
1008	count_var=count_var+1
1009
1010	if (prof3d .EQ. 0) then
1011	temp = f->$vNam(varn)$
1012	else
1013	do i=0,dimz-1
1014	do j=start_time_step,end_time_step
1015	temp = f->$vNam(varn)$(j,i,:,:)
1016	data(varn,j-start_time_step,i) = dim_avg_Wrap(dim_avg_Wrap(temp(:,:)))
1017	end do
1018	end do
1019	print(" ")
1020	print("Variable '"+vNam(varn)+"' is read")
1021	print(" ")
1022	unit(varn) = temp@units
1023	a=getvaratts(temp)
1024	b=dimsizes(a)
1025	end if
1026
1027	if (prof3d .EQ. 0) then
1028	z = f->$vNam(varn+1)$
1029	unit(varn) = temp@units
1030	data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1)
1031	else
1032	do i=0,b-1
1033	if (isStrSubset( a(i),"zu_3d" ))
1034	z_v(varn,:) = z_u
1035	z = z_v(varn,:)
1036	else
1037	if (isStrSubset( a(i),"zw_3d" ))
1038	z_v(varn,:) = z_w
1039	z = z_v(varn,:)
1040	end if
1041	end if
1042	end do
1043	end if
1044	if (over .EQ. 0) then
1045	res@gsnLeftString = vNam(varn)
1046	res@gsnRightString = unit(varn)
1047	res@trYMinF = min_z
1048	res@trYMaxF = max_z
1049	if (.not. isvar("xs")) then
1050	if (parameter(63) .NE. "x0") then
1051	res@trXMinF = stringtofloat(parameter(63))
1052	else
1053	res@trXMinF = min(data(varn,:,:))
1054	end if
1055	else
1056	res@trXMinF = xs
1057	end if
1058	if (.not. isvar("xe")) then
1059	if (parameter(65) .NE. "xdim") then
1060	res@trXMaxF = stringtofloat(parameter(65))
1061	else
1062	res@trXMaxF = max(data(varn,:,:))
1063	end if
1064	else
1065	res@trXMaxF = xe
1066	end if
1067	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res)
1068	end if
1069
1070	if (vNam(varn) .EQ. "u") then
1071	miniu=min(data(varn,:,:))
1072	maxiu=max(data(varn,:,:))
1073	if (over .EQ. 1) then
1074	res@xyDashPattern = 0
1075	plot_u = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1076	else
1077	res@gsnLeftString = vNam(varn)
1078	res@gsnRightString = unit(varn)
1079	if (.not. isvar("xs")) then
1080	if (parameter(63) .NE. "x0") then
1081	res@trXMinF = stringtofloat(parameter(63))
1082	else
1083	res@trXMinF = miniu
1084	end if
1085	else
1086	res@trXMinF = xs
1087	end if
1088	if (.not. isvar("xe")) then
1089	if (parameter(65) .NE. "xdim") then
1090	res@trXMaxF = stringtofloat(parameter(65))
1091	else
1092	res@trXMaxF = maxiu
1093	end if
1094	else
1095	res@trXMaxF = xe
1096	end if
1097	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1098	end if
1099	end if
1100	if (vNam(varn) .EQ. "v") then
1101	miniv=min(data(varn,:,:))
1102	maxiv=max(data(varn,:,:))
1103	if (over .EQ. 1) then
1104	res@xyDashPattern = 1
1105	plot_v = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1106	else
1107	res@gsnLeftString = vNam(varn)
1108	res@gsnRightString = unit(varn)
1109	if (.not. isvar("xs")) then
1110	if (parameter(63) .NE. "x0") then
1111	res@trXMinF = stringtofloat(parameter(63))
1112	else
1113	res@trXMinF = miniv
1114	end if
1115	else
1116	res@trXMinF = xs
1117	end if
1118	if (.not. isvar("xe")) then
1119	if (parameter(65) .NE. "xdim") then
1120	res@trXMaxF = stringtofloat(parameter(65))
1121	else
1122	res@trXMaxF = maxiv
1123	end if
1124	else
1125	res@trXMaxF = xe
1126	end if
1127	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1128	end if
1129	end if
1130	if (vNam(varn) .EQ. "w") then
1131	miniw=min(data(varn,:,:))
1132	maxiw=max(data(varn,:,:))
1133	if (over .EQ. 1) then
1134	res@xyDashPattern = 2
1135	plot_w = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1136	else
1137	res@gsnLeftString = vNam(varn)
1138	res@gsnRightString = unit(varn)
1139	if (.not. isvar("xs")) then
1140	if (parameter(63) .NE. "x0") then
1141	res@trXMinF = stringtofloat(parameter(63))
1142	else
1143	res@trXMinF = miniw
1144	end if
1145	else
1146	res@trXMinF = xs
1147	end if
1148	if (.not. isvar("xe")) then
1149	if (parameter(65) .NE. "xdim") then
1150	res@trXMaxF = stringtofloat(parameter(65))
1151	else
1152	res@trXMaxF = maxiw
1153	end if
1154	else
1155	res@trXMaxF = xe
1156	end if
1157	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1158	end if
1159	end if
1160
1161	if (vNam(varn) .EQ. "pt") then
1162	minipt=min(data(varn,:,:))
1163	maxipt=max(data(varn,:,:))
1164	if (over .EQ. 1) then
1165	res@xyDashPattern = 0
1166	plot_pt = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1167	else
1168	res@gsnLeftString = vNam(varn)
1169	res@gsnRightString = unit(varn)
1170	if (.not. isvar("xs")) then
1171	if (parameter(63) .NE. "x0") then
1172	res@trXMinF = stringtofloat(parameter(63))
1173	else
1174	res@trXMinF = minipt
1175	end if
1176	else
1177	res@trXMinF = xs
1178	end if
1179	if (.not. isvar("xe")) then
1180	if (parameter(65) .NE. "xdim") then
1181	res@trXMaxF = stringtofloat(parameter(65))
1182	else
1183	res@trXMaxF = maxipt
1184	end if
1185	else
1186	res@trXMaxF = xe
1187	end if
1188	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1189	end if
1190	end if
1191	if (vNam(varn) .EQ. "vpt") then
1192	minivpt=min(data(varn,:,:))
1193	maxivpt=max(data(varn,:,:))
1194	if (over .EQ. 1) then
1195	res@xyDashPattern = 1
1196	plot_vpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1197	else
1198	res@gsnLeftString = vNam(varn)
1199	res@gsnRightString = unit(varn)
1200	if (.not. isvar("xs")) then
1201	if (parameter(63) .NE. "x0") then
1202	res@trXMinF = stringtofloat(parameter(63))
1203	else
1204	res@trXMinF = minivpt
1205	end if
1206	else
1207	res@trXMinF = xs
1208	end if
1209	if (.not. isvar("xe")) then
1210	if (parameter(65) .NE. "xdim") then
1211	res@trXMaxF = stringtofloat(parameter(65))
1212	else
1213	res@trXMaxF = maxivpt
1214	end if
1215	else
1216	res@trXMaxF = xe
1217	end if
1218	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1219	end if
1220	end if
1221	if (vNam(varn) .EQ. "lpt") then
1222	minilpt=min(data(varn,:,:))
1223	maxilpt=max(data(varn,:,:))
1224	if (over .EQ. 1) then
1225	res@xyDashPattern = 2
1226	plot_lpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1227	else
1228	res@gsnLeftString = vNam(varn)
1229	res@gsnRightString = unit(varn)
1230	if (.not. isvar("xs")) then
1231	if (parameter(63) .NE. "x0") then
1232	res@trXMinF = stringtofloat(parameter(63))
1233	else
1234	res@trXMinF = minilpt
1235	end if
1236	else
1237	res@trXMinF = xs
1238	end if
1239	if (.not. isvar("xe")) then
1240	if (parameter(65) .NE. "xdim") then
1241	res@trXMaxF = stringtofloat(parameter(65))
1242	else
1243	res@trXMaxF = maxilpt
1244	end if
1245	else
1246	res@trXMaxF = xe
1247	end if
1248	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1249	end if
1250	end if
1251
1252	if (vNam(varn) .EQ. "q") then
1253	miniq=min(data(varn,:,:))
1254	maxiq=max(data(varn,:,:))
1255	if (over .EQ. 1) then
1256	res@xyDashPattern = 0
1257	plot_q = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1258	else
1259	res@gsnLeftString = vNam(varn)
1260	res@gsnRightString = unit(varn)
1261	if (.not. isvar("xs")) then
1262	if (parameter(63) .NE. "x0") then
1263	res@trXMinF = stringtofloat(parameter(63))
1264	else
1265	res@trXMinF = minilq
1266	end if
1267	else
1268	res@trXMinF = xs
1269	end if
1270	if (.not. isvar("xe")) then
1271	if (parameter(65) .NE. "xdim") then
1272	res@trXMaxF = stringtofloat(parameter(65))
1273	else
1274	res@trXMaxF = maxilq
1275	end if
1276	else
1277	res@trXMaxF = xe
1278	end if
1279	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1280	end if
1281	end if
1282	if (vNam(varn) .EQ. "qv") then
1283	miniqv=min(data(varn,:,:))
1284	maxiqv=max(data(varn,:,:))
1285	if (over .EQ. 1) then
1286	res@xyDashPattern = 1
1287	plot_qv = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1288	else
1289	res@gsnLeftString = vNam(varn)
1290	res@gsnRightString = unit(varn)
1291	if (.not. isvar("xs")) then
1292	if (parameter(63) .NE. "x0") then
1293	res@trXMinF = stringtofloat(parameter(63))
1294	else
1295	res@trXMinF = minilqv
1296	end if
1297	else
1298	res@trXMinF = xs
1299	end if
1300	if (.not. isvar("xe")) then
1301	if (parameter(65) .NE. "xdim") then
1302	res@trXMaxF = stringtofloat(parameter(65))
1303	else
1304	res@trXMaxF = maxilqv
1305	end if
1306	else
1307	res@trXMaxF = xe
1308	end if
1309	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1310	end if
1311	end if
1312	if (vNam(varn) .EQ. "ql") then
1313	miniql=min(data(varn,:,:))
1314	maxiql=max(data(varn,:,:))
1315	if (over .EQ. 1) then
1316	res@xyDashPattern = 2
1317	plot_ql = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1318	else
1319	res@gsnLeftString = vNam(varn)
1320	res@gsnRightString = unit(varn)
1321	if (.not. isvar("xs")) then
1322	if (parameter(63) .NE. "x0") then
1323	res@trXMinF = stringtofloat(parameter(63))
1324	else
1325	res@trXMinF = miniql
1326	end if
1327	else
1328	res@trXMinF = xs
1329	end if
1330	if (.not. isvar("xe")) then
1331	if (parameter(65) .NE. "xdim") then
1332	res@trXMaxF = stringtofloat(parameter(65))
1333	else
1334	res@trXMaxF = maxiql
1335	end if
1336	else
1337	res@trXMaxF = xe
1338	end if
1339	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1340	end if
1341	end if
1342
1343	if (vNam(varn) .EQ. "e") then
1344	minie=min(data(varn,:,:))
1345	maxie=max(data(varn,:,:))
1346	if (over .EQ. 1) then
1347	res@xyDashPattern = 0
1348	plot_e = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1349	else
1350	res@gsnLeftString = vNam(varn)
1351	res@gsnRightString = unit(varn)
1352	if (.not. isvar("xs")) then
1353	if (parameter(63) .NE. "x0") then
1354	res@trXMinF = stringtofloat(parameter(63))
1355	else
1356	res@trXMinF = minie
1357	end if
1358	else
1359	res@trXMinF = xs
1360	end if
1361	if (.not. isvar("xe")) then
1362	if (parameter(65) .NE. "xdim") then
1363	res@trXMaxF = stringtofloat(parameter(65))
1364	else
1365	res@trXMaxF = maxie
1366	end if
1367	else
1368	res@trXMaxF = xe
1369	end if
1370	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1371	end if
1372	end if
1373	if (vNam(varn) .EQ. "es") then
1374	minies=min(data(varn,:,:))
1375	maxies=max(data(varn,:,:))
1376	if (over .EQ. 1) then
1377	res@xyDashPattern = 1
1378	plot_es = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1379	else
1380	res@gsnLeftString = vNam(varn)
1381	res@gsnRightString = unit(varn)
1382	if (.not. isvar("xs")) then
1383	if (parameter(63) .NE. "x0") then
1384	res@trXMinF = stringtofloat(parameter(63))
1385	else
1386	res@trXMinF = minies
1387	end if
1388	else
1389	res@trXMinF = xs
1390	end if
1391	if (.not. isvar("xe")) then
1392	if (parameter(65) .NE. "xdim") then
1393	res@trXMaxF = stringtofloat(parameter(65))
1394	else
1395	res@trXMaxF = maxies
1396	end if
1397	else
1398	res@trXMaxF = xe
1399	end if
1400	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1401	end if
1402	end if
1403
1404	if (vNam(varn) .EQ. "km") then
1405	minikm=min(data(varn,:,:))
1406	maxikm=max(data(varn,:,:))
1407	if (over .EQ. 1) then
1408	res@xyDashPattern = 0
1409	plot_km = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1410	else
1411	res@gsnLeftString = vNam(varn)
1412	res@gsnRightString = unit(varn)
1413	if (.not. isvar("xs")) then
1414	if (parameter(63) .NE. "x0") then
1415	res@trXMinF = stringtofloat(parameter(63))
1416	else
1417	res@trXMinF = minikm
1418	end if
1419	else
1420	res@trXMinF = xs
1421	end if
1422	if (.not. isvar("xe")) then
1423	if (parameter(65) .NE. "xdim") then
1424	res@trXMaxF = stringtofloat(parameter(65))
1425	else
1426	res@trXMaxF = maxikm
1427	end if
1428	else
1429	res@trXMaxF = xe
1430	end if
1431	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1432	end if
1433	end if
1434	if (vNam(varn) .EQ. "kh") then
1435	minikh=min(data(varn,:,:))
1436	maxikh=max(data(varn,:,:))
1437	if (over .EQ. 1) then
1438	res@xyDashPattern = 1
1439	plot_kh = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1440	else
1441	res@gsnLeftString = vNam(varn)
1442	res@gsnRightString = unit(varn)
1443	if (.not. isvar("xs")) then
1444	if (parameter(63) .NE. "x0") then
1445	res@trXMinF = stringtofloat(parameter(63))
1446	else
1447	res@trXMinF = minikh
1448	end if
1449	else
1450	res@trXMinF = xs
1451	end if
1452	if (.not. isvar("xe")) then
1453	if (parameter(65) .NE. "xdim") then
1454	res@trXMaxF = stringtofloat(parameter(65))
1455	else
1456	res@trXMaxF = maxikh
1457	end if
1458	else
1459	res@trXMaxF = xe
1460	end if
1461	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1462	end if
1463	end if
1464
1465	if (vNam(varn) .EQ. "wpup") then
1466	miniwpup=min(data(varn,:,:))
1467	maxiwpup=max(data(varn,:,:))
1468	if (over .EQ. 1) then
1469	res@xyDashPattern = 0
1470	plot_wpup = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1471	else
1472	res@gsnLeftString = vNam(varn)
1473	res@gsnRightString = unit(varn)
1474	if (.not. isvar("xs")) then
1475	if (parameter(63) .NE. "x0") then
1476	res@trXMinF = stringtofloat(parameter(63))
1477	else
1478	res@trXMinF = miniwpup
1479	end if
1480	else
1481	res@trXMinF = xs
1482	end if
1483	if (.not. isvar("xe")) then
1484	if (parameter(65) .NE. "xdim") then
1485	res@trXMaxF = stringtofloat(parameter(65))
1486	else
1487	res@trXMaxF = maxiwpup
1488	end if
1489	else
1490	res@trXMaxF = xe
1491	end if
1492	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1493	end if
1494	end if
1495	if (vNam(varn) .EQ. "wsus") then
1496	miniwsus=min(data(varn,:,:))
1497	maxiwsus=max(data(varn,:,:))
1498	if (over .EQ. 1) then
1499	res@xyDashPattern = 1
1500	plot_wsus = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1501	else
1502	res@gsnLeftString = vNam(varn)
1503	res@gsnRightString = unit(varn)
1504	if (.not. isvar("xs")) then
1505	if (parameter(63) .NE. "x0") then
1506	res@trXMinF = stringtofloat(parameter(63))
1507	else
1508	res@trXMinF = miniwsus
1509	end if
1510	else
1511	res@trXMinF = xs
1512	end if
1513	if (.not. isvar("xe")) then
1514	if (parameter(65) .NE. "xdim") then
1515	res@trXMaxF = stringtofloat(parameter(65))
1516	else
1517	res@trXMaxF = maxiwsus
1518	end if
1519	else
1520	res@trXMaxF = xe
1521	end if
1522	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1523	end if
1524	end if
1525	if (vNam(varn) .EQ. "wu") then
1526	miniwu=min(data(varn,:,:))
1527	maxiwu=max(data(varn,:,:))
1528	if (over .EQ. 1) then
1529	res@xyDashPattern = 2
1530	plot_wu = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1531	else
1532	res@gsnLeftString = vNam(varn)
1533	res@gsnRightString = unit(varn)
1534	if (.not. isvar("xs")) then
1535	if (parameter(63) .NE. "x0") then
1536	res@trXMinF = stringtofloat(parameter(63))
1537	else
1538	res@trXMinF = miniwu
1539	end if
1540	else
1541	res@trXMinF = xs
1542	end if
1543	if (.not. isvar("xe")) then
1544	if (parameter(65) .NE. "xdim") then
1545	res@trXMaxF = stringtofloat(parameter(65))
1546	else
1547	res@trXMaxF = maxiwu
1548	end if
1549	else
1550	res@trXMaxF = xe
1551	end if
1552	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1553	end if
1554	end if
1555
1556	if (vNam(varn) .EQ. "wpvp") then
1557	miniwpvp=min(data(varn,:,:))
1558	maxiwpvp=max(data(varn,:,:))
1559	if (over .EQ. 1) then
1560	res@xyDashPattern = 0
1561	plot_wpvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1562	else
1563	res@gsnLeftString = vNam(varn)
1564	res@gsnRightString = unit(varn)
1565	if (.not. isvar("xs")) then
1566	if (parameter(63) .NE. "x0") then
1567	res@trXMinF = stringtofloat(parameter(63))
1568	else
1569	res@trXMinF = miniwpvp
1570	end if
1571	else
1572	res@trXMinF = xs
1573	end if
1574	if (.not. isvar("xe")) then
1575	if (parameter(65) .NE. "xdim") then
1576	res@trXMaxF = stringtofloat(parameter(65))
1577	else
1578	res@trXMaxF = maxiwpvp
1579	end if
1580	else
1581	res@trXMaxF = xe
1582	end if
1583	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1584	end if
1585	end if
1586	if (vNam(varn) .EQ. "wsvs") then
1587	miniwsvs=min(data(varn,:,:))
1588	maxiwsvs=max(data(varn,:,:))
1589	if (over .EQ. 1) then
1590	res@xyDashPattern = 1
1591	plot_wsvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1592	else
1593	res@gsnLeftString = vNam(varn)
1594	res@gsnRightString = unit(varn)
1595	if (.not. isvar("xs")) then
1596	if (parameter(63) .NE. "x0") then
1597	res@trXMinF = stringtofloat(parameter(63))
1598	else
1599	res@trXMinF = miniwsvs
1600	end if
1601	else
1602	res@trXMinF = xs
1603	end if
1604	if (.not. isvar("xe")) then
1605	if (parameter(65) .NE. "xdim") then
1606	res@trXMaxF = stringtofloat(parameter(65))
1607	else
1608	res@trXMaxF = maxiwsvs
1609	end if
1610	else
1611	res@trXMaxF = xe
1612	end if
1613	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1614	end if
1615	end if
1616	if (vNam(varn) .EQ. "wv") then
1617	miniwv=min(data(varn,:,:))
1618	maxiwv=max(data(varn,:,:))
1619	if (over .EQ. 1) then
1620	res@xyDashPattern = 2
1621	plot_wv = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1622	else
1623	res@gsnLeftString = vNam(varn)
1624	res@gsnRightString = unit(varn)
1625	if (.not. isvar("xs")) then
1626	if (parameter(63) .NE. "x0") then
1627	res@trXMinF = stringtofloat(parameter(63))
1628	else
1629	res@trXMinF = miniwv
1630	end if
1631	else
1632	res@trXMinF = xs
1633	end if
1634	if (.not. isvar("xe")) then
1635	if (parameter(65) .NE. "xdim") then
1636	res@trXMaxF = stringtofloat(parameter(65))
1637	else
1638	res@trXMaxF = maxiwv
1639	end if
1640	else
1641	res@trXMaxF = xe
1642	end if
1643	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1644	end if
1645	end if
1646
1647	if (vNam(varn) .EQ. "wpptp") then
1648	miniwpptp=min(data(varn,:,:))
1649	maxiwpptp=max(data(varn,:,:))
1650	if (over .EQ. 1) then
1651	res@xyDashPattern = 0
1652	plot_wpptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1653	else
1654	res@gsnLeftString = vNam(varn)
1655	res@gsnRightString = unit(varn)
1656	if (.not. isvar("xs")) then
1657	if (parameter(63) .NE. "x0") then
1658	res@trXMinF = stringtofloat(parameter(63))
1659	else
1660	res@trXMinF = miniwpptp
1661	end if
1662	else
1663	res@trXMinF = xs
1664	end if
1665	if (.not. isvar("xe")) then
1666	if (parameter(65) .NE. "xdim") then
1667	res@trXMaxF = stringtofloat(parameter(65))
1668	else
1669	res@trXMaxF = maxiwpptp
1670	end if
1671	else
1672	res@trXMaxF = xe
1673	end if
1674	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1675	end if
1676	end if
1677	if (vNam(varn) .EQ. "wspts") then
1678	miniwspts=min(data(varn,:,:))
1679	maxiwspts=max(data(varn,:,:))
1680	if (over .EQ. 1) then
1681	res@xyDashPattern = 1
1682	plot_wspts = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1683	else
1684	res@gsnLeftString = vNam(varn)
1685	res@gsnRightString = unit(varn)
1686	if (.not. isvar("xs")) then
1687	if (parameter(63) .NE. "x0") then
1688	res@trXMinF = stringtofloat(parameter(63))
1689	else
1690	res@trXMinF = miniwspts
1691	end if
1692	else
1693	res@trXMinF = xs
1694	end if
1695	if (.not. isvar("xe")) then
1696	if (parameter(65) .NE. "xdim") then
1697	res@trXMaxF = stringtofloat(parameter(65))
1698	else
1699	res@trXMaxF = maxiwspts
1700	end if
1701	else
1702	res@trXMaxF = xe
1703	end if
1704	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1705	end if
1706	end if
1707	if (vNam(varn) .EQ. "wpt") then
1708	miniwpt=min(data(varn,:,:))
1709	maxiwpt=max(data(varn,:,:))
1710	if (over .EQ. 1) then
1711	res@xyDashPattern = 2
1712	plot_wpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1713	else
1714	res@gsnLeftString = vNam(varn)
1715	res@gsnRightString = unit(varn)
1716	if (.not. isvar("xs")) then
1717	if (parameter(63) .NE. "x0") then
1718	res@trXMinF = stringtofloat(parameter(63))
1719	else
1720	res@trXMinF = miniwpt
1721	end if
1722	else
1723	res@trXMinF = xs
1724	end if
1725	if (.not. isvar("xe")) then
1726	if (parameter(65) .NE. "xdim") then
1727	res@trXMaxF = stringtofloat(parameter(65))
1728	else
1729	res@trXMaxF = maxiwpt
1730	end if
1731	else
1732	res@trXMaxF = xe
1733	end if
1734	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1735	end if
1736	end if
1737
1738	if (vNam(varn) .EQ. "wsptsBC") then
1739	miniwsptsBC=min(data(varn,:,:))
1740	maxiwsptsBC=max(data(varn,:,:))
1741	if (over .EQ. 1) then
1742	res@xyDashPattern = 0
1743	plot_wsptsBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1744	else
1745	res@gsnLeftString = vNam(varn)
1746	res@gsnRightString = unit(varn)
1747	if (.not. isvar("xs")) then
1748	if (parameter(63) .NE. "x0") then
1749	res@trXMinF = stringtofloat(parameter(63))
1750	else
1751	res@trXMinF = miniwsptsBC
1752	end if
1753	else
1754	res@trXMinF = xs
1755	end if
1756	if (.not. isvar("xe")) then
1757	if (parameter(65) .NE. "xdim") then
1758	res@trXMaxF = stringtofloat(parameter(65))
1759	else
1760	res@trXMaxF = maxiwsptsBC
1761	end if
1762	else
1763	res@trXMaxF = xe
1764	end if
1765	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1766	end if
1767	end if
1768	if (vNam(varn) .EQ. "wptBC") then
1769	miniwptBC=min(data(varn,:,:))
1770	maxiwptBC=max(data(varn,:,:))
1771	if (over .EQ. 1) then
1772	res@xyDashPattern = 1
1773	plot_wptBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1774	else
1775	res@gsnLeftString = vNam(varn)
1776	res@gsnRightString = unit(varn)
1777	if (.not. isvar("xs")) then
1778	if (parameter(63) .NE. "x0") then
1779	res@trXMinF = stringtofloat(parameter(63))
1780	else
1781	res@trXMinF = miniwptBC
1782	end if
1783	else
1784	res@trXMinF = xs
1785	end if
1786	if (.not. isvar("xe")) then
1787	if (parameter(65) .NE. "xdim") then
1788	res@trXMaxF = stringtofloat(parameter(65))
1789	else
1790	res@trXMaxF = maxiwptBC
1791	end if
1792	else
1793	res@trXMaxF = xe
1794	end if
1795	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1796	end if
1797	end if
1798
1799	if (vNam(varn) .EQ. "wpvptp") then
1800	miniwpvptp=min(data(varn,:,:))
1801	maxiwpvptp=max(data(varn,:,:))
1802	if (over .EQ. 1) then
1803	res@xyDashPattern = 0
1804	plot_wpvptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1805	else
1806	res@gsnLeftString = vNam(varn)
1807	res@gsnRightString = unit(varn)
1808	if (.not. isvar("xs")) then
1809	if (parameter(63) .NE. "x0") then
1810	res@trXMinF = stringtofloat(parameter(63))
1811	else
1812	res@trXMinF = miniwpvptp
1813	end if
1814	else
1815	res@trXMinF = xs
1816	end if
1817	if (.not. isvar("xe")) then
1818	if (parameter(65) .NE. "xdim") then
1819	res@trXMaxF = stringtofloat(parameter(65))
1820	else
1821	res@trXMaxF = maxiwpvptp
1822	end if
1823	else
1824	res@trXMaxF = xe
1825	end if
1826	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1827	end if
1828	end if
1829	if (vNam(varn) .EQ. "wsvpts") then
1830	miniwsvpts=min(data(varn,:,:))
1831	maxiwsvpts=max(data(varn,:,:))
1832	if (over .EQ. 1) then
1833	res@xyDashPattern = 1
1834	plot_wsvpts = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1835	else
1836	res@gsnLeftString = vNam(varn)
1837	res@gsnRightString = unit(varn)
1838	if (.not. isvar("xs")) then
1839	if (parameter(63) .NE. "x0") then
1840	res@trXMinF = stringtofloat(parameter(63))
1841	else
1842	res@trXMinF = miniwsvpts
1843	end if
1844	else
1845	res@trXMinF = xs
1846	end if
1847	if (.not. isvar("xe")) then
1848	if (parameter(65) .NE. "xdim") then
1849	res@trXMaxF = stringtofloat(parameter(65))
1850	else
1851	res@trXMaxF = maxiwsvpts
1852	end if
1853	else
1854	res@trXMaxF = xe
1855	end if
1856	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1857	end if
1858	end if
1859	if (vNam(varn) .EQ. "wvpt") then
1860	miniwvpt=min(data(varn,:,:))
1861	maxiwvpt=max(data(varn,:,:))
1862	if (over .EQ. 1) then
1863	res@xyDashPattern = 2
1864	plot_wvpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1865	else
1866	res@gsnLeftString = vNam(varn)
1867	res@gsnRightString = unit(varn)
1868	if (.not. isvar("xs")) then
1869	if (parameter(63) .NE. "x0") then
1870	res@trXMinF = stringtofloat(parameter(63))
1871	else
1872	res@trXMinF = miniwvpt
1873	end if
1874	else
1875	res@trXMinF = xs
1876	end if
1877	if (.not. isvar("xe")) then
1878	if (parameter(65) .NE. "xdim") then
1879	res@trXMaxF = stringtofloat(parameter(65))
1880	else
1881	res@trXMaxF = maxiwvpt
1882	end if
1883	else
1884	res@trXMaxF = xe
1885	end if
1886	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1887	end if
1888	end if
1889
1890	if (vNam(varn) .EQ. "wpqp") then
1891	miniwpqp=min(data(varn,:,:))
1892	maxiwpqp=max(data(varn,:,:))
1893	if (over .EQ. 1) then
1894	res@xyDashPattern = 0
1895	plot_wpqp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1896	else
1897	res@gsnLeftString = vNam(varn)
1898	res@gsnRightString = unit(varn)
1899	if (.not. isvar("xs")) then
1900	if (parameter(63) .NE. "x0") then
1901	res@trXMinF = stringtofloat(parameter(63))
1902	else
1903	res@trXMinF = miniwpqp
1904	end if
1905	else
1906	res@trXMinF = xs
1907	end if
1908	if (.not. isvar("xe")) then
1909	if (parameter(65) .NE. "xdim") then
1910	res@trXMaxF = stringtofloat(parameter(65))
1911	else
1912	res@trXMaxF = maxiwpqp
1913	end if
1914	else
1915	res@trXMaxF = xe
1916	end if
1917	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1918	end if
1919	end if
1920	if (vNam(varn) .EQ. "wsqs") then
1921	miniwsqs=min(data(varn,:,:))
1922	maxiwsqs=max(data(varn,:,:))
1923	if (over .EQ. 1) then
1924	res@xyDashPattern = 1
1925	plot_wsqs = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1926	else
1927	res@gsnLeftString = vNam(varn)
1928	res@gsnRightString = unit(varn)
1929	if (.not. isvar("xs")) then
1930	if (parameter(63) .NE. "x0") then
1931	res@trXMinF = stringtofloat(parameter(63))
1932	else
1933	res@trXMinF = miniwsqs
1934	end if
1935	else
1936	res@trXMinF = xs
1937	end if
1938	if (.not. isvar("xe")) then
1939	if (parameter(65) .NE. "xdim") then
1940	res@trXMaxF = stringtofloat(parameter(65))
1941	else
1942	res@trXMaxF = maxiwsqs
1943	end if
1944	else
1945	res@trXMaxF = xe
1946	end if
1947	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1948	end if
1949	end if
1950	if (vNam(varn) .EQ. "wq") then
1951	miniwq=min(data(varn,:,:))
1952	maxiwq=max(data(varn,:,:))
1953	if (over .EQ. 1) then
1954	res@xyDashPattern = 2
1955	plot_wq = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1956	else
1957	res@gsnLeftString = vNam(varn)
1958	res@gsnRightString = unit(varn)
1959	if (.not. isvar("xs")) then
1960	if (parameter(63) .NE. "x0") then
1961	res@trXMinF = stringtofloat(parameter(63))
1962	else
1963	res@trXMinF = miniwq
1964	end if
1965	else
1966	res@trXMinF = xs
1967	end if
1968	if (.not. isvar("xe")) then
1969	if (parameter(65) .NE. "xdim") then
1970	res@trXMaxF = stringtofloat(parameter(65))
1971	else
1972	res@trXMaxF = maxiwq
1973	end if
1974	else
1975	res@trXMaxF = xe
1976	end if
1977	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1978	end if
1979	end if
1980
1981	if (vNam(varn) .EQ. "wpqvp") then
1982	miniwpqvp=min(data(varn,:,:))
1983	maxiwpqvp=max(data(varn,:,:))
1984	if (over .EQ. 1) then
1985	res@xyDashPattern = 0
1986	plot_wpqvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1987	else
1988	res@gsnLeftString = vNam(varn)
1989	res@gsnRightString = unit(varn)
1990	if (.not. isvar("xs")) then
1991	if (parameter(63) .NE. "x0") then
1992	res@trXMinF = stringtofloat(parameter(63))
1993	else
1994	res@trXMinF = miniwpqvp
1995	end if
1996	else
1997	res@trXMinF = xs
1998	end if
1999	if (.not. isvar("xe")) then
2000	if (parameter(65) .NE. "xdim") then
2001	res@trXMaxF = stringtofloat(parameter(65))
2002	else
2003	res@trXMaxF = maxiwpqvp
2004	end if
2005	else
2006	res@trXMaxF = xe
2007	end if
2008	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2009	end if
2010	end if
2011	if (vNam(varn) .EQ. "wsqvs") then
2012	miniwsqvs=min(data(varn,:,:))
2013	maxiwsqvs=max(data(varn,:,:))
2014	if (over .EQ. 1) then
2015	res@xyDashPattern = 1
2016	plot_wsqvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2017	else
2018	res@gsnLeftString = vNam(varn)
2019	res@gsnRightString = unit(varn)
2020	if (.not. isvar("xs")) then
2021	if (parameter(63) .NE. "x0") then
2022	res@trXMinF = stringtofloat(parameter(63))
2023	else
2024	res@trXMinF = miniwsqvs
2025	end if
2026	else
2027	res@trXMinF = xs
2028	end if
2029	if (.not. isvar("xe")) then
2030	if (parameter(65) .NE. "xdim") then
2031	res@trXMaxF = stringtofloat(parameter(65))
2032	else
2033	res@trXMaxF = maxiwsqvs
2034	end if
2035	else
2036	res@trXMaxF = xe
2037	end if
2038	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2039	end if
2040	end if
2041	if (vNam(varn) .EQ. "wqv") then
2042	miniwqv=min(data(varn,:,:))
2043	maxiwqv=max(data(varn,:,:))
2044	if (over .EQ. 1) then
2045	res@xyDashPattern = 2
2046	plot_wqv = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2047	else
2048	res@gsnLeftString = vNam(varn)
2049	res@gsnRightString = unit(varn)
2050	if (.not. isvar("xs")) then
2051	if (parameter(63) .NE. "x0") then
2052	res@trXMinF = stringtofloat(parameter(63))
2053	else
2054	res@trXMinF = miniwqv
2055	end if
2056	else
2057	res@trXMinF = xs
2058	end if
2059	if (.not. isvar("xe")) then
2060	if (parameter(65) .NE. "xdim") then
2061	res@trXMaxF = stringtofloat(parameter(65))
2062	else
2063	res@trXMaxF = maxiwqv
2064	end if
2065	else
2066	res@trXMaxF = xe
2067	end if
2068	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2069	end if
2070	end if
2071
2072	if (vNam(varn) .EQ. "wpsp") then
2073	miniwpsp=min(data(varn,:,:))
2074	maxiwpsp=max(data(varn,:,:))
2075	if (over .EQ. 1) then
2076	res@xyDashPattern = 0
2077	plot_wpsp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2078	else
2079	res@gsnLeftString = vNam(varn)
2080	res@gsnRightString = unit(varn)
2081	if (.not. isvar("xs")) then
2082	if (parameter(63) .NE. "x0") then
2083	res@trXMinF = stringtofloat(parameter(63))
2084	else
2085	res@trXMinF = miniwpsp
2086	end if
2087	else
2088	res@trXMinF = xs
2089	end if
2090	if (.not. isvar("xe")) then
2091	if (parameter(65) .NE. "xdim") then
2092	res@trXMaxF = stringtofloat(parameter(65))
2093	else
2094	res@trXMaxF = maxiwpsp
2095	end if
2096	else
2097	res@trXMaxF = xe
2098	end if
2099	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2100	end if
2101	end if
2102	if (vNam(varn) .EQ. "wsss") then
2103	miniwsss=min(data(varn,:,:))
2104	maxiwsss=max(data(varn,:,:))
2105	if (over .EQ. 1) then
2106	res@xyDashPattern = 1
2107	plot_wsss = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2108	else
2109	res@gsnLeftString = vNam(varn)
2110	res@gsnRightString = unit(varn)
2111	if (.not. isvar("xs")) then
2112	if (parameter(63) .NE. "x0") then
2113	res@trXMinF = stringtofloat(parameter(63))
2114	else
2115	res@trXMinF = miniwsss
2116	end if
2117	else
2118	res@trXMinF = xs
2119	end if
2120	if (.not. isvar("xe")) then
2121	if (parameter(65) .NE. "xdim") then
2122	res@trXMaxF = stringtofloat(parameter(65))
2123	else
2124	res@trXMaxF = maxiwsss
2125	end if
2126	else
2127	res@trXMaxF = xe
2128	end if
2129	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2130	end if
2131	end if
2132	if (vNam(varn) .EQ. "ws") then
2133	miniws=min(data(varn,:,:))
2134	maxiws=max(data(varn,:,:))
2135	if (over .EQ. 1) then
2136	res@xyDashPattern = 2
2137	plot_ws = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2138	else
2139	res@gsnLeftString = vNam(varn)
2140	res@gsnRightString = unit(varn)
2141	if (.not. isvar("xs")) then
2142	if (parameter(63) .NE. "x0") then
2143	res@trXMinF = stringtofloat(parameter(63))
2144	else
2145	res@trXMinF = miniws
2146	end if
2147	else
2148	res@trXMinF = xs
2149	end if
2150	if (.not. isvar("xe")) then
2151	if (parameter(65) .NE. "xdim") then
2152	res@trXMaxF = stringtofloat(parameter(65))
2153	else
2154	res@trXMaxF = maxiws
2155	end if
2156	else
2157	res@trXMaxF = xe
2158	end if
2159	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2160	end if
2161	end if
2162
2163	if (vNam(varn) .EQ. "wpsap") then
2164	miniwpsap=min(data(varn,:,:))
2165	maxiwpsap=max(data(varn,:,:))
2166	if (over .EQ. 1) then
2167	res@xyDashPattern = 0
2168	plot_wpsap = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2169	else
2170	res@gsnLeftString = vNam(varn)
2171	res@gsnRightString = unit(varn)
2172	if (.not. isvar("xs")) then
2173	if (parameter(63) .NE. "x0") then
2174	res@trXMinF = stringtofloat(parameter(63))
2175	else
2176	res@trXMinF = miniwpsap
2177	end if
2178	else
2179	res@trXMinF = xs
2180	end if
2181	if (.not. isvar("xe")) then
2182	if (parameter(65) .NE. "xdim") then
2183	res@trXMaxF = stringtofloat(parameter(65))
2184	else
2185	res@trXMaxF = maxiwpsap
2186	end if
2187	else
2188	res@trXMaxF = xe
2189	end if
2190	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2191	end if
2192	end if
2193	if (vNam(varn) .EQ. "wssas") then
2194	miniwssas=min(data(varn,:,:))
2195	maxiwssas=max(data(varn,:,:))
2196	if (over .EQ. 1) then
2197	res@xyDashPattern = 1
2198	plot_wssas = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2199	else
2200	res@gsnLeftString = vNam(varn)
2201	res@gsnRightString = unit(varn)
2202	if (.not. isvar("xs")) then
2203	if (parameter(63) .NE. "x0") then
2204	res@trXMinF = stringtofloat(parameter(63))
2205	else
2206	res@trXMinF = miniwssas
2207	end if
2208	else
2209	res@trXMinF = xs
2210	end if
2211	if (.not. isvar("xe")) then
2212	if (parameter(65) .NE. "xdim") then
2213	res@trXMaxF = stringtofloat(parameter(65))
2214	else
2215	res@trXMaxF = maxiwssas
2216	end if
2217	else
2218	res@trXMaxF = xe
2219	end if
2220	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2221	end if
2222	end if
2223	if (vNam(varn) .EQ. "wsa") then
2224	miniwsa=min(data(varn,:,:))
2225	maxiwsa=max(data(varn,:,:))
2226	if (over .EQ. 1) then
2227	res@xyDashPattern = 2
2228	plot_wsa = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2229	else
2230	res@gsnLeftString = vNam(varn)
2231	res@gsnRightString = unit(varn)
2232	if (.not. isvar("xs")) then
2233	if (parameter(63) .NE. "x0") then
2234	res@trXMinF = stringtofloat(parameter(63))
2235	else
2236	res@trXMinF = miniwsa
2237	end if
2238	else
2239	res@trXMinF = xs
2240	end if
2241	if (.not. isvar("xe")) then
2242	if (parameter(65) .NE. "xdim") then
2243	res@trXMaxF = stringtofloat(parameter(65))
2244	else
2245	res@trXMaxF = maxiwsa
2246	end if
2247	else
2248	res@trXMaxF = xe
2249	end if
2250	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2251	end if
2252	end if
2253
2254	if (vNam(varn) .EQ. "us2") then
2255	minius2=min(data(varn,:,:))
2256	maxius2=max(data(varn,:,:))
2257	if (over .EQ. 1) then
2258	res@xyDashPattern = 0
2259	plot_us2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2260	else
2261	res@gsnLeftString = vNam(varn)
2262	res@gsnRightString = unit(varn)
2263	if (.not. isvar("xs")) then
2264	if (parameter(63) .NE. "x0") then
2265	res@trXMinF = stringtofloat(parameter(63))
2266	else
2267	res@trXMinF = minius2
2268	end if
2269	else
2270	res@trXMinF = xs
2271	end if
2272	if (.not. isvar("xe")) then
2273	if (parameter(65) .NE. "xdim") then
2274	res@trXMaxF = stringtofloat(parameter(65))
2275	else
2276	res@trXMaxF = maxius2
2277	end if
2278	else
2279	res@trXMaxF = xe
2280	end if
2281	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2282	end if
2283	end if
2284	if (vNam(varn) .EQ. "vs2") then
2285	minivs2=min(data(varn,:,:))
2286	maxivs2=max(data(varn,:,:))
2287	if (over .EQ. 1) then
2288	res@xyDashPattern = 1
2289	plot_vs2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2290	else
2291	res@gsnLeftString = vNam(varn)
2292	res@gsnRightString = unit(varn)
2293	if (.not. isvar("xs")) then
2294	if (parameter(63) .NE. "x0") then
2295	res@trXMinF = stringtofloat(parameter(63))
2296	else
2297	res@trXMinF = minivs2
2298	end if
2299	else
2300	res@trXMinF = xs
2301	end if
2302	if (.not. isvar("xe")) then
2303	if (parameter(65) .NE. "xdim") then
2304	res@trXMaxF = stringtofloat(parameter(65))
2305	else
2306	res@trXMaxF = maxivs2
2307	end if
2308	else
2309	res@trXMaxF = xe
2310	end if
2311	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2312	end if
2313	end if
2314	if (vNam(varn) .EQ. "ws2") then
2315	miniws2=min(data(varn,:,:))
2316	maxiws2=max(data(varn,:,:))
2317	if (over .EQ. 1) then
2318	res@xyDashPattern = 2
2319	plot_ws2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2320	else
2321	res@gsnLeftString = vNam(varn)
2322	res@gsnRightString = unit(varn)
2323	if (.not. isvar("xs")) then
2324	if (parameter(63) .NE. "x0") then
2325	res@trXMinF = stringtofloat(parameter(63))
2326	else
2327	res@trXMinF = miniws2
2328	end if
2329	else
2330	res@trXMinF = xs
2331	end if
2332	if (.not. isvar("xe")) then
2333	if (parameter(65) .NE. "xdim") then
2334	res@trXMaxF = stringtofloat(parameter(65))
2335	else
2336	res@trXMaxF = maxiws2
2337	end if
2338	else
2339	res@trXMaxF = xe
2340	end if
2341	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2342	end if
2343	end if
2344
2345	if (vNam(varn) .EQ. "wsususodz") then
2346	miniwsususodz=min(data(varn,:,:))
2347	maxiwsususodz=max(data(varn,:,:))
2348	if (over .EQ. 1) then
2349	res@xyDashPattern = 0
2350	plot_wsususodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2351	else
2352	res@gsnLeftString = vNam(varn)
2353	res@gsnRightString = unit(varn)
2354	if (.not. isvar("xs")) then
2355	if (parameter(63) .NE. "x0") then
2356	res@trXMinF = stringtofloat(parameter(63))
2357	else
2358	res@trXMinF = miniwsususodz
2359	end if
2360	else
2361	res@trXMinF = xs
2362	end if
2363	if (.not. isvar("xe")) then
2364	if (parameter(65) .NE. "xdim") then
2365	res@trXMaxF = stringtofloat(parameter(65))
2366	else
2367	res@trXMaxF = maxiwsususodz
2368	end if
2369	else
2370	res@trXMaxF = xe
2371	end if
2372	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2373	end if
2374	end if
2375	if (vNam(varn) .EQ. "wspsodz") then
2376	miniwspsodz=min(data(varn,:,:))
2377	maxiwspsodz=max(data(varn,:,:))
2378	if (over .EQ. 1) then
2379	res@xyDashPattern = 1
2380	plot_wspsodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2381	else
2382	res@gsnLeftString = vNam(varn)
2383	res@gsnRightString = unit(varn)
2384	if (.not. isvar("xs")) then
2385	if (parameter(63) .NE. "x0") then
2386	res@trXMinF = stringtofloat(parameter(63))
2387	else
2388	res@trXMinF = miniwspsodz
2389	end if
2390	else
2391	res@trXMinF = xs
2392	end if
2393	if (.not. isvar("xe")) then
2394	if (parameter(65) .NE. "xdim") then
2395	res@trXMaxF = stringtofloat(parameter(65))
2396	else
2397	res@trXMaxF = maxiwspsodz
2398	end if
2399	else
2400	res@trXMaxF = xe
2401	end if
2402	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2403	end if
2404	end if
2405	if (vNam(varn) .EQ. "wpeodz") then
2406	miniwpeodz=min(data(varn,:,:))
2407	maxiwpeodz=max(data(varn,:,:))
2408	if (over .EQ. 1) then
2409	res@xyDashPattern = 2
2410	plot_wpeodz = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2411	else
2412	res@gsnLeftString = vNam(varn)
2413	res@gsnRightString = unit(varn)
2414	if (.not. isvar("xs")) then
2415	if (parameter(63) .NE. "x0") then
2416	res@trXMinF = stringtofloat(parameter(63))
2417	else
2418	res@trXMinF = miniwpeodz
2419	end if
2420	else
2421	res@trXMinF = xs
2422	end if
2423	if (.not. isvar("xe")) then
2424	if (parameter(65) .NE. "xdim") then
2425	res@trXMaxF = stringtofloat(parameter(65))
2426	else
2427	res@trXMaxF = maxiwpeodz
2428	end if
2429	else
2430	res@trXMaxF = xe
2431	end if
2432	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2433	end if
2434	end if
2435	if (over .EQ. 0) then
2436	n=n+1
2437	end if
2438	if (prof3d .EQ. 0)then
2439	varn=varn+1
2440	end if
2441	delete(temp)
2442	end if
2443	end do
2444
2445	if (count_var .EQ. 0) then
2446	print(" ")
2447	print("Please select a variable 'var=' or use the default value")
2448	print(" ")
2449	print("Your selection '"+var+"' does not exist on the input file")
2450	print(" ")
2451	exit
2452	end if
2453
2454	if (over .EQ. 1 ) then
2455
2456	overlay(plot_u,plot_v)
2457	overlay(plot_u,plot_w)
2458	u=0
2459	overlay(plot_pt,plot_vpt)
2460	overlay(plot_pt,plot_lpt)
2461	pt=0
2462	overlay(plot_q,plot_qv)
2463	overlay(plot_q,plot_ql)
2464	q=0
2465	overlay(plot_e,plot_es)
2466	e=0
2467	overlay(plot_km,plot_kh)
2468	km=0
2469	overlay(plot_wpup,plot_wsus)
2470	overlay(plot_wpup,plot_wu)
2471	wpup=0
2472	overlay(plot_wpvp,plot_wsvs)
2473	overlay(plot_wpvp,plot_wv)
2474	wpvp=0
2475	overlay(plot_wpptp,plot_wspts)
2476	overlay(plot_wpptp,plot_wpt)
2477	wpptp=0
2478	overlay(plot_wsptsBC,plot_wptBC)
2479	wsptsBC=0
2480	overlay(plot_wpvptp,plot_wsvpts)
2481	overlay(plot_wpvptp,plot_wvpt)
2482	wpvptp=0
2483	overlay(plot_wpqp,plot_wsqs)
2484	overlay(plot_wpqp,plot_wq)
2485	wpqp=0
2486	overlay(plot_wpqvp,plot_wsqvs)
2487	overlay(plot_wpqvp,plot_wqv)
2488	wpqvp=0
2489	overlay(plot_wpsp,plot_wsss)
2490	overlay(plot_wpsp,plot_ws)
2491	wpsp=0
2492	overlay(plot_wpsap,plot_wssas)
2493	overlay(plot_wpsap,plot_wsa)
2494	wpsap=0
2495	overlay(plot_us2,plot_vs2)
2496	overlay(plot_us2,plot_ws2)
2497	us2=0
2498	overlay(plot_wsususodz,plot_wspsodz)
2499	overlay(plot_wsususodz,plot_wpeodz)
2500	wsususodz=0
2501
2502	end if
2503
2504	if (over .EQ. 1) then
2505
2506	do varn = 0,dim-1
2507
2508	check = True
2509
2510	if (prof3d .EQ. 0) then
2511	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
2512	check = False
2513	end if
2514	else
2515	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "zusi") .OR. isStrSubset( vNam(varn), "zwwi") .OR. isStrSubset( vNam(varn), "x") .OR. isStrSubset( vNam(varn), "xu") .OR. isStrSubset( vNam(varn), "y") .OR. isStrSubset( vNam(varn), "yv") .OR. isStrSubset( vNam(varn), "zu_3d") .OR. isStrSubset( vNam(varn), "zw_3d")) then
2516	check = False
2517	end if
2518	end if
2519
2520	if (.not. isvar("var")) then
2521	if (parameter(21) .NE. "variables") then
2522	var=parameter(21)
2523	check = isStrSubset( var,","+vNam(varn)+"," )
2524	end if
2525	else
2526	check = isStrSubset( var,","+vNam(varn)+"," )
2527	end if
2528
2529	if (check)
2530
2531	if (prof3d .EQ. 1) then
2532
2533	end if
2534
2535	if (prof3d .EQ. 0) then
2536	z = f->$vNam(varn+1)$
2537	else
2538	z = z_v(varn,:)
2539	end if
2540
2541	res@gsnLeftString = vNam(varn)
2542	res@gsnRightString = unit(varn)
2543	res@trYMinF = min_z
2544	res@trYMaxF = max_z
2545	if (.not. isvar("xs")) then
2546	if (parameter(63) .NE. "x0") then
2547	res@trXMinF = stringtofloat(parameter(63))
2548	else
2549	res@trXMinF = min(data(varn,:,:))
2550	end if
2551	else
2552	res@trXMinF = xs
2553	end if
2554	if (.not. isvar("xe")) then
2555	if (parameter(65) .NE. "xdim") then
2556	res@trXMaxF = stringtofloat(parameter(65))
2557	else
2558	res@trXMaxF = max(data(varn,:,:))
2559	end if
2560	else
2561	res@trXMaxF = xe
2562	end if
2563	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res)
2564
2565	if (vNam(varn) .EQ. "u" .OR. vNam(varn) .EQ. "v" .OR. vNam(varn) .EQ. "w") then
2566	if (u .EQ. 0) then
2567	res@gsnLeftString = "u, v and w"
2568	res@gsnRightString = unit(varn)
2569	if (.not. isvar("xs")) then
2570	if (parameter(63) .NE. "x0") then
2571	res@trXMinF = stringtofloat(parameter(63))
2572	else
2573	res@trXMinF = min((/miniu,miniv,miniw/))
2574	end if
2575	else
2576	res@trXMinF = xs
2577	end if
2578	if (.not. isvar("xe")) then
2579	if (parameter(65) .NE. "xdim") then
2580	res@trXMaxF = stringtofloat(parameter(65))
2581	else
2582	res@trXMaxF = max((/maxiu,maxiv,maxiw/))
2583	end if
2584	else
2585	res@trXMaxF = xe
2586	end if
2587	if (vNam(varn) .EQ. "w") then
2588	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2589	else
2590	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2591	end if
2592
2593	; ***************************************************
2594	; legend for combined plot
2595	; ***************************************************
2596
2597	lgres = True
2598	lgMonoDashIndex = False
2599	lgres@lgLabelFont = "helvetica"
2600	lgres@lgLabelFontHeightF = .1
2601	lgres@vpWidthF = 0.12
2602	lgres@vpHeightF = 0.1
2603	lgres@lgDashIndexes = (/0,1,2/)
2604	lbid = gsn_create_legend(wks,3,(/"u","v","w"/),lgres)
2605
2606	amres = True
2607	amres@amParallelPosF = 0.65
2608	amres@amOrthogonalPosF = -0.2
2609	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2610	overlay(plot(n),plot_u)
2611	u=1
2612	else
2613	if (prof3d .EQ. 0)then
2614	varn=varn+1
2615	end if
2616	continue
2617	end if
2618	end if
2619
2620	if (vNam(varn) .EQ. "pt" .OR. vNam(varn) .EQ. "vpt" .OR. vNam(varn) .EQ. "lpt") then
2621	if (pt .EQ. 0) then
2622	res@gsnLeftString = "pt, vpt and lpt"
2623	res@gsnRightString = unit(varn)
2624	if (.not. isvar("xs")) then
2625	if (parameter(63) .NE. "x0") then
2626	res@trXMinF = stringtofloat(parameter(63))
2627	else
2628	res@trXMinF = min((/minipt,minivpt,minilpt/))
2629	end if
2630	else
2631	res@trXMinF = xs
2632	end if
2633	if (.not. isvar("xe")) then
2634	if (parameter(65) .NE. "xdim") then
2635	res@trXMaxF = stringtofloat(parameter(65))
2636	else
2637	res@trXMaxF = max((/maxipt,maxivpt,maxilpt/))
2638	end if
2639	else
2640	res@trXMaxF = xe
2641	end if
2642
2643	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2644
2645	; ***************************************************
2646	; legend for combined plot
2647	; ***************************************************
2648
2649	lgres = True
2650	lgMonoDashIndex = False
2651	lgres@lgLabelFont = "helvetica"
2652	lgres@lgLabelFontHeightF = .1
2653	lgres@vpWidthF = 0.12
2654	lgres@vpHeightF = 0.1
2655	lgres@lgDashIndexes = (/0,1,2/)
2656	lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres)
2657
2658	amres = True
2659	amres@amParallelPosF = 0.65
2660	amres@amOrthogonalPosF = -0.2
2661	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2662	overlay(plot(n),plot_pt)
2663	pt=1
2664	else
2665	if (prof3d .EQ. 0)then
2666	varn=varn+1
2667	end if
2668	continue
2669	end if
2670	end if
2671	if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then
2672	if (q .EQ. 0) then
2673	res@gsnLeftString = "q, qv and ql"
2674	res@gsnRightString = unit(varn)
2675	if (.not. isvar("xs")) then
2676	if (parameter(63) .NE. "x0") then
2677	res@trXMinF = stringtofloat(parameter(63))
2678	else
2679	res@trXMinF = min((/miniq,miniqv,miniql/))
2680	end if
2681	else
2682	res@trXMinF = xs
2683	end if
2684	if (.not. isvar("xe")) then
2685	if (parameter(65) .NE. "xdim") then
2686	res@trXMaxF = stringtofloat(parameter(65))
2687	else
2688	res@trXMaxF = max((/maxiq,maxiqv,maxiql/))
2689	end if
2690	else
2691	res@trXMaxF = xe
2692	end if
2693
2694	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2695
2696	; ***************************************************
2697	; legend for combined plot
2698	; ***************************************************
2699
2700	lgres = True
2701	lgMonoDashIndex = False
2702	lgres@lgLabelFont = "helvetica"
2703	lgres@lgLabelFontHeightF = .1
2704	lgres@vpWidthF = 0.12
2705	lgres@vpHeightF = 0.1
2706	lgres@lgDashIndexes = (/0,1,2/)
2707	lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres)
2708
2709	amres = True
2710	amres@amParallelPosF = 0.65
2711	amres@amOrthogonalPosF = -0.2
2712	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2713	overlay(plot(n),plot_q)
2714	q=1
2715	else
2716	if (prof3d .EQ. 0)then
2717	varn=varn+1
2718	end if
2719	continue
2720	end if
2721	end if
2722
2723	if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then
2724	if (e .EQ. 0) then
2725	res@gsnLeftString = "e and es"
2726	res@gsnRightString = unit(varn)
2727	if (.not. isvar("xs")) then
2728	if (parameter(63) .NE. "x0") then
2729	res@trXMinF = stringtofloat(parameter(63))
2730	else
2731	res@trXMinF = min((/minie,minies/))
2732	end if
2733	else
2734	res@trXMinF = xs
2735	end if
2736	if (.not. isvar("xe")) then
2737	if (parameter(65) .NE. "xdim") then
2738	res@trXMaxF = stringtofloat(parameter(65))
2739	else
2740	res@trXMaxF = max((/maxie,maxies/))
2741	end if
2742	else
2743	res@trXMaxF = xe
2744	end if
2745
2746	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2747
2748	; ***************************************************
2749	; legend for combined plot
2750	; ***************************************************
2751
2752	lgres = True
2753	lgMonoDashIndex = False
2754	lgres@lgLabelFont = "helvetica"
2755	lgres@lgLabelFontHeightF = .1
2756	lgres@vpWidthF = 0.12
2757	lgres@vpHeightF = 0.1
2758	lgres@lgDashIndexes = (/0,1,2/)
2759	lbid = gsn_create_legend(wks,3,(/"e","es"/),lgres)
2760
2761	amres = True
2762	amres@amParallelPosF = 0.65
2763	amres@amOrthogonalPosF = -0.2
2764	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2765	overlay(plot(n),plot_e)
2766	e=1
2767	else
2768	if (prof3d .EQ. 0)then
2769	varn=varn+1
2770	end if
2771	continue
2772	end if
2773	end if
2774	if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then
2775	if (km .EQ. 0) then
2776	res@gsnLeftString = "km and kh"
2777	res@gsnRightString = unit(varn)
2778	if (.not. isvar("xs")) then
2779	if (parameter(63) .NE. "x0") then
2780	res@trXMinF = stringtofloat(parameter(63))
2781	else
2782	res@trXMinF = min((/minikm,minikh/))
2783	end if
2784	else
2785	res@trXMinF = xs
2786	end if
2787	if (.not. isvar("xe")) then
2788	if (parameter(65) .NE. "xdim") then
2789	res@trXMaxF = stringtofloat(parameter(65))
2790	else
2791	res@trXMaxF = max((/maxikm,maxikh/))
2792	end if
2793	else
2794	res@trXMaxF = xe
2795	end if
2796
2797	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2798
2799	; ***************************************************
2800	; legend for combined plot
2801	; ***************************************************
2802
2803	lgres = True
2804	lgMonoDashIndex = False
2805	lgres@lgLabelFont = "helvetica"
2806	lgres@lgLabelFontHeightF = .1
2807	lgres@vpWidthF = 0.12
2808	lgres@vpHeightF = 0.1
2809	lgres@lgDashIndexes = (/0,1,2/)
2810	lbid = gsn_create_legend(wks,3,(/"km","kh"/),lgres)
2811
2812	amres = True
2813	amres@amParallelPosF = 0.65
2814	amres@amOrthogonalPosF = -0.2
2815	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2816	overlay(plot(n),plot_km)
2817	km=1
2818	else
2819	if (prof3d .EQ. 0)then
2820	varn=varn+1
2821	end if
2822	continue
2823	end if
2824	end if
2825
2826	if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then
2827	if (wpup .EQ. 0) then
2828	res@gsnLeftString = "wpup, wsus and wu"
2829	res@gsnRightString = unit(varn)
2830	if (.not. isvar("xs")) then
2831	if (parameter(63) .NE. "x0") then
2832	res@trXMinF = stringtofloat(parameter(63))
2833	else
2834	res@trXMinF = min((/miniwpup,miniwsus,miniwu/))
2835	end if
2836	else
2837	res@trXMinF = xs
2838	end if
2839	if (.not. isvar("xe")) then
2840	if (parameter(65) .NE. "xdim") then
2841	res@trXMaxF = stringtofloat(parameter(65))
2842	else
2843	res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/))
2844	end if
2845	else
2846	res@trXMaxF = xe
2847	end if
2848
2849	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2850
2851	; ***************************************************
2852	; legend for combined plot
2853	; ***************************************************
2854
2855	lgres = True
2856	lgMonoDashIndex = False
2857	lgres@lgLabelFont = "helvetica"
2858	lgres@lgLabelFontHeightF = .1
2859	lgres@vpWidthF = 0.12
2860	lgres@vpHeightF = 0.1
2861	lgres@lgDashIndexes = (/0,1,2/)
2862	lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres)
2863
2864	amres = True
2865	amres@amParallelPosF = 0.65
2866	amres@amOrthogonalPosF = -0.2
2867	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2868	overlay(plot(n),plot_wpup)
2869	wpup=1
2870	else
2871	if (prof3d .EQ. 0)then
2872	varn=varn+1
2873	end if
2874	continue
2875	end if
2876	end if
2877	if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then
2878	if (wpvp .EQ. 0) then
2879	res@gsnLeftString = "wpvp, wsus and wv"
2880	res@gsnRightString = unit(varn)
2881	if (.not. isvar("xs")) then
2882	if (parameter(63) .NE. "x0") then
2883	res@trXMinF = stringtofloat(parameter(63))
2884	else
2885	res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/))
2886	end if
2887	else
2888	res@trXMinF = xs
2889	end if
2890	if (.not. isvar("xe")) then
2891	if (parameter(65) .NE. "xdim") then
2892	res@trXMaxF = stringtofloat(parameter(65))
2893	else
2894	res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/))
2895	end if
2896	else
2897	res@trXMaxF = xe
2898	end if
2899
2900	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2901
2902	; ***************************************************
2903	; legend for combined plot
2904	; ***************************************************
2905
2906	lgres = True
2907	lgMonoDashIndex = False
2908	lgres@lgLabelFont = "helvetica"
2909	lgres@lgLabelFontHeightF = .1
2910	lgres@vpWidthF = 0.12
2911	lgres@vpHeightF = 0.1
2912	lgres@lgDashIndexes = (/0,1,2/)
2913	lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres)
2914
2915	amres = True
2916	amres@amParallelPosF = 0.65
2917	amres@amOrthogonalPosF = -0.2
2918	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2919	overlay(plot(n),plot_wpvp)
2920	wpup=0
2921	else
2922	if (prof3d .EQ. 0)then
2923	varn=varn+1
2924	end if
2925	continue
2926	end if
2927	end if
2928	if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then
2929	if (wpptp .EQ. 0) then
2930	res@gsnLeftString = "wpptp, wspts and wv"
2931	res@gsnRightString = unit(varn)
2932	if (.not. isvar("xs")) then
2933	if (parameter(63) .NE. "x0") then
2934	res@trXMinF = stringtofloat(parameter(63))
2935	else
2936	res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/))
2937	end if
2938	else
2939	res@trXMinF = xs
2940	end if
2941	if (.not. isvar("xe")) then
2942	if (parameter(65) .NE. "xdim") then
2943	res@trXMaxF = stringtofloat(parameter(65))
2944	else
2945	res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/))
2946	end if
2947	else
2948	res@trXMaxF = xe
2949	end if
2950
2951	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2952
2953	; ***************************************************
2954	; legend for combined plot
2955	; ***************************************************
2956
2957	lgres = True
2958	lgMonoDashIndex = False
2959	lgres@lgLabelFont = "helvetica"
2960	lgres@lgLabelFontHeightF = .1
2961	lgres@vpWidthF = 0.12
2962	lgres@vpHeightF = 0.1
2963	lgres@lgDashIndexes = (/0,1,2/)
2964	lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres)
2965
2966	amres = True
2967	amres@amParallelPosF = 0.65
2968	amres@amOrthogonalPosF = -0.2
2969	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
2970	overlay(plot(n),plot_wpptp)
2971	wpptp=1
2972	else
2973	if (prof3d .EQ. 0)then
2974	varn=varn+1
2975	end if
2976	continue
2977	end if
2978	end if
2979	if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then
2980	if (wsptsBC .EQ. 0) then
2981	res@gsnLeftString = "wsptsBC and wptBC"
2982	res@gsnRightString = unit(varn)
2983	if (.not. isvar("xs")) then
2984	if (parameter(63) .NE. "x0") then
2985	res@trXMinF = stringtofloat(parameter(63))
2986	else
2987	res@trXMinF = min((/miniwsptsBC,miniwptBC/))
2988	end if
2989	else
2990	res@trXMinF = xs
2991	end if
2992	if (.not. isvar("xe")) then
2993	if (parameter(65) .NE. "xdim") then
2994	res@trXMaxF = stringtofloat(parameter(65))
2995	else
2996	res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/))
2997	end if
2998	else
2999	res@trXMaxF = xe
3000	end if
3001
3002	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3003
3004	; ***************************************************
3005	; legend for combined plot
3006	; ***************************************************
3007
3008	lgres = True
3009	lgMonoDashIndex = False
3010	lgres@lgLabelFont = "helvetica"
3011	lgres@lgLabelFontHeightF = .1
3012	lgres@vpWidthF = 0.12
3013	lgres@vpHeightF = 0.1
3014	lgres@lgDashIndexes = (/0,1,2/)
3015	lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres)
3016
3017	amres = True
3018	amres@amParallelPosF = 0.65
3019	amres@amOrthogonalPosF = -0.2
3020	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3021	overlay(plot(n),plot_wsptsBC)
3022	wsptsBC=1
3023	else
3024	if (prof3d .EQ. 0)then
3025	varn=varn+1
3026	end if
3027	continue
3028	end if
3029	end if
3030	if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then
3031	if (wpvptp .EQ. 0) then
3032	res@gsnLeftString = "wpvptp, wsvpts and wv"
3033	res@gsnRightString = unit(varn)
3034	if (.not. isvar("xs")) then
3035	if (parameter(63) .NE. "x0") then
3036	res@trXMinF = stringtofloat(parameter(63))
3037	else
3038	res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/))
3039	end if
3040	else
3041	res@trXMinF = xs
3042	end if
3043	if (.not. isvar("xe")) then
3044	if (parameter(65) .NE. "xdim") then
3045	res@trXMaxF = stringtofloat(parameter(65))
3046	else
3047	res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/))
3048	end if
3049	else
3050	res@trXMaxF = xe
3051	end if
3052
3053	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3054
3055	; ***************************************************
3056	; legend for combined plot
3057	; ***************************************************
3058
3059	lgres = True
3060	lgMonoDashIndex = False
3061	lgres@lgLabelFont = "helvetica"
3062	lgres@lgLabelFontHeightF = .1
3063	lgres@vpWidthF = 0.12
3064	lgres@vpHeightF = 0.1
3065	lgres@lgDashIndexes = (/0,1,2/)
3066	lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres)
3067
3068	amres = True
3069	amres@amParallelPosF = 0.65
3070	amres@amOrthogonalPosF = -0.2
3071	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3072	overlay(plot(n),plot_wpvptp)
3073	wpvptp=1
3074	else
3075	if (prof3d .EQ. 0)then
3076	varn=varn+1
3077	end if
3078	continue
3079	end if
3080	end if
3081	if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then
3082	if (wpqp .EQ. 0) then
3083	res@gsnLeftString = "wpqp, wsqs and wq"
3084	res@gsnRightString = unit(varn)
3085	if (.not. isvar("xs")) then
3086	if (parameter(63) .NE. "x0") then
3087	res@trXMinF = stringtofloat(parameter(63))
3088	else
3089	res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/))
3090	end if
3091	else
3092	res@trXMinF = xs
3093	end if
3094	if (.not. isvar("xe")) then
3095	if (parameter(65) .NE. "xdim") then
3096	res@trXMaxF = stringtofloat(parameter(65))
3097	else
3098	res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/))
3099	end if
3100	else
3101	res@trXMaxF = xe
3102	end if
3103
3104	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3105
3106	; ***************************************************
3107	; legend for combined plot
3108	; ***************************************************
3109
3110	lgres = True
3111	lgMonoDashIndex = False
3112	lgres@lgLabelFont = "helvetica"
3113	lgres@lgLabelFontHeightF = .1
3114	lgres@vpWidthF = 0.12
3115	lgres@vpHeightF = 0.1
3116	lgres@lgDashIndexes = (/0,1,2/)
3117	lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres)
3118
3119	amres = True
3120	amres@amParallelPosF = 0.65
3121	amres@amOrthogonalPosF = -0.2
3122	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3123	overlay(plot(n),plot_wpqp)
3124	wpqp=1
3125	else
3126	if (prof3d .EQ. 0)then
3127	varn=varn+1
3128	end if
3129	continue
3130	end if
3131	end if
3132	if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then
3133	if (wpqvp .EQ. 0) then
3134	res@gsnLeftString = "wpqvp, wsqvs and wqv"
3135	res@gsnRightString = unit(varn)
3136	if (.not. isvar("xs")) then
3137	if (parameter(63) .NE. "x0") then
3138	res@trXMinF = stringtofloat(parameter(63))
3139	else
3140	res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/))
3141	end if
3142	else
3143	res@trXMinF = xs
3144	end if
3145	if (.not. isvar("xe")) then
3146	if (parameter(65) .NE. "xdim") then
3147	res@trXMaxF = stringtofloat(parameter(65))
3148	else
3149	res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/))
3150	end if
3151	else
3152	res@trXMaxF = xe
3153	end if
3154
3155	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3156
3157	; ***************************************************
3158	; legend for combined plot
3159	; ***************************************************
3160
3161	lgres = True
3162	lgMonoDashIndex = False
3163	lgres@lgLabelFont = "helvetica"
3164	lgres@lgLabelFontHeightF = .1
3165	lgres@vpWidthF = 0.12
3166	lgres@vpHeightF = 0.1
3167	lgres@lgDashIndexes = (/0,1,2/)
3168	lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres)
3169
3170	amres = True
3171	amres@amParallelPosF = 0.65
3172	amres@amOrthogonalPosF = -0.2
3173	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3174	overlay(plot(n),plot_wpqvp)
3175	wpqvp=1
3176	else
3177	if (prof3d .EQ. 0)then
3178	varn=varn+1
3179	end if
3180	continue
3181	end if
3182	end if
3183	if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then
3184	if (wpsp .EQ. 0) then
3185	res@gsnLeftString = "wpsp, wsss and ws"
3186	res@gsnRightString = unit(varn)
3187	if (.not. isvar("xs")) then
3188	if (parameter(63) .NE. "x0") then
3189	res@trXMinF = stringtofloat(parameter(63))
3190	else
3191	res@trXMinF = min((/miniwpsp,miniwsss,miniws/))
3192	end if
3193	else
3194	res@trXMinF = xs
3195	end if
3196	if (.not. isvar("xe")) then
3197	if (parameter(65) .NE. "xdim") then
3198	res@trXMaxF = stringtofloat(parameter(65))
3199	else
3200	res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/))
3201	end if
3202	else
3203	res@trXMaxF = xe
3204	end if
3205
3206	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3207
3208	; ***************************************************
3209	; legend for combined plot
3210	; ***************************************************
3211
3212	lgres = True
3213	lgMonoDashIndex = False
3214	lgres@lgLabelFont = "helvetica"
3215	lgres@lgLabelFontHeightF = .1
3216	lgres@vpWidthF = 0.12
3217	lgres@vpHeightF = 0.1
3218	lgres@lgDashIndexes = (/0,1,2/)
3219	lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres)
3220
3221	amres = True
3222	amres@amParallelPosF = 0.65
3223	amres@amOrthogonalPosF = -0.2
3224	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3225	overlay(plot(n),plot_wpsp)
3226	wpsp=1
3227	else
3228	if (prof3d .EQ. 0)then
3229	varn=varn+1
3230	end if
3231	continue
3232	end if
3233	end if
3234	if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then
3235	if (wpsap .EQ. 0) then
3236	res@gsnLeftString = "wpsap, wssas and wsa"
3237	res@gsnRightString = unit(varn)
3238	if (.not. isvar("xs")) then
3239	if (parameter(63) .NE. "x0") then
3240	res@trXMinF = stringtofloat(parameter(63))
3241	else
3242	res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/))
3243	end if
3244	else
3245	res@trXMinF = xs
3246	end if
3247	if (.not. isvar("xe")) then
3248	if (parameter(65) .NE. "xdim") then
3249	res@trXMaxF = stringtofloat(parameter(65))
3250	else
3251	res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/))
3252	end if
3253	else
3254	res@trXMaxF = xe
3255	end if
3256
3257	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3258
3259	; ***************************************************
3260	; legend for combined plot
3261	; ***************************************************
3262
3263	lgres = True
3264	lgMonoDashIndex = False
3265	lgres@lgLabelFont = "helvetica"
3266	lgres@lgLabelFontHeightF = .1
3267	lgres@vpWidthF = 0.12
3268	lgres@vpHeightF = 0.1
3269	lgres@lgDashIndexes = (/0,1,2/)
3270	lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres)
3271
3272	amres = True
3273	amres@amParallelPosF = 0.65
3274	amres@amOrthogonalPosF = -0.2
3275	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3276	overlay(plot(n),plot_wpsap)
3277	wpsap=1
3278	else
3279	if (prof3d .EQ. 0)then
3280	varn=varn+1
3281	end if
3282	continue
3283	end if
3284	end if
3285
3286	if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then
3287	if (us2 .EQ. 0) then
3288	res@gsnLeftString = "us2, vs2 and ws2"
3289	res@gsnRightString = unit(varn)
3290	if (.not. isvar("xs")) then
3291	if (parameter(63) .NE. "x0") then
3292	res@trXMinF = stringtofloat(parameter(63))
3293	else
3294	res@trXMinF = min((/minius2,minivs2,miniws2/))
3295	end if
3296	else
3297	res@trXMinF = xs
3298	end if
3299	if (.not. isvar("xe")) then
3300	if (parameter(65) .NE. "xdim") then
3301	res@trXMaxF = stringtofloat(parameter(65))
3302	else
3303	res@trXMaxF = max((/maxius2,maxivs2,maxiws2/))
3304	end if
3305	else
3306	res@trXMaxF = xe
3307	end if
3308
3309	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3310
3311	; ***************************************************
3312	; legend for combined plot
3313	; ***************************************************
3314
3315	lgres = True
3316	lgMonoDashIndex = False
3317	lgres@lgLabelFont = "helvetica"
3318	lgres@lgLabelFontHeightF = .1
3319	lgres@vpWidthF = 0.12
3320	lgres@vpHeightF = 0.1
3321	lgres@lgDashIndexes = (/0,1,2/)
3322	lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres)
3323
3324	amres = True
3325	amres@amParallelPosF = 0.65
3326	amres@amOrthogonalPosF = -0.2
3327	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3328	overlay(plot(n),plot_us2)
3329	us2=1
3330	else
3331	if (prof3d .EQ. 0)then
3332	varn=varn+1
3333	end if
3334	continue
3335	end if
3336	end if
3337
3338	if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then
3339	if (wsususodz .EQ. 0) then
3340	res@gsnLeftString = "wsususodz, wspsodz and ws2"
3341	res@gsnRightString = unit(varn)
3342	if (.not. isvar("xs")) then
3343	if (parameter(63) .NE. "x0") then
3344	res@trXMinF = stringtofloat(parameter(63))
3345	else
3346	res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/))
3347	end if
3348	else
3349	res@trXMinF = xs
3350	end if
3351	if (.not. isvar("xe")) then
3352	if (parameter(65) .NE. "xdim") then
3353	res@trXMaxF = stringtofloat(parameter(65))
3354	else
3355	res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/))
3356	end if
3357	else
3358	res@trXMaxF = xe
3359	end if
3360
3361	plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3362
3363	; ***************************************************
3364	; legend for combined plot
3365	; ***************************************************
3366
3367	lgres = True
3368	lgMonoDashIndex = False
3369	lgres@lgLabelFont = "helvetica"
3370	lgres@lgLabelFontHeightF = .1
3371	lgres@vpWidthF = 0.12
3372	lgres@vpHeightF = 0.1
3373	lgres@lgDashIndexes = (/0,1,2/)
3374	lbid = gsn_create_legend(wks,3,(/"wsususodz","wspsodz","wpeodz"/),lgres)
3375
3376	amres = True
3377	amres@amParallelPosF = 0.65
3378	amres@amOrthogonalPosF = -0.2
3379	annoid1 = gsn_add_annotation(plot(n),lbid,amres)
3380	overlay(plot(n),plot_wsususodz)
3381	wsususodz=1
3382	else
3383	if (prof3d .EQ. 0)then
3384	varn=varn+1
3385	end if
3386	continue
3387	end if
3388	end if
3389	n=n+1
3390	if (prof3d .EQ. 0)then
3391	varn=varn+1
3392	end if
3393	end if
3394	end do
3395	end if
3396
3397	if (combine .EQ. 1) then
3398	co=0
3399	n_o=0
3400	do varn = 0,dim-1
3401
3402	check = True
3403
3404	if (prof3d .EQ. 0) then
3405	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
3406	check = False
3407	end if
3408	else
3409	if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "zusi") .OR. isStrSubset( vNam(varn), "zwwi") .OR. isStrSubset( vNam(varn), "x") .OR. isStrSubset( vNam(varn), "xu") .OR. isStrSubset( vNam(varn), "y") .OR. isStrSubset( vNam(varn), "yv") .OR. isStrSubset( vNam(varn), "zu_3d") .OR. isStrSubset( vNam(varn), "zw_3d")) then
3410	check = False
3411	end if
3412	end if
3413
3414	if (.not. isvar("var")) then
3415	if (parameter(21) .NE. "variables") then
3416	var=parameter(21)
3417	check = isStrSubset( var,","+vNam(varn)+"," )
3418	end if
3419	else
3420	check = isStrSubset( var,","+vNam(varn)+"," )
3421	end if
3422
3423	if (check)
3424
3425	if (prof3d .EQ. 0) then
3426	z = f->$vNam(varn+1)$
3427	else
3428	z = z_v(varn,:)
3429	end if
3430
3431	com=isStrSubset( c_var,","+vNam(varn)+"," )
3432
3433	if (com)
3434	co = co+1
3435	if (n_o .EQ. 1) then
3436	res@xyDashPattern = 1
3437
3438	else
3439	if (n_o .EQ. 2) then
3440	res@xyDashPattern = 2
3441	else
3442	res@xyDashPattern = 0
3443	res@gsnLeftString = c_var
3444	res@gsnRightString = unit(varn)
3445	if (.not. isvar("xs")) then
3446	if (parameter(63) .NE. "x0") then
3447	res@trXMinF = stringtofloat(parameter(63))
3448	else
3449	res@trXMinF = min(mini)
3450	end if
3451	else
3452	res@trXMinF = xs
3453	end if
3454	if (.not. isvar("xe")) then
3455	if (parameter(65) .NE. "xdim") then
3456	res@trXMaxF = stringtofloat(parameter(65))
3457	else
3458	res@trXMaxF = max(maxi)
3459	end if
3460	else
3461	res@trXMaxF = xe
3462	end if
3463	end if
3464	end if
3465	label(n_o)=vNam(varn)
3466	color_o(n_o)=237
3467	plot_o(n_o)=gsn_csm_xy(wks,data(varn,:,:),z,res)
3468	n_o=n_o+1
3469	end if
3470	if (prof3d .EQ. 0)then
3471	varn=varn+1
3472	end if
3473	end if
3474	end do
3475
3476	if(number_comb .EQ. 2)then
3477	if (co .EQ. 2)then
3478	overlay(plot_o(0),plot_o(1))
3479	else
3480	print("'c_var' is not conform with existing variables on the input file")
3481	end if
3482	end if
3483	if(number_comb .EQ. 3)then
3484	if (co .EQ. 3)then
3485	overlay(plot_o(0),plot_o(1))
3486	overlay(plot_o(0),plot_o(2))
3487	else
3488	print("'c_var' is not conform with existing variables on the input file")
3489	end if
3490	end if
3491
3492	; ***************************************************
3493	; legend for combined plot
3494	; ***************************************************
3495
3496	lgres = True
3497	lgMonoDashIndex = False
3498	lgres@lgDashIndexes = (/0,1,2/)
3499	lgres@lgLabelFont = "helvetica"
3500	lgres@lgLabelFontHeightF = .1
3501	lgres@vpWidthF = 0.12
3502	lgres@vpHeightF = 0.1
3503
3504	lbid = gsn_create_legend(wks,number_comb,label,lgres)
3505
3506	amres = True
3507	amres@amParallelPosF = 0.65
3508	amres@amOrthogonalPosF = -0.2
3509	annoid1 = gsn_add_annotation(plot_o(0),lbid,amres)
3510
3511	plot(0) = plot_o(0)
3512
3513	end if
3514
3515	; ***************************************************
3516	; merge plots onto one page
3517	; ***************************************************
3518
3519	do m=0,n-1
3520	plot_(m)=plot(n-1-m)
3521	end do
3522
3523	if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
3524	gsn_panel(wks,plot_(0:n-1),(/n,1/),resP)
3525	else
3526	do i = 0,n-1, no_lines*no_columns
3527	if( (i+no_lines*no_columns) .gt. (n-1)) then
3528	gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP)
3529	else
3530	gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
3531	end if
3532	end do
3533	end if
3534
3535	print(" ")
3536	print("Output to: " + file_out +"."+ format_out)
3537	print(" ")
3538
3539	end

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

Download in other formats:

Original Format