Home

Context Navigation

timeseries.ncl @ 624

Last change on this file since 624 was 585, checked in by heinze, 14 years ago
Bugfix: enable plot of data if it is of kind double instead of kind float
File size: 50.6 KB

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

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

Download in other formats:

Original Format