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

Last change on this file since 3856 was 2837, checked in by gronemeier, 7 years ago

palmplot: config file of ncl scripts can be chosen in shell command

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