Home

Context Navigation

profiles.ncl @ 844

Last change on this file since 844 was 827, checked in by heinze, 13 years ago
allow plotting of data with very small time increments
File size: 171.1 KB

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

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

Download in other formats:

Original Format