Home

Context Navigation

profiles.ncl @ 1534

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

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

Download in other formats:

Original Format