Home

Context Navigation

timeseries.ncl @ 188

Last change on this file since 188 was 175, checked in by steinfeld, 17 years ago
NCL scripts for spectra added. Bug fix concerning spectra in netcdf.f90
File size: 51.0 KB

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

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

Download in other formats:

Original Format