Home

Context Navigation

profiles.ncl @ 622

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

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

Download in other formats:

Original Format