Home

Context Navigation

timeseries.ncl @ 482

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

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

Download in other formats:

Original Format