Home

Context Navigation

timeseries.ncl @ 532

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

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

Download in other formats:

Original Format