Home

Context Navigation

profiles.ncl @ 190

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

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

Download in other formats:

Original Format