source: palm/trunk/SCRIPTS/NCL/timeseries.ncl @ 405

Last change on this file since 405 was 317, checked in by heinze, 15 years ago

Bugfix in timeseries.ncl

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