Home

Context Navigation

timeseries.ncl @ 424

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

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

Download in other formats:

Original Format