Home

Context Navigation

profiles.ncl @ 958

Last change on this file since 958 was 958, checked in by hoffmann, 12 years ago
minor changes in palmplot pr
File size: 140.6 KB

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

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

Download in other formats:

Original Format