Home

Context Navigation

timeseries.ncl @ 169

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

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

Download in other formats:

Original Format