Home

Context Navigation

profiles.ncl @ 243

Last change on this file since 243 was 236, checked in by letzel, 15 years ago
NCL scripts in trunk/SCRIPTS/NCL updated
File size: 143.1 KB

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

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

Download in other formats:

Original Format