Home

Context Navigation

profiles.ncl @ 389

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

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

Download in other formats:

Original Format