Home

Context Navigation

profiles.ncl @ 328

Last change on this file since 328 was 324, checked in by weinreis, 16 years ago
modified spectra.ncl and profiles.ncl
File size: 149.1 KB

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

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

Download in other formats:

Original Format