Home

Context Navigation

profiles.ncl @ 742

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

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

Download in other formats:

Original Format