Home

Context Navigation

profiles.ncl @ 241

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

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

Download in other formats:

Original Format