Home

Context Navigation

timeseries.ncl @ 566

Last change on this file since 566 was 566, checked in by heinze, 14 years ago
Formatting of all NCL scripts. Items of .ncl.config are re-sorted and xyc, xzc and yzc are no parameters any more. Parameters start_f_1/end_f_1 are renamed to start_f/end_f in profiles.ncl. Bugfix in cross_sections: enable vector plot if var is set explicitly. Deletion of NCL user guide because guide is no available online.
File size: 49.6 KB

Rev	Line
[194]	1	load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
[154]	2	load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
	3	load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
	4	load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
	5
[190]	6	;***************************************************
[418]	7	; load .ncl.config or .ncl.config.default
[190]	8	;***************************************************
[194]	9
	10	function which_script()
	11	local script
	12	begin
	13	script="timeseries"
	14	return(script)
	15	end
[418]	16
	17	if (isfilepresent("$PALM_BIN/../../.ncl.config")) then
	18	loadscript("$PALM_BIN/../../.ncl.config")
[190]	19	else
[418]	20	if (isfilepresent("$PALM_BIN/NCL/.ncl.config.default")) then
	21	loadscript( "$PALM_BIN/NCL/.ncl.config.default")
[190]	22	else
[418]	23	palm_bin_path = getenv("PALM_BIN")
[190]	24	print(" ")
[418]	25	print("Neither the personal configuration file '.ncl.config' exists in")
	26	print("~/palm/current_version")
[566]	27	print("nor the default configuration file '.ncl.config.default' "+\
	28	"exists in")
[418]	29	print(palm_bin_path + "/NCL")
[190]	30	print(" ")
	31	exit
[154]	32	end if
[418]	33	end if
[190]	34
	35	begin
[154]	36
[190]	37	;***************************************************
[532]	38	; Retrieving the NCL version used
	39	;***************************************************
	40
	41	ncl_version_ch = systemfunc("ncl -V")
	42	ncl_version = stringtofloat(ncl_version_ch)
	43
	44	;***************************************************
[526]	45	; Retrieving the double quote character
	46	;***************************************************
	47
	48	dq=str_get_dq()
	49
	50	;***************************************************
[190]	51	; set up default parameter values and strings
	52	;***************************************************
[218]	53
[190]	54	if (file_1 .EQ. "File in") then
	55	print(" ")
[418]	56	print("Declare input file 'file_1=' in '.ncl.config' or prompt")
[190]	57	print(" ")
	58	exit
[175]	59	else
	60	file_in = file_1
[154]	61	end if
[175]	62
[566]	63	if (format_out .NE. "x11" .AND. format_out .NE. "pdf" .AND. \
	64	format_out .NE. "eps" .AND. format_out .NE. "ps" .AND. \
	65	format_out .NE. "epsi" .AND. format_out .NE. "ncgm" .AND. \
	66	format_out .NE. "png")then
[190]	67	print(" ")
	68	print("'format_out = "+format_out+"' is invalid and set to'x11'")
	69	print(" ")
	70	format_out="x11"
[154]	71	end if
[175]	72
[532]	73	if (ncl_version .LE. 5.1 .AND. format_out .EQ. "png") then
	74	print(" ")
	75	print("Output of png files not available")
	76	print("png output is avaiable with NCL version 5.2.0 and higher ")
	77	print("NCL version used: " + ncl_version_ch)
	78	print(" ")
	79	exit
	80	end if
	81
[190]	82	if (over .NE. 0 .AND. over .NE. 1) then
	83	print(" ")
	84	print("'over'= "+over+" is invalid and set to 0")
	85	print(" ")
[161]	86	over = 0
[218]	87	end if
	88
	89	if (norm_t .EQ. 0) then
	90	print(" ")
	91	print("Normalising with 0 is not allowed, 'norm_t' is set to 1.0")
	92	print(" ")
	93	norm_t = 1.0
	94	end if
[190]	95
[175]	96
[190]	97	;***************************************************
[154]	98	; open input file
[190]	99	;***************************************************
[218]	100
	101	file_in_1 = False
	102	if (isStrSubset(file_in, ".nc"))then
	103	start_f = -2
	104	end_f = -2
	105	file_in_1 = True
	106	end if
[154]	107
[218]	108	if (start_f .EQ. -1)then
	109	print(" ")
[566]	110	print("'start_f' must be one of the cyclic numbers (at least 0) of "+\
	111	"your input file(s)")
[218]	112	print(" ")
	113	exit
	114	end if
	115	if (end_f .EQ. -1)then
	116	print(" ")
[566]	117	print("'end_f' must be one of the cyclic numbers (at least 0) of "+\
	118	"your input file(s)")
[218]	119	print(" ")
	120	exit
	121	end if
[154]	122
[218]	123	files=new(end_f-start_f+1,string)
	124
	125	if (file_in_1)then
	126	if (isfilepresent(file_in))then
	127	files(0)=file_in
	128	else
	129	print(" ")
	130	print("1st input file: '"+file_in+"' does not exist")
	131	print(" ")
	132	exit
	133	end if
	134	else
	135	if (start_f .EQ. 0)then
	136	if (isfilepresent(file_in+".nc"))then
	137	files(0)=file_in+".nc"
	138	do i=1,end_f
	139	if (isfilepresent(file_in+"."+i+".nc"))then
	140	files(i)=file_in+"."+i+".nc"
	141	else
	142	print(" ")
	143	print("Input file: '"+file_in+"."+i+".nc' does not exist")
	144	print(" ")
	145	exit
	146	end if
	147	end do
	148	else
	149	print(" ")
	150	print("Input file: '"+file_in+".nc' does not exist")
	151	print(" ")
	152	exit
	153	end if
	154	else
	155	do i=start_f,end_f
	156	if (isfilepresent(file_in+"."+i+".nc"))then
	157	files(i-start_f)=file_in+"."+i+".nc"
	158	else
	159	print(" ")
	160	print("Input file: '"+file_in+"."+i+".nc' does not exist")
	161	print(" ")
	162	exit
	163	end if
	164	end do
	165	end if
	166	end if
	167
	168	f=addfiles(files,"r")
	169	f_att=addfile(files(0),"r")
	170	ListSetType(f,"cat")
	171
	172	vNam=getfilevarnames(f_att)
[310]	173	vNam_static=vNam
[218]	174
[154]	175	print(" ")
[190]	176	print("Variables in input file:")
	177	print("- "+ vNam)
[154]	178	print(" ")
[161]	179	dim = dimsizes(vNam)
	180	if (dim .EQ. 0) then
	181	print(" ")
[310]	182	print("There is no data on file")
[161]	183	print(" ")
	184	end if
	185
[310]	186	t_all = f[:]->time
	187	nt = dimsizes(t_all)
	188	delta_t = t_all(nt-1)/nt
[218]	189
[190]	190	;****************************************************
[154]	191	; start of time step and different types of mistakes that could be done
[190]	192	;****************************************************
[154]	193
[190]	194	if (start_time_step .EQ. -1.) then
	195	start_time_step=t_all(0)/3600
[154]	196	else
[162]	197	if (start_time_step .GE. t_all(nt-1)/3600)
[154]	198	print(" ")
[566]	199	print("'start_time_step' = "+ start_time_step +\
	200	"h is equal or greater than last time step = " + \
	201	t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
[154]	202	print(" ")
[218]	203	print("Select another 'start_time_step'")
[162]	204	print(" ")
[154]	205	exit
	206	end if
[162]	207	if (start_time_step .LT. t_all(0)/3600)
[154]	208	print(" ")
[566]	209	print("'start_time_step' = "+ start_time_step +\
	210	"h is lower than first time step = " + t_all(0)+"s = "+\
	211	t_all(0)/3600+"h")
[154]	212	print(" ")
[218]	213	print("Select another 'start_time_step'")
[162]	214	print(" ")
[154]	215	exit
	216	end if
	217	end if
[162]	218	do i=0,nt-2
[566]	219	if (start_time_step .GE. (t_all(i)-delta_t/2)/3600 .AND. \
	220	start_time_step .LT. (t_all(i)+delta_t/2)/3600)then
[162]	221	st=i
	222	break
	223	end if
	224	end do
[566]	225	if (start_time_step .GE. t_all(nt-1)-delta_t/2 .AND. \
	226	start_time_step .LT. t_all(nt-1)) then
[162]	227	st=nt-2
	228	end if
[190]	229	if (.not. isvar("st"))then
	230	print(" ")
	231	print("'start_time_step' = "+ start_time_step +"h is invalid")
	232	print(" ")
[218]	233	print("Select another 'start_time_step'")
[190]	234	print(" ")
	235	exit
	236	end if
[162]	237
[190]	238	;****************************************************
[154]	239	; end of time step and different types of mistakes that could be done
[190]	240	;****************************************************
[154]	241
[190]	242	if (end_time_step .EQ. -1.) then
	243	end_time_step = t_all(nt-1)/3600
[154]	244	else
[162]	245	if (end_time_step .LE. t_all(0)/3600)
[154]	246	print(" ")
[566]	247	print("'end_time_step' = "+end_time_step+ \
	248	"h is lower or equal than first time step = " + \
	249	t_all(0)+"s = "+t_all(0)/3600+"h")
[154]	250	print(" ")
[218]	251	print("Select another 'end_time_step'")
[162]	252	print(" ")
[154]	253	exit
	254	end if
[162]	255	if (end_time_step .GT. t_all(nt-1)/3600)
[154]	256	print(" ")
[566]	257	print("'end_time_step' = "+ end_time_step +\
	258	"h is greater than last time step = " + t_all(nt-1)+"s = "+\
	259	t_all(nt-1)/3600+"h")
[154]	260	print(" ")
[218]	261	print("Select another 'end_time_step'")
[162]	262	print(" ")
[154]	263	exit
	264	end if
[162]	265	if (end_time_step .LE. start_time_step/3600)
[161]	266	print(" ")
[566]	267	print("'end_time_step' = "+ end_time_step +\
	268	"h is equal or lower than 'start_time_step' = "+\
	269	start_time_step+"h")
[161]	270	print(" ")
[218]	271	print("Select another 'start_time_step' or 'end_time_step'")
[162]	272	print(" ")
[161]	273	exit
[154]	274	end if
[162]	275	end if
	276	do i=0,nt-1
[566]	277	if (end_time_step .GE. (t_all(i)-delta_t/2)/3600 .AND. \
	278	end_time_step .LT. (t_all(i)+delta_t/2)/3600)then
[162]	279	et=i
	280	break
	281	end if
	282	end do
[190]	283
	284	if (.not. isvar("et"))then
	285	print(" ")
	286	print("'end_time_step' = "+ end_time_step +"h is invalid")
	287	print(" ")
[218]	288	print("Select another 'end_time_step'")
[190]	289	print(" ")
	290	exit
	291	end if
[162]	292
	293	delete(start_time_step)
	294	start_time_step=round(st,3)
	295	delete(end_time_step)
	296	end_time_step=round(et,3)
[154]	297
[175]	298	print(" ")
[566]	299	print("Output of time steps from "+t_all(start_time_step)/3600+" h = "+\
	300	t_all(start_time_step)+" s => index = "+start_time_step)
	301	print(" till "+t_all(end_time_step)/3600+" h = "+\
	302	t_all(end_time_step)+" s => index = "+end_time_step)
[175]	303	print(" ")
	304
[218]	305	t = t_all(start_time_step:end_time_step)/norm_t
[154]	306
	307	; ***************************************************
	308	; set up recourses
	309	; ***************************************************
	310
	311	res = True
	312	res@gsnDraw = False
	313	res@gsnFrame = False
	314	res@tmXBMode = True
	315	res@tmYLMode = True
	316	res@txFont = "helvetica"
	317	res@tiMainFont = "helvetica"
	318	res@tiXAxisFont = "helvetica"
	319	res@tiYAxisFont = "helvetica"
	320	res@tmXBLabelFont = "helvetica"
[161]	321	res@tmYLLabelFont = "helvetica"
	322	res@xyLineColors = (/237/)
[566]	323	res@trXMaxF = t_all(start_time_step)
	324	res@trXMinF = t_all(end_time_step)
[161]	325
[218]	326	res@lgLabelFontHeightF = 0.02
[161]	327
[154]	328	resP = True
	329	resP@txFont = "helvetica"
[218]	330	resP@txString = f_att@title+" time series "
[154]	331	resP@txFuncCode = "~"
	332	resP@txFontHeightF = 0.015
	333
[218]	334	res@tmXBMinorPerMajor = 4
	335	res@tmYLMinorPerMajor = 4
	336
[154]	337	res@vpWidthF=4
	338
[161]	339	txres = True
[154]	340
	341	; ***************************************************
	342	; read data and create plots
	343	; ***************************************************
	344
[532]	345	if ( format_out .EQ. "pdf" .OR. format_out .EQ. "ps" ) then
	346	format_out@wkPaperSize = "A4"
	347	end if
	348	if ( format_out .EQ. "png" ) then
	349	format_out@wkWidth = 1000
	350	format_out@wkHeight = 1000
	351	end if
	352
[154]	353	wks_ps = gsn_open_wks(format_out,file_out)
	354	gsn_define_colormap(wks_ps,"rainbow+white")
[161]	355	plot_ps=new(dim,graphic)
	356
[154]	357	n=0
[161]	358	minE=1.E27
	359	maxE=-1.E27
	360	minus=1.E27
	361	maxus=-1.E27
	362	minu=1.E27
	363	maxu=-1.E27
	364	minz=1.E27
	365	maxz=-1.E27
	366	minw=1.E27
	367	maxw=-1.E27
	368	minp=1.E27
	369	maxp=-1.E27
	370	mins=1.E27
	371	maxs=-1.E27
	372
	373	data = new((/dim,(end_time_step-start_time_step)+1/),float)
	374	unit = new(dim,string)
	375	data_0 = new((end_time_step-start_time_step)+1,float)
	376	data_0 = 0.0
	377	mini = new(dim,float)
	378	maxi = new(dim,float)
	379
	380	if (over .EQ. 1) then
	381	plot_E = gsn_csm_xy(wks_ps,t,data_0(:),res)
	382	plot_Es = gsn_csm_xy(wks_ps,t,data_0(:),res)
	383	plot_us = gsn_csm_xy(wks_ps,t,data_0(:),res)
	384	plot_ws = gsn_csm_xy(wks_ps,t,data_0(:),res)
	385	plot_umax = gsn_csm_xy(wks_ps,t,data_0(:),res)
	386	plot_vmax = gsn_csm_xy(wks_ps,t,data_0(:),res)
	387	plot_wmax = gsn_csm_xy(wks_ps,t,data_0(:),res)
	388	plot_z_i_wpt = gsn_csm_xy(wks_ps,t,data_0(:),res)
	389	plot_z_i_pt = gsn_csm_xy(wks_ps,t,data_0(:),res)
	390	plot_wpptp0 = gsn_csm_xy(wks_ps,t,data_0(:),res)
	391	plot_wpptp = gsn_csm_xy(wks_ps,t,data_0(:),res)
	392	plot_wpt = gsn_csm_xy(wks_ps,t,data_0(:),res)
	393	plot_pt_0_ = gsn_csm_xy(wks_ps,t,data_0(:),res)
	394	plot_pt_zp_ = gsn_csm_xy(wks_ps,t,data_0(:),res)
	395	end if
[162]	396
[161]	397
[310]	398	if (var .NE. "all") then
	399	vNam_temp = new((/dim/),string)
	400	comma = 0
	401	var_char = stringtocharacter(var)
	402	no_char = dimsizes(var_char)
	403
	404	do j = 0, no_char-1
	405	if(var_char(j) .eq. ",")
	406	comma = comma + 1
	407	end if
	408	end do
	409
[317]	410	if(comma .eq. 0 .or. comma .eq. 1)
	411	print(" ")
	412	print("The variables 'var="+var+"'" )
	413	print("do not exist on your input file;")
	414	print("be sure to have one comma before and after each variable")
	415	print(" ")
	416	exit
	417	end if
	418
[310]	419	indices = new((/comma/), integer)
	420	comma = 0
	421
	422	do j = 0, no_char-1
	423	if(var_char(j) .eq. ",")
	424	indices(comma) = j
	425	comma = comma + 1
	426	end if
	427	end do
	428
	429	do j = 0, comma -2
[566]	430	vNam_temp(j) = charactertostring(\
	431	var_char(indices(j)+1:indices(j+1)-1))
[310]	432	end do
	433
	434	delete(vNam)
	435	vNam = new((/comma-1/),string)
	436	do varn = 0, comma - 2
	437	vNam(varn) = vNam_temp(comma - 2 - varn)
	438	end do
	439
	440	do j = 0, comma-2
	441	count_check = 0
	442	do varn = 0, dim-1
	443	if(vNam(j) .ne. vNam_static(varn))
	444	count_check=count_check+1
	445	end if
[317]	446	end do
	447
[310]	448	if(count_check .eq. dim)
	449	print(" ")
	450	print("The variables 'var="+var+"'" )
	451	print("do not exist on your input file;")
	452	print("be sure to have one comma before and after each variable")
	453	print(" ")
	454	exit
	455	end if
	456	end do
	457
	458	dim = comma-1
	459	end if
	460
	461
	462	do varn = dim-1, 0, 1
	463
[161]	464	if( isStrSubset (vNam(varn), "time") )
	465	check = False
	466	else
	467	check = True
[310]	468	end if
[161]	469
	470	if(check) then
[162]	471
[218]	472	data_all = f[:]->$vNam(varn)$
	473	data_att = f_att->$vNam(varn)$
	474	unit(varn) = data_att@units
[310]	475
	476
[161]	477	data(varn,:)=data_all(start_time_step:end_time_step)
	478
	479	if (over .EQ. 1) then
	480
	481	mini(varn) = min(data(varn,:))
	482	maxi(varn) = max(data(varn,:))
	483
[566]	484	if (vNam(varn) .EQ. "E" .OR. vNam(varn) .EQ. "Es" .OR. \
	485	vNam(varn) .EQ. "E*") then
[161]	486	if (mini(varn) .EQ. maxi(varn)) then
	487	if (min(data(varn,:)) .EQ. 0)then
[218]	488	mini(varn)= mini(varn)-0.1
	489	maxi(varn)= maxi(varn)+0.1
[161]	490	end if
	491	if (min(data(varn,:)) .LT. 0)then
	492	mini(varn)= mini(varn)-1.+(mini(varn))/2
	493	maxi(varn)= maxi(varn)+1.-(maxi(varn))/2
	494	end if
	495	if (min(data(varn,:)) .GT. 0)then
	496	mini(varn)= mini(varn)-1.-(mini(varn))/2
	497	maxi(varn)= maxi(varn)+1.+(maxi(varn))/2
	498	end if
	499	end if
	500	minE=min((/minE,mini(varn)/))
	501	maxE=max((/maxE,maxi(varn)/))
	502	end if
	503
[526]	504	if (vNam(varn) .EQ. "us" .OR. vNam(varn) .EQ. "ws"\
	505	.OR. vNam(varn) .EQ. "u" .OR. vNam(varn) .EQ. "w") then
[161]	506	if (mini(varn) .EQ. maxi(varn)) then
	507	if (min(data(varn,:)) .EQ. 0)then
[218]	508	mini(varn)= mini(varn)-0.1
	509	maxi(varn)= maxi(varn)+0.1
[161]	510	end if
	511	if (min(data(varn,:)) .LT. 0)then
	512	mini(varn)= mini(varn)-1.+(mini(varn))/2
	513	maxi(varn)= maxi(varn)+1.-(maxi(varn))/2
	514	end if
	515	if (min(data(varn,:)) .GT. 0)then
	516	mini(varn)= mini(varn)-1.-(mini(varn))/2
	517	maxi(varn)= maxi(varn)+1.+(maxi(varn))/2
	518	end if
	519	end if
	520	minus=min((/minus,mini(varn)/))
	521	maxus=max((/maxus,maxi(varn)/))
	522	end if
	523
[566]	524	if (vNam(varn) .EQ. "umax" .OR. vNam(varn) .EQ. "vmax" .OR. \
	525	vNam(varn) .EQ. "wmax") then
[161]	526	if (mini(varn) .EQ. maxi(varn)) then
	527	if (mini(varn) .EQ. 0)then
[218]	528	mini(varn)= mini(varn)-0.1
	529	maxi(varn)= maxi(varn)+0.1
[161]	530	end if
	531	if (mini(varn) .LT. 0)then
	532	mini(varn)= mini(varn)-1.+(mini(varn))/2
	533	maxi(varn)= maxi(varn)+1.-(maxi(varn))/2
	534	end if
	535	if (mini(varn) .GT. 0)then
	536	mini(varn)= mini(varn)-1.-(mini(varn))/2
	537	maxi(varn)= maxi(varn)+1.+(maxi(varn))/2
	538	end if
	539	end if
	540	minu=min((/minu,mini(varn)/))
	541	maxu=max((/maxu,maxi(varn)/))
	542	end if
	543
	544	if (vNam(varn) .EQ. "z_i_wpt" .OR. vNam(varn) .EQ. "z_i_pt") then
	545	if (mini(varn) .EQ. maxi(varn)) then
	546	if (min(data(varn,:)) .EQ. 0)then
[218]	547	mini(varn)= mini(varn)-0.1
	548	maxi(varn)= maxi(varn)+0.1
[161]	549	end if
	550	if (min(data(varn,:)) .LT. 0)then
	551	mini(varn)= mini(varn)-1.+(mini(varn))/2
	552	maxi(varn)= maxi(varn)+1.-(maxi(varn))/2
	553	end if
	554	if (min(data(varn,:)) .GT. 0)then
	555	mini(varn)= mini(varn)-1.-(mini(varn))/2
	556	maxi(varn)= maxi(varn)+1.+(maxi(varn))/2
	557	end if
	558	end if
	559	minz=min((/minz,mini(varn)/))
	560	maxz=max((/maxz,maxi(varn)/))
	561	end if
	562
[566]	563	if (vNam(varn) .EQ. "wpptp0" .OR. vNam(varn) .EQ. "wpptp" .OR. \
	564	vNam(varn) .EQ. "wpt" .OR. vNam(varn) .EQ. \
	565	"w"+dq+"pt"+dq+"0" .OR. vNam(varn) .EQ. "w"+dq+"pt"+dq ) then
[161]	566	if (mini(varn) .EQ. maxi(varn)) then
	567	if (min(data(varn,:)) .EQ. 0)then
[218]	568	mini(varn)= mini(varn)-0.1
	569	maxi(varn)= maxi(varn)+0.1
[161]	570	end if
	571	if (min(data(varn,:)) .LT. 0)then
	572	mini(varn)= mini(varn)-1.+(mini(varn))/2
	573	maxi(varn)= maxi(varn)+1.-(maxi(varn))/2
	574	end if
	575	if (min(data(varn,:)) .GT. 0)then
	576	mini(varn)= mini(varn)-1.-(mini(varn))/2
	577	maxi(varn)= maxi(varn)+1.+(maxi(varn))/2
	578	end if
	579	end if
	580	minw=min((/minw,mini(varn)/))
	581	maxw=max((/maxw,maxi(varn)/))
	582	end if
	583
[526]	584	if (vNam(varn) .EQ. "pt_0_" .OR. vNam(varn) .EQ. "pt_zp_"\
	585	.OR. vNam(varn) .EQ. "pt(0)" .OR. vNam(varn) .EQ. "pt(zp)") then
[161]	586	if (mini(varn) .EQ. maxi(varn)) then
	587	if (min(data(varn,:)) .EQ. 0)then
[218]	588	mini(varn)= mini(varn)-0.1
	589	maxi(varn)= maxi(varn)+0.1
[161]	590	end if
	591	if (min(data(varn,:)) .LT. 0)then
	592	mini(varn)= mini(varn)-1.+(mini(varn))/2
	593	maxi(varn)= maxi(varn)+1.-(maxi(varn))/2
	594	end if
	595	if (min(data(varn,:)) .GT. 0)then
	596	mini(varn)= mini(varn)-1.-(mini(varn))/2
	597	maxi(varn)= maxi(varn)+1.+(maxi(varn))/2
	598	end if
	599	end if
	600	minp=min((/minp,mini(varn)/))
	601	maxp=max((/maxp,maxi(varn)/))
	602	end if
	603
	604	end if
	605	end if
	606	end do
[310]	607
	608	if (isStrSubset(data@long_name," SR " ))then
[162]	609	print(" ")
[566]	610	print("If you have outputs of statistic regions you cannot overlay "+\
	611	"variables;")
[310]	612	print("'over' is set to 0" )
[162]	613	print(" ")
[310]	614	over = 0
[175]	615	end if
[162]	616
[310]	617	do varn = dim-1, 0, 1
[161]	618
[154]	619	if( isStrSubset (vNam(varn), "time") )
	620	check = False
[161]	621	else
	622	check = True
	623	end if
[190]	624
	625	if (var .NE. "all") then
[161]	626	check = isStrSubset( var,","+vNam(varn)+"," )
	627	end if
[154]	628
	629	if(check) then
[310]	630
[161]	631	if (over .EQ. 1) then
[154]	632
[161]	633	res@gsnLeftString = "overlayed plot"
	634	res@gsnRightString = unit(varn)
[218]	635	if (norm_t .NE. 1.)then
[534]	636	res@tiXAxisString = "t ("+unit_t+")"
[218]	637	else
[534]	638	res@tiXAxisString = "t (s)"
[218]	639	end if
	640	res@tiYAxisString = " "
	641	res@tiXAxisFontHeightF = font_size
	642	res@txFontHeightF = font_size
	643	res@tiYAxisFontHeightF = font_size
[154]	644
[161]	645	if (vNam(varn) .EQ. "E")
	646	E=0
	647	res@xyLineColors = (/237/)
	648	res@xyLineLabelFontHeightF = 0.05
	649	res@xyLineLabelFontColor = 237
[310]	650	res@trYMaxF = maxE
	651	res@trYMinF = minE
[161]	652	plot_E = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	653	end if
[526]	654	if (vNam(varn) .EQ. "Es" .OR. vNam(varn) .EQ. "E*")
[161]	655	Es=0
	656	res@xyLineColors = (/144/)
	657	res@xyLineLabelFontHeightF = 0.05
	658	res@xyLineLabelFontColor = 144
	659	plot_Es = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	660	end if
	661
[526]	662	if (vNam(varn) .EQ. "us" .OR. vNam(varn) .EQ. "u*")
[161]	663	us=0
	664	res@xyLineColors = (/237/)
	665	res@xyLineLabelFontHeightF = 0.05
	666	res@xyLineLabelFontColor = 237
[310]	667	res@trYMaxF = maxus
	668	res@trYMinF = minus
[161]	669	plot_us = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	670	end if
[526]	671	if (vNam(varn) .EQ. "ws" .OR. vNam(varn) .EQ. "w*")
[161]	672	ws=0
	673	res@xyLineColors = (/144/)
	674	res@xyLineLabelFontHeightF = 0.05
	675	res@xyLineLabelFontColor = 144
	676	plot_ws = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	677	end if
	678
	679	if (vNam(varn) .EQ. "umax")
	680	u=0
	681	res@xyLineColors = (/237/)
	682	res@xyLineLabelFontHeightF = 0.05
	683	res@xyLineLabelFontColor = 237
[310]	684	res@trYMaxF = maxu
	685	res@trYMinF = minu
[161]	686	plot_umax = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	687	end if
	688	if (vNam(varn) .EQ. "vmax")
	689	v=0
	690	res@xyLineColors = (/144/)
	691	res@xyLineLabelFontHeightF = 0.05
	692	res@xyLineLabelFontColor = 144
	693	plot_vmax = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	694	end if
	695	if (vNam(varn) .EQ. "wmax")
	696	w=0
	697	res@xyLineColors = (/80/)
	698	res@xyLineLabelFontHeightF = 0.05
	699	res@xyLineLabelFontColor = 80
	700	plot_wmax = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	701	end if
[154]	702
[161]	703	if (vNam(varn) .EQ. "z_i_wpt")
	704	zw=0
	705	res@xyLineColors = (/237/)
	706	res@xyLineLabelFontHeightF = 0.05
	707	res@xyLineLabelFontColor = 237
[310]	708	res@trYMaxF = maxz
	709	res@trYMinF = minz
[161]	710	plot_z_i_wpt = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	711	end if
	712	if (vNam(varn) .EQ. "z_i_pt")
	713	z=0
	714	res@xyLineColors = (/144/)
	715	res@xyLineLabelFontHeightF = 0.05
	716	res@xyLineLabelFontColor = 144
	717	plot_z_i_pt = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	718	end if
	719
[566]	720	if (vNam(varn) .EQ. "wpptp0" .OR. \
	721	vNam(varn) .EQ. "w"+dq+"pt"+dq+"0" )
[161]	722	w0=0
	723	res@xyLineColors = (/237/)
	724	res@xyLineLabelFontHeightF = 0.05
	725	res@xyLineLabelFontColor = 237
[310]	726	res@trYMaxF = maxw
	727	res@trYMinF = minw
[161]	728	plot_wpptp0 = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	729	end if
[526]	730	if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "w"+dq+"pt"+dq)
[161]	731	wp=0
	732	res@xyLineColors = (/144/)
	733	res@xyLineLabelFontHeightF = 0.05
	734	res@xyLineLabelFontColor = 144
	735	plot_wpptp = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	736	end if
	737	if (vNam(varn) .EQ. "wpt")
	738	wt=0
	739	res@xyLineColors = (/80/)
	740	res@xyLineLabelFontHeightF = 0.05
	741	res@xyLineLabelFontColor = 80
	742	plot_wpt = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	743	end if
	744
[526]	745	if (vNam(varn) .EQ. "pt_0_" .OR. vNam(varn) .EQ. "pt(0)")
[161]	746	p=0
	747	res@xyLineColors = (/237/)
	748	res@xyLineLabelFontHeightF = 0.05
	749	res@xyLineLabelFontColor = 237
[310]	750	res@trYMaxF = maxp
	751	res@trYMinF = minp
[161]	752	plot_pt_0_ = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	753	end if
[526]	754	if (vNam(varn) .EQ. "pt_zp_" .OR. vNam(varn) .EQ. "pt(zp)")
[161]	755	pz=0
	756	res@xyLineColors = (/144/)
	757	res@xyLineLabelFontHeightF = 0.05
	758	res@xyLineLabelFontColor = 144
	759	plot_pt_zp_ = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	760	end if
	761
	762	end if
[154]	763	end if
[161]	764	end do
[310]	765
	766
	767	do varn = dim-1, 0, 1
[161]	768
	769	if( isStrSubset (vNam(varn), "time") )
	770	check = False
	771	else
	772	check = True
	773	end if
[190]	774
	775	if (var.NE. "all") then
[161]	776	check = isStrSubset( var,","+vNam(varn)+"," )
	777	end if
	778
	779	if(check) then
	780
	781	if (over .EQ. 1) then
	782
	783	if (vNam(varn) .EQ. "E" .AND. Es .NE. 1) then
	784	E=1
	785	overlay(plot_E,plot_Es)
	786	n=n+1
[162]	787	plot_ps(n) = plot_E
	788
	789	; ***************************************************
[566]	790	; legend for overlaid plot
[162]	791	; ***************************************************
	792
	793	lgres = True
	794	lgMonoDashIndex = False
	795	lgres@lgLabelFont = "helvetica"
	796	lgres@lgLabelFontHeightF = .1
	797	lgres@vpWidthF = 0.4
	798	lgres@vpHeightF = 0.4
	799	lgres@lgDashIndexes = (/0,0,0/)
	800	lgres@lgLineColors = (/237,144,80/)
[526]	801	lbid = gsn_create_legend(wks_ps,2,(/"E","E*"/),lgres)
[162]	802
	803	amres = True
	804	amres@amParallelPosF = 0.6
	805	amres@amOrthogonalPosF = -0.2
	806	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	807	end if
[566]	808	if ((vNam(varn) .EQ. "Es" .OR. vNam(varn) .EQ. "E*") .AND. \
	809	E .NE. 1) then
[161]	810	Es=1
	811	overlay(plot_E,plot_Es)
	812	n=n+1
[162]	813	plot_ps(n) = plot_E
	814
	815	; ***************************************************
[566]	816	; legend for overlaid plot
[162]	817	; ***************************************************
	818
	819	lgres = True
	820	lgMonoDashIndex = False
	821	lgres@lgLabelFont = "helvetica"
	822	lgres@lgLabelFontHeightF = .1
	823	lgres@vpWidthF = 0.4
	824	lgres@vpHeightF = 0.4
	825	lgres@lgDashIndexes = (/0,0,0/)
	826	lgres@lgLineColors = (/237,144,80/)
[526]	827	lbid = gsn_create_legend(wks_ps,2,(/"E","E*"/),lgres)
[162]	828
	829	amres = True
	830	amres@amParallelPosF = 0.6
	831	amres@amOrthogonalPosF = -0.2
	832	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	833	end if
	834
[566]	835	if ((vNam(varn) .EQ. "us" .OR. vNam(varn) .EQ. "u*") .AND. \
	836	ws .NE. 1) then
[161]	837	us=1
	838	overlay(plot_us,plot_ws)
	839	n=n+1
	840	plot_ps(n) = plot_us
[162]	841
	842	; ***************************************************
[566]	843	; legend for overlaid plot
[162]	844	; ***************************************************
	845
	846	lgres = True
	847	lgMonoDashIndex = False
	848	lgres@lgLabelFont = "helvetica"
	849	lgres@lgLabelFontHeightF = .1
	850	lgres@vpWidthF = 0.4
	851	lgres@vpHeightF = 0.4
	852	lgres@lgDashIndexes = (/0,0,0/)
	853	lgres@lgLineColors = (/237,144,80/)
[526]	854	lbid = gsn_create_legend(wks_ps,2,(/"u","w"/),lgres)
[162]	855
	856	amres = True
	857	amres@amParallelPosF = 0.6
	858	amres@amOrthogonalPosF = -0.2
	859	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	860	end if
[566]	861	if ((vNam(varn) .EQ. "ws" .OR. vNam(varn) .EQ. "w*") .AND. \
	862	us .NE. 1) then
[161]	863	ws=1
	864	overlay(plot_us,plot_ws)
	865	n=n+1
	866	plot_ps(n) = plot_us
[162]	867
	868	; ***************************************************
[566]	869	; legend for overlaid plot
[162]	870	; ***************************************************
	871
	872	lgres = True
	873	lgMonoDashIndex = False
	874	lgres@lgLabelFont = "helvetica"
	875	lgres@lgLabelFontHeightF = .1
	876	lgres@vpWidthF = 0.4
	877	lgres@vpHeightF = 0.4
	878	lgres@lgDashIndexes = (/0,0,0/)
	879	lgres@lgLineColors = (/237,144,80/)
[526]	880	lbid = gsn_create_legend(wks_ps,2,(/"u","w"/),lgres)
[162]	881
	882	amres = True
	883	amres@amParallelPosF = 0.6
	884	amres@amOrthogonalPosF = -0.2
	885	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	886	end if
	887
	888	if (vNam(varn) .EQ. "umax" .AND. v .NE. 1)
	889	if (w .NE. 1) then
	890	u=1
	891	overlay(plot_umax,plot_vmax)
	892	overlay(plot_umax,plot_wmax)
[310]	893
[161]	894	n=n+1
	895	plot_ps(n) = plot_umax
[162]	896
	897	; ***************************************************
[566]	898	; legend for overlaid plot
[162]	899	; ***************************************************
	900
	901	lgres = True
	902	lgMonoDashIndex = False
	903	lgres@lgLabelFont = "helvetica"
	904	lgres@lgLabelFontHeightF = .1
	905	lgres@vpWidthF = 0.4
	906	lgres@vpHeightF = 0.4
	907	lgres@lgDashIndexes = (/0,0,0/)
	908	lgres@lgLineColors = (/237,144,80/)
[566]	909	lbid = gsn_create_legend(\
	910	wks_ps,3,(/"umax","vmax","wmax"/),lgres)
[162]	911
	912	amres = True
	913	amres@amParallelPosF = 0.6
	914	amres@amOrthogonalPosF = -0.2
	915	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	916	end if
	917	end if
	918	if (vNam(varn) .EQ. "vmax" .AND. u .NE. 1)
	919	if (w .NE. 1) then
	920	v=1
	921	overlay(plot_umax,plot_vmax)
	922	overlay(plot_umax,plot_wmax)
	923	n=n+1
	924	plot_ps(n) = plot_umax
[162]	925
	926	; ***************************************************
[566]	927	; legend for overlaid plot
[162]	928	; ***************************************************
	929
	930	lgres = True
	931	lgMonoDashIndex = False
	932	lgres@lgLabelFont = "helvetica"
	933	lgres@lgLabelFontHeightF = .1
	934	lgres@vpWidthF = 0.4
	935	lgres@vpHeightF = 0.4
	936	lgres@lgDashIndexes = (/0,0,0/)
	937	lgres@lgLineColors = (/237,144,80/)
[566]	938	lbid = gsn_create_legend(\
	939	wks_ps,3,(/"umax","vmax","wmax"/),lgres)
[162]	940	amres = True
	941	amres@amParallelPosF = 0.6
	942	amres@amOrthogonalPosF = -0.2
	943	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	944	end if
	945	end if
	946	if (vNam(varn) .EQ. "wmax" .AND. v .NE. 1)
	947	if(u .NE. 1) then
	948	w=1
	949	overlay(plot_umax,plot_vmax)
	950	overlay(plot_umax,plot_wmax)
	951	n=n+1
	952	plot_ps(n) = plot_umax
[162]	953
	954	; ***************************************************
[566]	955	; legend for overlaid plot
[162]	956	; ***************************************************
	957
	958	lgres = True
	959	lgMonoDashIndex = False
	960	lgres@lgLabelFont = "helvetica"
	961	lgres@lgLabelFontHeightF = .1
	962	lgres@vpWidthF = 0.4
	963	lgres@vpHeightF = 0.4
	964	lgres@lgDashIndexes = (/0,0,0/)
	965	lgres@lgLineColors = (/237,144,80/)
[566]	966	lbid = gsn_create_legend(\
	967	wks_ps,3,(/"umax","vmax","wmax"/),lgres)
[162]	968	amres = True
	969	amres@amParallelPosF = 0.6
	970	amres@amOrthogonalPosF = -0.2
	971	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	972	end if
	973	end if
	974
	975	if (vNam(varn) .EQ. "z_i_wpt" .AND. z .NE. 1) then
	976	zw=1
	977	overlay(plot_z_i_wpt,plot_z_i_pt)
	978	n=n+1
[162]	979	plot_ps(n) = plot_z_i_wpt
	980
	981	; ***************************************************
[566]	982	; legend for overlaid plot
[162]	983	; ***************************************************
	984
	985	lgres = True
	986	lgMonoDashIndex = False
	987	lgres@lgLabelFont = "helvetica"
	988	lgres@lgLabelFontHeightF = .1
	989	lgres@vpWidthF = 0.4
	990	lgres@vpHeightF = 0.4
	991	lgres@lgDashIndexes = (/0,0,0/)
	992	lgres@lgLineColors = (/237,144,80/)
[566]	993	lbid = gsn_create_legend(wks_ps,2,(/"z_i_wpt","z_i_pt"/),lgres)
[162]	994
	995	amres = True
	996	amres@amParallelPosF = 0.6
	997	amres@amOrthogonalPosF = -0.2
	998	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	999	end if
	1000	if (vNam(varn) .EQ. "z_i_pt" .AND. zw .NE. 1) then
	1001	z=1
	1002	overlay(plot_z_i_wpt,plot_z_i_pt)
	1003	n=n+1
[162]	1004	plot_ps(n) = plot_z_i_wpt
	1005
	1006	; ***************************************************
[566]	1007	; legend for overlaid plot
[162]	1008	; ***************************************************
	1009
	1010	lgres = True
	1011	lgMonoDashIndex = False
	1012	lgres@lgLabelFont = "helvetica"
	1013	lgres@lgLabelFontHeightF = .1
	1014	lgres@vpWidthF = 0.4
	1015	lgres@vpHeightF = 0.4
	1016	lgres@lgDashIndexes = (/0,0,0/)
	1017	lgres@lgLineColors = (/237,144,80/)
[566]	1018	lbid = gsn_create_legend(wks_ps,2,(/"z_i_wpt","z_i_pt"/),lgres)
[162]	1019
	1020	amres = True
	1021	amres@amParallelPosF = 0.6
	1022	amres@amOrthogonalPosF = -0.2
	1023	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1024	end if
	1025
[526]	1026	if ((vNam(varn) .EQ. "wpptp0" .OR. vNam(varn) .EQ. "w"+dq+"pt"+dq+"0") .AND. wp .NE. 1)
[161]	1027	if (wt .NE. 1) then
	1028	w0=1
	1029	overlay(plot_wpptp0,plot_wpptp)
	1030	overlay(plot_wpptp0,plot_wpt)
	1031	n=n+1
	1032	plot_ps(n) = plot_wpptp0
[162]	1033
	1034	; ***************************************************
[566]	1035	; legend for overlaid plot
[162]	1036	; ***************************************************
	1037
	1038	lgres = True
	1039	lgMonoDashIndex = False
	1040	lgres@lgLabelFont = "helvetica"
	1041	lgres@lgLabelFontHeightF = .1
	1042	lgres@vpWidthF = 0.4
	1043	lgres@vpHeightF = 0.4
	1044	lgres@lgDashIndexes = (/0,0,0/)
	1045	lgres@lgLineColors = (/237,144,80/)
[566]	1046	lbid = gsn_create_legend(\
	1047	wks_ps,3,(/"wpptp0","wpptp","wpt"/),lgres)
[162]	1048	amres = True
	1049	amres@amParallelPosF = 0.6
	1050	amres@amOrthogonalPosF = -0.2
	1051	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1052	end if
	1053	end if
[566]	1054	if ((vNam(varn) .EQ. "wpptp" .OR. \
	1055	vNam(varn) .EQ. "w"+dq+"pt"+dq) .AND. w0 .NE. 1)
[161]	1056	if (wt .NE. 1) then
	1057	wp=1
	1058	overlay(plot_wpptp0,plot_wpptp)
	1059	overlay(plot_wpptp0,plot_wpt)
	1060	n=n+1
[162]	1061	plot_ps(n) = plot_wpptp0
	1062
	1063	; ***************************************************
[566]	1064	; legend for overlaid plot
[162]	1065	; ***************************************************
	1066
	1067	lgres = True
	1068	lgMonoDashIndex = False
	1069	lgres@lgLabelFont = "helvetica"
	1070	lgres@lgLabelFontHeightF = .1
	1071	lgres@vpWidthF = 0.4
	1072	lgres@vpHeightF = 0.4
	1073	lgres@lgDashIndexes = (/0,0,0/)
	1074	lgres@lgLineColors = (/237,144,80/)
[566]	1075	lbid = gsn_create_legend(\
	1076	wks_ps,3,(/"wpptp0","wpptp","wpt"/),lgres)
[162]	1077
	1078	amres = True
	1079	amres@amParallelPosF = 0.6
	1080	amres@amOrthogonalPosF = -0.2
	1081	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1082	end if
	1083	end if
	1084	if (vNam(varn) .EQ. "wpt" .AND. wp .NE. 1)
	1085	if (w0 .NE. 1) then
	1086	wt=1
	1087	overlay(plot_wpptp0,plot_wpptp)
	1088	overlay(plot_wpptp0,plot_wpt)
	1089	n=n+1
[162]	1090	plot_ps(n) = plot_wpptp0
	1091
	1092	; ***************************************************
[566]	1093	; legend for overlaid plot
[162]	1094	; ***************************************************
	1095
	1096	lgres = True
	1097	lgMonoDashIndex = False
	1098	lgres@lgLabelFont = "helvetica"
	1099	lgres@lgLabelFontHeightF = .1
	1100	lgres@vpWidthF = 0.4
	1101	lgres@vpHeightF = 0.4
	1102	lgres@lgDashIndexes = (/0,0,0/)
	1103	lgres@lgLineColors = (/237,144,80/)
[566]	1104	lbid = gsn_create_legend(\
	1105	wks_ps,3,(/"wpptp0","wpptp","wpt"/),lgres)
[162]	1106	amres = True
	1107	amres@amParallelPosF = 0.6
	1108	amres@amOrthogonalPosF = -0.2
	1109	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1110	end if
	1111	end if
	1112
[566]	1113	if ((vNam(varn) .EQ. "pt_0_" .OR. vNam(varn) .EQ. "pt(0)") .AND. \
	1114	pz .NE. 1) then
[161]	1115	p=1
	1116	overlay(plot_pt_0_,plot_pt_zp_)
	1117	n=n+1
[162]	1118	plot_ps(n) = plot_pt_0_
	1119
	1120	; ***************************************************
[566]	1121	; legend for overlaid plot
[162]	1122	; ***************************************************
	1123
	1124	lgres = True
	1125	lgMonoDashIndex = False
	1126	lgres@lgLabelFont = "helvetica"
	1127	lgres@lgLabelFontHeightF = .1
	1128	lgres@vpWidthF = 0.4
	1129	lgres@vpHeightF = 0.4
	1130	lgres@lgDashIndexes = (/0,0,0/)
	1131	lgres@lgLineColors = (/237,144,80/)
[566]	1132	lbid = gsn_create_legend(wks_ps,2,(/"pt_0_","pt_zp_"/),lgres)
[162]	1133
	1134	amres = True
	1135	amres@amParallelPosF = 0.6
	1136	amres@amOrthogonalPosF = -0.2
	1137	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1138	end if
[566]	1139	if ((vNam(varn) .EQ. "pt_zp_" .OR. vNam(varn) .EQ. "pt(zp)") .AND.\
	1140	p .NE. 1) then
[161]	1141	pz=1
	1142	overlay(plot_pt_0_,plot_pt_zp_)
	1143	n=n+1
[162]	1144	plot_ps(n) = plot_pt_0_
	1145
	1146	; ***************************************************
[566]	1147	; legend for overlaid plot
[162]	1148	; ***************************************************
	1149
	1150	lgres = True
	1151	lgMonoDashIndex = False
	1152	lgres@lgLabelFont = "helvetica"
	1153	lgres@lgLabelFontHeightF = .1
	1154	lgres@vpWidthF = 0.4
	1155	lgres@vpHeightF = 0.4
	1156	lgres@lgDashIndexes = (/0,0,0/)
	1157	lgres@lgLineColors = (/237,144,80/)
[566]	1158	lbid = gsn_create_legend(wks_ps,2,(/"pt_0_","pt_zp_"/),lgres)
[162]	1159
	1160	amres = True
	1161	amres@amParallelPosF = 0.6
	1162	amres@amOrthogonalPosF = -0.2
	1163	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1164	end if
[162]	1165
[526]	1166	if(vNam(varn) .NE. "pt_zp_" .AND. vNam(varn) .NE. "pt(zp)" .AND. \
[566]	1167	vNam(varn) .NE. "pt_0_" .AND. vNam(varn) .NE. "pt(0)" .AND. \
	1168	vNam(varn) .NE. "wpt" .AND. vNam(varn) .NE. "wpptp" .AND. \
	1169	vNam(varn) .NE. "w"+dq+"pt"+dq .AND. \
	1170	vNam(varn) .NE. "wpptp0" .AND. \
	1171	vNam(varn) .NE. "w"+dq+"pt"+dq+"0" .AND. \
	1172	vNam(varn) .NE. "z_i_pt" .AND. \
	1173	vNam(varn) .NE. "z_i_wpt" .AND. vNam(varn) .NE. "wmax" .AND. \
	1174	vNam(varn) .NE. "vmax" .AND. vNam(varn) .NE. "umax" .AND. \
	1175	vNam(varn) .NE. "ws" .AND. vNam(varn) .NE. "w*" .AND. \
	1176	vNam(varn) .NE. "us" .AND. vNam(varn) .NE. "u*" .AND. \
	1177	vNam(varn) .NE. "Es" .AND. vNam(varn) .NE. "E*" .AND. \
[526]	1178	vNam(varn) .NE. "E") then
[162]	1179
[161]	1180	n=n+1
	1181	res@xyLineColors = (/237/)
	1182	res@xyLabelMode = False
	1183	res@gsnLeftString = vNam(varn)
	1184	res@gsnRightString = unit(varn)
[310]	1185	res@trYMaxF = max(data(varn,:))
	1186	res@trYMinF = min(data(varn,:))
[161]	1187	if (min(data(varn,:)) .EQ. max(data(varn,:))) then
	1188	if (min(data(varn,:)) .EQ. 0)then
[310]	1189	res@trYMaxF = max(data(varn,:))-0.1
	1190	res@trYMinF = min(data(varn,:))+0.1
[161]	1191	end if
	1192	if (min(data(varn,:)) .LT. 0)then
[310]	1193	res@trYMaxF = max(data(varn,:))+(min(data(varn,:)))/2
	1194	res@trYMinF = min(data(varn,:))-(max(data(varn,:)))/2
[161]	1195	end if
	1196	if (min(data(varn,:)) .GT. 0)then
[310]	1197	res@trYMaxF = max(data(varn,:))-(min(data(varn,:)))/2
	1198	res@trYMinF = min(data(varn,:))+(max(data(varn,:)))/2
[161]	1199	end if
	1200	end if
	1201	plot_ps(n) = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	1202	end if
[310]	1203
	1204	; ***************************************************
	1205	; transparent legend as placeholder for single plots
	1206	; ***************************************************
	1207	lgresmono = True
	1208	lgMonoDashIndex = False
	1209	lgresmono@lgMonoLineColor = True
	1210	lgresmono@lgLabelFont = "helvetica"
	1211	lgresmono@lgLabelFontHeightF = .1
	1212	lgresmono@vpWidthF = 0.4
	1213	lgresmono@vpHeightF = 0.4
	1214	lgresmono@lgLineColor = 0
	1215	lgresmono@lgPerimOn = False
	1216	lbid = gsn_create_legend(wks_ps,1,(/""/),lgresmono)
	1217
	1218	amres = True
	1219	amres@amParallelPosF = 0.6
	1220	amres@amOrthogonalPosF = -0.2
	1221	annoid1 = gsn_add_annotation(plot_ps(n),lbid,amres)
[161]	1222
	1223	else
	1224
	1225	print("plot of " + vNam(varn))
	1226
	1227	n=n+1
	1228	res@xyLineColors = (/237/)
	1229	res@gsnLeftString = vNam(varn)
	1230	res@gsnRightString = unit(varn)
[310]	1231
[218]	1232	if (norm_t .NE. 1.)then
[534]	1233	res@tiXAxisString = "t ("+unit_t+")"
[218]	1234	else
[534]	1235	res@tiXAxisString = "t (s)"
[218]	1236	end if
	1237	res@tiYAxisString = " "
	1238	res@tiXAxisFontHeightF = font_size
	1239	res@txFontHeightF = font_size
	1240	res@tiYAxisFontHeightF = font_size
	1241	res@tmXBLabelFontHeightF = font_size
	1242	res@tmYLLabelFontHeightF = font_size
[310]	1243	res@trYMaxF = max(data(varn,:))
	1244	res@trYMinF = min(data(varn,:))
[161]	1245	if (min(data(varn,:)) .EQ. max(data(varn,:))) then
	1246	if (min(data(varn,:)) .EQ. 0)then
[310]	1247	res@trYMaxF = max(data(varn,:))-0.1
	1248	res@trYMinF = min(data(varn,:))+0.1
[161]	1249	end if
	1250	if (min(data(varn,:)) .LT. 0)then
[310]	1251	res@trYMaxF = max(data(varn,:))+(min(data(varn,:)))/2
	1252	res@trYMinF = min(data(varn,:))-(max(data(varn,:)))/2
[161]	1253	end if
	1254	if (min(data(varn,:)) .GT. 0)then
[310]	1255	res@trYMaxF = max(data(varn,:))-(min(data(varn,:)))/2
	1256	res@trYMinF = min(data(varn,:))+(max(data(varn,:)))/2
[161]	1257	end if
	1258	end if
	1259	plot_ps(n) = gsn_csm_xy(wks_ps,t,data(varn,:),res)
	1260
	1261	end if
	1262	end if
[154]	1263	end do
[161]	1264
[154]	1265	; ***************************************************
	1266	; merge plots onto one page
	1267	; ***************************************************
[532]	1268
	1269	no_frames = 0
[162]	1270
[566]	1271	if ((format_out .EQ. "eps" .OR. format_out .EQ. "epsi") .AND. \
	1272	n .gt. no_rows*no_columns) then
[161]	1273	gsn_panel(wks_ps,plot_ps(1:n),(/n,1/),resP)
[250]	1274	print(" ")
	1275	print("Outputs to .eps or .epsi have only one frame")
	1276	print(" ")
[161]	1277	else
[218]	1278	do np = 1,n,no_rows*no_columns
	1279	if ( np + no_rows*no_columns .gt. n) then
	1280	gsn_panel(wks_ps, plot_ps(np:n),(/no_rows,no_columns/),resP)
[532]	1281	no_frames = no_frames + 1
[161]	1282	else
[566]	1283	gsn_panel(wks_ps, plot_ps(np:np+no_rows*no_columns-1),\
	1284	(/no_rows,no_columns/),resP)
[532]	1285	no_frames = no_frames + 1
[161]	1286	end if
	1287	end do
	1288	end if
[154]	1289
[532]	1290	if (format_out .EQ. "png" ) then
	1291	png_output = new((/no_frames/), string)
	1292	j = 0
	1293	do i=0, no_frames-1
	1294	j = i + 1
	1295	if (j .LE. 9) then
	1296	png_output(i) = file_out+".00000"+j+".png"
	1297	end if
	1298	if (j .GT. 9 .AND. j .LE. 99) then
	1299	png_output(i) = file_out+".0000"+j+".png"
	1300	end if
	1301	if (j .GT. 99 .AND. j .LE. 999) then
	1302	png_output(i) = file_out+".000"+j+".png"
	1303	end if
	1304	if (j .GT. 999) then
	1305	png_output(i) = file_out+".00"+j+".png"
	1306	end if
	1307
	1308	;using imagemagick's convert for reducing the white
	1309	;space around the plot
	1310	cmd = "convert -geometry 1000x1000 -density 300 -trim " + \
	1311	png_output(i) + " " + png_output(i)
	1312	system(cmd)
	1313	end do
	1314
	1315	print(" ")
	1316	print("Output to: "+ png_output)
	1317	print(" ")
	1318	else
	1319	print(" ")
	1320	print("Output to: " + file_out +"."+ format_out)
	1321	print(" ")
	1322	end if
[154]	1323
[532]	1324
[154]	1325	end

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

Download in other formats:

Original Format