Home

Context Navigation

profiles.ncl @ 515

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

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

Download in other formats:

Original Format