Home

Context Navigation

profiles.ncl @ 981

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

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

Download in other formats:

Original Format