source: palm/trunk/SCRIPTS/NCL/profiles.ncl @ 162

Last change on this file since 162 was 162, checked in by letzel, 14 years ago
  • NCL scripts in trunk/SCRIPTS/NCL updated
File size: 137.4 KB
Line 
1load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
2load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
3load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
4 
5begin
6
7   ; ***************************************************
8   ; read parameter_list
9   ; ***************************************************
10
11   if (isfilepresent("~/.ncl_preferences")) then
12      parameter = asciiread("~/.ncl_preferences",75,"string")
13      delete(parameter@_FillValue)
14   else
15      print(" ")
16      print("Please copy '.ncl_preferences' into your $home dircetory")
17      print(" ")
18      exit
19   end if
20 
21   ; ***************************************************
22   ; set up default parameter values and strings if not assigned in prompt or parameter list
23   ; ***************************************************
24 
25   if ( .not. isvar("file_in") ) then                   ; path+name of input file     
26      if (parameter(7) .EQ. "input file") then
27         print(" ")
28         print("Please provide input file 'file_in = ' either in prompt or parameter_list")
29         print(" ")
30         exit
31      else
32         file_in = parameter(7)
33      end if     
34   end if
35   if ( .not. isvar("format_out") ) then                ; format of output file
36      format_out = "x11"
37      if (parameter(9) .NE. "x11") then
38         format_out = parameter(9) 
39      end if
40   end if
41   if ( .not. isvar("file_out") ) then                  ; path+name of output file
42      file_out = "test"
43      if (parameter(11) .NE. "test_ts") then
44         file_out = parameter(11) 
45     end if     
46   end if
47   if ( .not. isvar("no_columns") ) then                ; number of plots in one row
48      no_columns = 1
49      if (parameter(17) .NE. "1") then
50         no_columns = stringtointeger(parameter(17)) 
51      end if
52   end if
53   if ( .not. isvar("no_lines") ) then                  ; number of plot-lines on one sheet
54      no_lines = 2
55      if (parameter(19) .NE. "2") then
56         no_lines = stringtointeger(parameter(19)) 
57      end if
58   end if
59   if ( .not. isvar("combine") ) then                   ; color of lines
60      combine = 0
61      if (parameter(23) .NE. "0") then
62         combine = stringtointeger(parameter(23))
63         if (stringtointeger(parameter(23)) .NE. 1) then
64            print(" ")
65            print("Please set 'combine' to 0 or 1")
66            print(" ")
67            exit
68         end if
69      end if
70   end if
71   if (combine .EQ. 1) then
72      if( .not. isvar("c_var") ) then
73         if (parameter(27) .EQ. "c_variables") then
74            print(" ")
75            print("Please select variables for overlaying ('c_var')")
76            print(" ")
77            exit
78         else
79            c_var=parameter(27)
80         end if
81      end if
82   end if
83   if ( .not. isvar("black") ) then                     ; color of lines
84      black = 0
85      if (parameter(31) .NE. "0") then
86         black = stringtointeger(parameter(31))
87         if (stringtointeger(parameter(31)) .NE. 1) then
88            print(" ")
89            print("Please set 'black' to 0 or 1")
90            print(" ")
91            exit
92         end if
93      end if
94   end if
95   if ( .not. isvar("dash") ) then                      ; pattern of lines
96      dash = 0
97      if (parameter(29) .NE. "0") then
98         dash = stringtointeger(parameter(29))
99         if (stringtointeger(parameter(29)) .NE. 1) then
100            print(" ")
101            print("Please set 'dash' to 0 or 1")
102            print(" ")
103            exit
104         end if 
105      end if
106   end if
107   if ( .not. isvar("over") ) then                      ; switches overlaying plots on
108      over = 0
109      if (stringtointeger(parameter(37)) .NE. 0) then
110         over = stringtointeger(parameter(37))
111         if (stringtointeger(parameter(37)) .NE. 1) then
112            print(" ")
113            print("Please set 'over' to 0 or 1")
114            print(" ")
115            exit
116         end if   
117      end if
118   end if
119 
120   ; ***************************************************
121   ; open input file
122   ; ***************************************************
123
124   f=addfile( file_in,"r")
125   
126   vNam = getfilevarnames(f)
127   print(" ")
128   print("Variable on netCDF file: " + vNam)
129   print(" ")
130   dim = dimsizes(vNam)
131   if (dim .EQ. 0) then
132      print(" ")
133      print("There are no data on file")
134      print(" ")
135   end if
136
137   do varn = dim-1,0,1
138      if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
139         varn=varn-1
140         continue
141      end if
142      if (vNam(varn) .EQ. "u")
143         z_u = f->zu
144         break
145      else
146         if (vNam(varn) .EQ. "v")then
147            z_u = f->zv
148            break
149         else
150            if(vNam(varn) .EQ. "pt")then
151               z_u = f->zpt
152               break
153            else
154               if(vNam(varn) .EQ. "vpt")then
155                  z_u = f->zvpt
156                  break
157               else
158                  if(vNam(varn) .EQ. "lpt")then
159                     z_u = f->zlpt
160                     break
161                  else   
162                     if(vNam(varn) .EQ. "q")then
163                        z_u = f->zq
164                        break
165                     else
166                        if(vNam(varn) .EQ. "qv")then
167                           z_u = f->zqv
168                           break
169                        else
170                           if(vNam(varn) .EQ. "ql")then
171                              z_u = f->zql
172                              break
173                           else
174                              if(vNam(varn) .EQ. "rho")then
175                                 z_u = f->zrho
176                                 break
177                              else
178                                 if(vNam(varn) .EQ. "s")then
179                                    z_u = f->zs
180                                    break
181                                 else
182                                    if(vNam(varn) .EQ. "sa")then
183                                       z_u = f->zsa
184                                       break
185                                    else
186                                       if(vNam(varn) .EQ. "e")then
187                                          z_u = f->ze
188                                          break
189                                       else
190                                          if(vNam(varn) .EQ. "es")then
191                                             z_u = f->zes
192                                             break
193                                          else
194                                             if(vNam(varn) .EQ. "km")then
195                                                z_u = f->zkm
196                                                break
197                                             else
198                                                if(vNam(varn) .EQ. "kh")then
199                                                   z_u = f->zkh
200                                                   break
201                                                else
202                                                   if(vNam(varn) .EQ. "l")then
203                                                      z_u = f->zl
204                                                      break
205                                                   else
206                                                      if(vNam(varn) .EQ. "us2")then
207                                                         z_u = f->zus2
208                                                         break
209                                                      else
210                                                         if(vNam(varn) .EQ. "vs2")then
211                                                            z_u = f->zvs2
212                                                            break
213                                                         else
214                                                            if(vNam(varn) .EQ. "pts2")then
215                                                               z_u = f->zpts2
216                                                               break
217                                                            else
218                                                               if(vNam(varn) .EQ. "wsususodz")then
219                                                                  z_u = f->zwsususodz
220                                                                  break
221                                                               else
222                                                                  if(vNam(varn) .EQ. "wspsodz")then
223                                                                     z_u = f->zwspsodz
224                                                                     break
225                                                                  else
226                                                                     if(vNam(varn) .EQ. "wpeodz")then
227                                                                        z_u = f->zwpeodz
228                                                                        break                                                                       
229                                                                     end if
230                                                                  end if
231                                                               end if
232                                                            end if
233                                                         end if
234                                                      end if
235                                                   end if
236                                                end if
237                                             end if
238                                          end if
239                                       end if
240                                    end if
241                                 end if
242                              end if
243                           end if
244                        end if
245                     end if
246                  end if
247               end if
248            end if
249         end if
250      end if
251      varn=varn-1
252   end do
253
254   do varn=dim-1,0,1
255      if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
256         varn=varn-1
257         continue
258      end if
259      if (vNam(varn) .EQ. "w")
260         z_w = f->zw
261         break
262      else
263         if (vNam(varn) .EQ. "wpup")then
264            z_w = f->zwpup
265            break
266         else
267            if(vNam(varn) .EQ. "wsus")then
268               z_w = f->zwsus
269               break
270            else
271               if(vNam(varn) .EQ. "wu")then
272                  z_w = f->zwu
273                  break
274               else
275                  if(vNam(varn) .EQ. "wpvp")then
276                     z_w = f->zwpvp
277                     break
278                  else   
279                     if(vNam(varn) .EQ. "wsvs")then
280                        z_w = f->zwsvs
281                        break
282                     else
283                        if(vNam(varn) .EQ. "wv")then
284                           z_w = f->zwv
285                           break
286                        else
287                           if(vNam(varn) .EQ. "wptpp")then
288                              z_w = f->zwptpp
289                              break
290                           else
291                              if(vNam(varn) .EQ. "wspts")then
292                                 z_w = f->zwspts
293                                 break
294                              else
295                                 if(vNam(varn) .EQ. "wpt")then
296                                    z_w = f->zwpt
297                                    break
298                                 else
299                                    if(vNam(varn) .EQ. "wsptsBC")then
300                                       z_w = f->zwsptsBC
301                                       break
302                                    else
303                                       if(vNam(varn) .EQ. "wptBC")then
304                                          z_w = f->zwptBC
305                                          break
306                                       else
307                                          if(vNam(varn) .EQ. "wpvptp")then
308                                             z_w = f->zwpvptp
309                                             break
310                                          else
311                                             if(vNam(varn) .EQ. "wsvpts")then
312                                                z_w = f->zwsvpts
313                                                break
314                                             else
315                                                if(vNam(varn) .EQ. "wvpt")then
316                                                   z_w = f->zwvpt
317                                                   break
318                                                else
319                                                   if(vNam(varn) .EQ. "wpqp")then
320                                                      z_w = f->zwpqp
321                                                      break
322                                                   else
323                                                      if(vNam(varn) .EQ. "wsqs")then
324                                                         z_w = f->zwsqs
325                                                         break
326                                                      else
327                                                         if(vNam(varn) .EQ. "wq")then
328                                                            z_w = f->zwq
329                                                            break
330                                                         else
331                                                            if(vNam(varn) .EQ. "wpqvp")then
332                                                               z_w = f->zwpqvp
333                                                               break
334                                                            else
335                                                               if(vNam(varn) .EQ. "wsqvs")then
336                                                                  z_w = f->zwsqvs
337                                                                  break
338                                                               else
339                                                                  if(vNam(varn) .EQ. "wqv")then
340                                                                     z_w = f->zwqv
341                                                                     break
342                                                                  else
343                                                                     if(vNam(varn) .EQ. "wpsp")then
344                                                                        z_w = f->zwpsp
345                                                                        break
346                                                                     else
347                                                                        if(vNam(varn) .EQ. "wsss")then
348                                                                           z_w = f->zwsss
349                                                                           break
350                                                                        else
351                                                                           if(vNam(varn) .EQ. "ws")then
352                                                                              z_w = f->zws
353                                                                              break
354                                                                           else
355                                                                              if(vNam(varn) .EQ. "wpsap")then
356                                                                                 z_w = f->zwpsap
357                                                                                 break
358                                                                              else
359                                                                                 if(vNam(varn) .EQ. "wssas")then
360                                                                                    z_w = f->zwssas
361                                                                                    break
362                                                                                 else
363                                                                                    if(vNam(varn) .EQ. "wsa")then
364                                                                                       z_w = f->zwsa
365                                                                                       break
366                                                                                    else
367                                                                                       if(vNam(varn) .EQ. "wses")then
368                                                                                          z_w = f->zwses
369                                                                                          break
370                                                                                       else
371                                                                                          if(vNam(varn) .EQ. "ws2")then
372                                                                                             z_w = f->zws2
373                                                                                             break
374                                                                                          else
375                                                                                             if(vNam(varn) .EQ. "ws3")then
376                                                                                                z_w = f->zws3
377                                                                                                break
378                                                                                             else
379                                                                                                if(vNam(varn) .EQ. "Sw")then
380                                                                                                   z_w = f->zSw
381                                                                                                   break
382                                                                                                else
383                                                                                                   if(vNam(varn) .EQ. "ws2pts")then
384                                                                                                      z_w = f->zws2pts
385                                                                                                      break
386                                                                                                   else
387                                                                                                      if(vNam(varn) .EQ. "wspts2")then
388                                                                                                         z_w = f->zwspts2
389                                                                                                         break                                           
390                                                                                                      end if
391                                                                                                   end if
392                                                                                                end if
393                                                                                             end if
394                                                                                          end if
395                                                                                       end if
396                                                                                    end if
397                                                                                 end if
398                                                                              end if
399                                                                           end if
400                                                                        end if   
401                                                                     end if
402                                                                  end if
403                                                               end if
404                                                            end if
405                                                         end if
406                                                      end if
407                                                   end if
408                                                end if
409                                             end if
410                                          end if
411                                       end if
412                                    end if
413                                 end if
414                              end if
415                           end if
416                        end if
417                     end if
418                  end if
419               end if
420            end if
421         end if
422      end if
423      varn=varn-1
424   end do
425   if ( .not. isvar("z_u") ) then
426      zu = 0
427   else
428      zu = 1
429   end if
430   if ( .not. isvar("z_w") ) then
431      zw = 0
432   else
433      zw = 1
434   end if
435
436   if (zu .EQ. 0 .AND. zw .EQ. 0) then
437      do varn=0,dim-1     
438         if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
439            check = False
440         else
441            if (.not. isvar("var")) then
442               check = True
443               if (parameter(21) .NE. "variables") then
444                  var=parameter(21)
445                  check = isStrSubset( var,","+vNam(varn)+"," )
446               end if
447            else         
448               check = isStrSubset( var,","+vNam(varn)+"," )
449            end if
450         end if
451         if (check)
452            z = f->$vNam(varn+1)$
453            dimz = dimsizes(z)
454            break
455         end if   
456      end do
457   end if
458   if (isvar("z_u") ) then
459      dimz  = dimsizes(z_u)
460   else
461      if (isvar("z_w"))then
462         dimz  = dimsizes(z_w)
463      end if
464   end if
465   t_all = f->time
466   nt    = dimsizes(t_all)
467   delta_t=t_all(nt-1)/nt
468 
469   ; ****************************************************       
470   ; start of time step and different types of mistakes that could be done
471   ; ****************************************************
472
473   if ( .not. isvar("start_time_step") ) then           
474      start_time_step=t_all(1)/3600
475      if (parameter(13) .NE. "t(0)") then
476         if (stringtodouble(parameter(13)) .GT. t_all(nt-1)/3600)
477            print(" ")
478            print("'start_time_step' = "+ parameter(13) +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
479            print(" ")
480            print("Please select another 'start_time_step'")
481            print(" ")
482            exit
483         end if
484         if (stringtofloat(parameter(13)) .LE. t_all(0)/3600)
485            print(" ")
486            print("'start_time_step' = "+ parameter(13) +"h is equal or lower than first time step = " + t_all(0)+"s = "+t_all(0)/3600+"h")
487            print(" ")
488            print("Please start profiles at least with second time step = "+t_all(1)+"s = "+t_all(1)/3600+"h")
489            print(" ")
490            exit
491         end if
492         start_time_step=stringtodouble(parameter(13))
493      end if
494   else
495      if (start_time_step .GT. t_all(nt-1)/3600)
496         print(" ")
497         print("'start_time_step' = "+ start_time_step +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
498         print(" ")
499         print("Please select another 'start_time_step'")
500         print(" ")
501         exit
502      end if
503      if (start_time_step .LE. t_all(0)/3600)
504         print(" ")
505         print("'start_time_step' = "+ start_time_step +"h is equal or lower than first time step = " + t_all(0)+"s = "+t_all(0)/3600+"h")
506         print(" ")
507         print("Please start profiles at least with second time step = "+t_all(1)+"s = "+t_all(1)/3600+"h")
508         print(" ")
509         exit
510      end if
511   end if
512   start_time_step = start_time_step*3600
513   if (start_time_step .GT. t_all(0) .AND. start_time_step .LT. t_all(1)+delta_t/2)then
514      st=1
515   end if
516   do i=2,nt-1     
517      if (start_time_step .GE. t_all(i)-delta_t/2 .AND. start_time_step .LT. t_all(i)+delta_t/2)then
518         st=i
519         break
520      end if
521   end do
522 
523   ; ****************************************************
524   ; end of time step and different types of mistakes that could be done
525   ; ****************************************************
526
527   if ( .not. isvar("end_time_step") ) then             
528      end_time_step = t_all(nt-1)/3600
529      if (parameter(15) .NE. "t(end)") then
530         if (stringtodouble(parameter(15)) .GT. t_all(nt-1)/3600)
531            print(" ")
532            print("'end_time_step' = "+ parameter(15) +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h")
533            print(" ")
534            print("Please select another 'end_time_step'") 
535            print(" ")
536            exit
537         end if
538         if (stringtodouble(parameter(15)) .LT. start_time_step/3600)
539            print(" ")
540            print("'end_time_step' = "+ parameter(15) +"h is lower than 'start_time_step' = "+parameter(13)+"h")
541            print(" ")
542            print("Please select another 'start_time_step' or 'end_time_step'")
543            print(" ")
544            exit
545         end if
546         end_time_step = stringtodouble(parameter(15))
547      end if   
548   else
549      if (end_time_step .GT. t_all(nt-1)/3600)
550         print(" ")
551         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")
552         print(" ")
553         print("Please select another 'end_time_step'") 
554         print(" ")
555         exit
556      end if
557      if (end_time_step .LT. start_time_step/3600)
558         print(" ")
559         print("'end_time_step' = "+ end_time_step +"h is lower than 'start_time_step' = "+start_time_step+"h")
560         print(" ")
561         print("Please select another 'start_time_step' or 'end_time_step'")
562         print(" ")
563         exit
564      end if
565   end if
566   end_time_step = end_time_step*3600
567   if (end_time_step .GT. t_all(0) .AND. end_time_step .LT. t_all(1)+delta_t/2)then
568      et=1
569   end if
570   do i=2,nt-1     
571      if (end_time_step .GE. t_all(i)-delta_t/2 .AND. end_time_step .LT. t_all(i)+delta_t/2)then
572         et=i
573         break
574      end if
575   end do
576
577   delete(start_time_step)
578   start_time_step=round(st,3)
579   delete(end_time_step)
580   end_time_step=round(et,3)
581
582   ; ****************************************************
583   ; set up legend and colors
584   ; ****************************************************
585   
586   legend_label=new(nt,double)
587   do p=start_time_step,end_time_step
588      legend_label(p-start_time_step)=t_all(p)
589   end do
590   
591   np = end_time_step-start_time_step+1 
592   if ( black .eq. 0 ) then
593      res=True     
594      res@xyLineColors = ispan(2,237,235/np)
595   end if
596
597   ; ***************************************************
598   ; set up recourses
599   ; ***************************************************
600
601   res@gsnDraw                 = False
602   res@gsnFrame                = False
603   res@gsnPaperOrientation     = "portrait"
604   res@gsnPaperWidth           = 8.27
605   res@gsnPaperHeight          = 11.69
606   res@gsnPaperMargin          = 0.79
607   res@txFont                  = "helvetica"
608   res@tiMainFont              = "helvetica"
609   res@tiXAxisFont             = "helvetica"
610   res@tiYAxisFont             = "helvetica"
611   res@tmXBLabelFont           = "helvetica"
612   res@tmYLLabelFont           = "helvetica"
613   res@lgLabelFont             = "helvetica"
614   res@tmLabelAutoStride       = True
615   res@pmLegendDisplayMode     = "Always"
616   res@pmLegendSide            = "Top"
617   res@xyExplicitLegendLabels  = legend_label
618   res@pmLegendParallelPosF    = 1.15
619   res@pmLegendOrthogonalPosF  = -1.0
620   res@pmLegendWidthF          = 0.12
621   res@pmLegendHeightF         = 0.3
622   res@lgLabelFontHeightF     = .02
623   res@txFontHeightF      = 0.02
624   res@tiXAxisFontHeightF = 0.02
625   res@tiYAxisFontHeightF = 0.02
626   res@tiXAxisString      = " "   
627   res@tiYAxisString      = "Height [z]"
628   
629   if ( dash .eq. 0 ) then
630      res@xyMonoDashPattern       = True 
631   end if
632
633   resP                        = True
634   resP@txFont                 = "helvetica"
635   resP@txString               = f@title
636   resP@txFuncCode             = "~"
637   resP@txFontHeightF          = 0.014
638
639   ; ***************************************************
640   ; set up graphics for plot
641   ; ***************************************************
642
643   plot = new(dim,graphic)
644   plot_ = new(dim,graphic)
645
646   if (combine .EQ. 1) then
647      if ( .not. isvar("number_comb") ) then           
648         if (parameter(25) .EQ. "0") then
649            print(" ")
650            print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot")
651            print(" ")
652            exit
653         else
654            number_comb=stringtointeger(parameter(25))
655            plot_o = new(number_comb,graphic)   
656         end if
657      else
658         if(number_comb .EQ. 2 .OR. number_comb .EQ. 3) then
659            plot_o = new(number_comb,graphic)     
660         else
661            print(" ")
662            print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot")
663            print(" ")
664            exit
665         end if   
666      end if
667      label=new(number_comb,string)
668      color_o=new(number_comb,integer)
669      mini=new(number_comb,float)
670      maxi=new(number_comb,float)
671   end if
672
673   wks=gsn_open_wks(format_out,file_out)
674   gsn_define_colormap(wks,"rainbow+white")
675 
676   ; ***************************************************
677   ; indicate plot number
678   ; ***************************************************
679   
680   if (combine .EQ. 1) then
681      n = 1
682   else
683      n = 0
684   end if
685
686   ; ***************************************************
687   ; set up minimum and maximum height
688   ; ***************************************************
689
690   if (.not. isvar("min_z"))
691      min_z=0
692      if (stringtointeger(parameter(33)) .NE. 0) then
693         if (stringtointeger(parameter(33)) .GE. max(z_u) ) then
694            print(" ")
695            print("Minimum of height ('min_z'="+stringtointeger(parameter(33))+") is greater than available heights (="+max(z_u)+")")
696            print(" ")
697            exit
698         end if
699         if (stringtointeger(parameter(33)) .LT. 0 ) then
700            print(" ")
701            print("Begin minimum of height 'min_z' with 0")
702            print(" ")
703            exit
704         end if
705         min_z=stringtointeger(parameter(33))
706      end if
707   else
708      if (min_z .GE. max(z_u) ) then
709         print(" ")
710         print("Minimum of height ('min_z'="+min_z+") is greater than available heights (="+max(z_u)+")")
711         print(" ")
712         exit
713      end if
714      if (min_z .LT. 0 ) then
715         print(" ")
716         print("Begin minimum of height 'min_z' with 0")
717         print(" ")
718         exit
719      end if
720   end if
721
722   if (.not. isvar("max_z"))
723      max_z=max(z_u)
724      if ((parameter(35)) .NE. "max(z_u)") then
725         if (stringtofloat(parameter(35)) .GE. max(z_u) ) then
726            print(" ")
727            print("Maximum of height ('max_z'="+parameter(35)+") is greater than available heights (="+max(z_u)+")")
728            print(" ")
729            exit
730         end if
731         if (stringtointeger(parameter(35)) .LE. 0 ) then
732            print(" ")
733            print("Maximum of height 'max_z' should be at least 1")
734            print(" ")
735            exit
736         end if
737         max_z=stringtointeger(parameter(35))
738      end if
739   else
740      if (max_z .GE. max(z_pr) ) then
741         print(" ")
742         print("Maximum of height ('max_z'="+max_z+") is greater than available heights (="+max(z_u)+")")
743         print(" ")
744         exit
745      end if
746      if (max_z .LE. 0 ) then
747         print(" ")
748         print("Maximum of height 'max_z' should be at least 1")
749         print(" ")
750         exit
751      end if
752   end if
753
754   ; ***************************************************
755   ; read data and create plots
756   ; ***************************************************
757     
758   do ti = start_time_step, end_time_step
759      if( t_all(ti) .lt. 10^36) then
760         start_time_step = ti
761         break
762      end if
763   end do 
764 
765   data   = new((/dim,(end_time_step-start_time_step)+1,dimz/),float)
766   data_0 = new((/(end_time_step-start_time_step)+1,dimz/),float)
767   data_0 = 0.0
768   t      = new((/(end_time_step-start_time_step)+1,dimz/),float)
769   t      = 0.0
770   unit   = new(dim,string)
771
772   if (over .EQ. 1) then
773      plot_u         = gsn_csm_xy(wks,t,data_0(:,:),res)
774      miniu = 1.E27
775      maxiu =-1.E27
776      plot_v         = gsn_csm_xy(wks,t,data_0(:,:),res)
777      miniv = 1.E27
778      maxiv =-1.E27
779      plot_w         = gsn_csm_xy(wks,t,data_0(:,:),res)
780      miniw = 1.E27
781      maxiw =-1.E27
782      plot_pt        = gsn_csm_xy(wks,t,data_0(:,:),res)
783      minipt = 1.E27
784      maxipt =-1.E27
785      plot_vpt       = gsn_csm_xy(wks,t,data_0(:,:),res)
786      minivpt = 1.E27
787      maxivpt =-1.E27
788      plot_lpt       = gsn_csm_xy(wks,t,data_0(:,:),res)
789      minilpt = 1.E27
790      maxilpt =-1.E27
791      plot_q         = gsn_csm_xy(wks,t,data_0(:,:),res)
792      miniq = 1.E27
793      maxiq =-1.E27
794      plot_qv        = gsn_csm_xy(wks,t,data_0(:,:),res)
795      miniqv = 1.E27
796      maxiqv =-1.E27
797      plot_ql        = gsn_csm_xy(wks,t,data_0(:,:),res)
798      miniql = 1.E27
799      maxiql =-1.E27
800      plot_rho       = gsn_csm_xy(wks,t,data_0(:,:),res)
801      plot_s         = gsn_csm_xy(wks,t,data_0(:,:),res)
802      plot_sa        = gsn_csm_xy(wks,t,data_0(:,:),res)
803      plot_e         = gsn_csm_xy(wks,t,data_0(:,:),res)
804      minie = 1.E27
805      maxie =-1.E27
806      plot_es        = gsn_csm_xy(wks,t,data_0(:,:),res)
807      minie = 1.E27
808      maxie =-1.E27
809      plot_km        = gsn_csm_xy(wks,t,data_0(:,:),res)
810      minie = 1.E27
811      maxie =-1.E27
812      plot_kh        = gsn_csm_xy(wks,t,data_0(:,:),res)
813      minie = 1.E27
814      maxie =-1.E27
815      plot_l         = gsn_csm_xy(wks,t,data_0(:,:),res)     
816      plot_wpup      = gsn_csm_xy(wks,t,data_0(:,:),res)
817      miniwpup = 1.E27
818      maxiwpup =-1.E27
819      plot_wsus      = gsn_csm_xy(wks,t,data_0(:,:),res)
820      miniwsus = 1.E27
821      maxiwsus =-1.E27
822      plot_wu        = gsn_csm_xy(wks,t,data_0(:,:),res)
823      miniwu = 1.E27
824      maxiwu =-1.E27
825      plot_wpvp      = gsn_csm_xy(wks,t,data_0(:,:),res)
826      miniwpvp = 1.E27
827      maxiwpvp =-1.E27
828      plot_wsvs      = gsn_csm_xy(wks,t,data_0(:,:),res)
829      miniwsvs = 1.E27
830      maxiwsvs =-1.E27
831      plot_wv        = gsn_csm_xy(wks,t,data_0(:,:),res)
832      miniwv = 1.E27
833      maxiwv =-1.E27
834      plot_wpptp     = gsn_csm_xy(wks,t,data_0(:,:),res)
835      miniwpptp = 1.E27
836      maxiwpptp =-1.E27
837      plot_wspts     = gsn_csm_xy(wks,t,data_0(:,:),res)
838      miniwspts = 1.E27
839      maxiwspts =-1.E27
840      plot_wpt       = gsn_csm_xy(wks,t,data_0(:,:),res)
841      miniwpt = 1.E27
842      maxiwpt =-1.E27
843      plot_wsptsBC   = gsn_csm_xy(wks,t,data_0(:,:),res)
844      miniwsptsBC = 1.E27
845      maxiwsptsBC =-1.E27
846      plot_wptBC     = gsn_csm_xy(wks,t,data_0(:,:),res)
847      miniwptBC = 1.E27
848      maxiwptBC =-1.E27
849      plot_wpvptp    = gsn_csm_xy(wks,t,data_0(:,:),res)
850      miniwpvptp = 1.E27
851      maxiwpvptp =-1.E27
852      plot_wsvpts    = gsn_csm_xy(wks,t,data_0(:,:),res)
853      miniwsvpts = 1.E27
854      maxiewsvpts=-1.E27
855      plot_wvpt      = gsn_csm_xy(wks,t,data_0(:,:),res)
856      miniwvpt = 1.E27
857      maxiwvpt =-1.E27
858      plot_wpqp      = gsn_csm_xy(wks,t,data_0(:,:),res)
859      miniwpqp = 1.E27
860      maxiwpqp =-1.E27
861      plot_wsqs      = gsn_csm_xy(wks,t,data_0(:,:),res)
862      miniwsqs = 1.E27
863      maxiwsqs =-1.E27
864      plot_wq        = gsn_csm_xy(wks,t,data_0(:,:),res)
865      miniwq = 1.E27
866      maxiwq =-1.E27
867      plot_wpqvp     = gsn_csm_xy(wks,t,data_0(:,:),res)
868      miniwpqvp = 1.E27
869      maxiwpqvp =-1.E27
870      plot_wsqvs     = gsn_csm_xy(wks,t,data_0(:,:),res)
871      miniwsqvs = 1.E27
872      maxiwsqvs =-1.E27
873      plot_wqv       = gsn_csm_xy(wks,t,data_0(:,:),res)
874      miniwqv = 1.E27
875      maxiwqv =-1.E27
876      plot_wpsp      = gsn_csm_xy(wks,t,data_0(:,:),res)
877      miniwpsp = 1.E27
878      maxiwpsp =-1.E27
879      plot_wsss      = gsn_csm_xy(wks,t,data_0(:,:),res)
880      miniwsss = 1.E27
881      maxiwsss =-1.E27
882      plot_ws        = gsn_csm_xy(wks,t,data_0(:,:),res)
883      miniws = 1.E27
884      maxiws =-1.E27
885      plot_wpsap     = gsn_csm_xy(wks,t,data_0(:,:),res)
886      miniwpsap = 1.E27
887      maxiwpsap =-1.E27
888      plot_wssas     = gsn_csm_xy(wks,t,data_0(:,:),res)
889      miniwssas = 1.E27
890      maxiwssas =-1.E27
891      plot_wsa       = gsn_csm_xy(wks,t,data_0(:,:),res)
892      miniwsa = 1.E27
893      maxiwsa =-1.E27
894      plot_wses      = gsn_csm_xy(wks,t,data_0(:,:),res)
895      plot_us2       = gsn_csm_xy(wks,t,data_0(:,:),res)
896      minius2 = 1.E27
897      maxius2 =-1.E27
898      plot_vs2       = gsn_csm_xy(wks,t,data_0(:,:),res)
899      minivs2 = 1.E27
900      maxivs2 =-1.E27
901      plot_ws2       = gsn_csm_xy(wks,t,data_0(:,:),res)
902      miniws2 = 1.E27
903      maxiws2 =-1.E27
904      plot_pts2      = gsn_csm_xy(wks,t,data_0(:,:),res)
905      plot_ws3       = gsn_csm_xy(wks,t,data_0(:,:),res)
906      plot_Sw        = gsn_csm_xy(wks,t,data_0(:,:),res)
907      plot_ws2pts    = gsn_csm_xy(wks,t,data_0(:,:),res)
908      plot_wspts2    = gsn_csm_xy(wks,t,data_0(:,:),res)
909      plot_wsususodz = gsn_csm_xy(wks,t,data_0(:,:),res)
910      miniwsususodz = 1.E27
911      maxiwsususodz =-1.E27
912      plot_wspsodz   = gsn_csm_xy(wks,t,data_0(:,:),res)
913      miniwspsodz = 1.E27
914      maxiwspsodz =-1.E27
915      plot_wpeodz    = gsn_csm_xy(wks,t,data_0(:,:),res)
916      miniwpeodz = 1.E27
917      maxiwpeodz =-1.E27
918   end if
919 
920   n_o=0
921   count_var=0
922   do varn = 0,dim-1
923       
924      temp = f->$vNam(varn)$
925     
926      if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
927         check = False
928      else
929         if (.not. isvar("var")) then
930            check = True
931            if (parameter(21) .NE. "variables") then
932               var=parameter(21)
933               check = isStrSubset( var,","+vNam(varn)+"," )
934            end if
935         else         
936            check = isStrSubset( var,","+vNam(varn)+"," )
937         end if
938      end if
939
940      if (combine .EQ. 1) then         
941         com=isStrSubset(c_var,","+vNam(varn)+"," )     
942         if (com)                       
943            data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1)
944            unit(varn) = temp@units
945            if (n_o .GT. number_comb-1) then
946               print(" ")
947               print("Please set 'number_comb' to the right number of overlaying variables ('c_var')")
948               print(" ")
949               exit
950            end if
951            mini(n_o)=min(data(varn,:,:))
952            maxi(n_o)=max(data(varn,:,:))
953            n_o=n_o+1
954         end if
955      end if
956
957      if(check) then
958         count_var=count_var+1
959     
960         z = f->$vNam(varn+1)$
961 
962         unit(varn) = temp@units
963         data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1)
964 
965         if (over .EQ. 0) then 
966            res@gsnLeftString      = vNam(varn)
967            res@gsnRightString     = unit(varn)
968            res@trYMinF            = min_z
969            res@trYMaxF            = max_z
970            if (.not. isvar("xs")) then
971               if (parameter(63) .NE. "x0") then
972                  res@trXMinF = stringtofloat(parameter(63))
973               else
974                  res@trXMinF            = min(data(varn,:,:))
975               end if
976            else
977               res@trXMinF            = xs
978            end if
979            if (.not. isvar("xe")) then
980               if (parameter(65) .NE. "xdim") then
981                  res@trXMaxF = stringtofloat(parameter(65))
982               else
983                  res@trXMaxF            = max(data(varn,:,:))
984               end if
985            else
986               res@trXMaxF            = xe 
987            end if                     
988            plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 
989         end if       
990         
991            if (vNam(varn) .EQ. "u") then
992               miniu=min(data(varn,:,:))
993               maxiu=max(data(varn,:,:))
994               if (over .EQ. 1) then
995                  res@xyDashPattern  = 0
996                  plot_u = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
997               else
998                  res@gsnLeftString      = vNam(varn)
999                  res@gsnRightString     = unit(varn)
1000                  if (.not. isvar("xs")) then
1001                     if (parameter(63) .NE. "x0") then
1002                        res@trXMinF = stringtofloat(parameter(63))
1003                     else
1004                        res@trXMinF            = miniu
1005                     end if
1006                  else
1007                     res@trXMinF            = xs
1008                  end if
1009                  if (.not. isvar("xe")) then
1010                     if (parameter(65) .NE. "xdim") then
1011                        res@trXMaxF = stringtofloat(parameter(65))
1012                     else
1013                        res@trXMaxF            = maxiu
1014                     end if
1015                  else
1016                     res@trXMaxF            = xe 
1017                  end if               
1018                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 
1019               end if
1020            end if
1021            if (vNam(varn) .EQ. "v") then
1022               miniv=min(data(varn,:,:))
1023               maxiv=max(data(varn,:,:))
1024               if (over .EQ. 1) then
1025                  res@xyDashPattern  = 1
1026                  plot_v = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1027               else
1028                  res@gsnLeftString      = vNam(varn)
1029                  res@gsnRightString     = unit(varn)
1030                  if (.not. isvar("xs")) then
1031                     if (parameter(63) .NE. "x0") then
1032                        res@trXMinF = stringtofloat(parameter(63))
1033                     else
1034                        res@trXMinF            = miniv
1035                     end if
1036                  else
1037                     res@trXMinF            = xs
1038                  end if
1039                  if (.not. isvar("xe")) then
1040                     if (parameter(65) .NE. "xdim") then
1041                        res@trXMaxF = stringtofloat(parameter(65))
1042                     else
1043                        res@trXMaxF            = maxiv
1044                     end if
1045                  else
1046                     res@trXMaxF            = xe 
1047                  end if               
1048                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1049               end if 
1050            end if
1051            if (vNam(varn) .EQ. "w") then
1052               miniw=min(data(varn,:,:))
1053               maxiw=max(data(varn,:,:))
1054               if (over .EQ. 1) then
1055                  res@xyDashPattern  = 2
1056                  plot_w = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1057               else
1058                  res@gsnLeftString      = vNam(varn)
1059                  res@gsnRightString     = unit(varn)
1060                  if (.not. isvar("xs")) then
1061                     if (parameter(63) .NE. "x0") then
1062                        res@trXMinF = stringtofloat(parameter(63))
1063                     else
1064                        res@trXMinF            = miniw
1065                     end if
1066                  else
1067                     res@trXMinF            = xs
1068                  end if
1069                  if (.not. isvar("xe")) then
1070                     if (parameter(65) .NE. "xdim") then
1071                        res@trXMaxF = stringtofloat(parameter(65))
1072                     else
1073                        res@trXMaxF            = maxiw
1074                     end if
1075                  else
1076                     res@trXMaxF            = xe 
1077                  end if           
1078                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1079               end if   
1080            end if
1081
1082            if (vNam(varn) .EQ. "pt") then
1083               minipt=min(data(varn,:,:))
1084               maxipt=max(data(varn,:,:))
1085               if (over .EQ. 1) then
1086                  res@xyDashPattern  = 0
1087                  plot_pt = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1088               else
1089                  res@gsnLeftString      = vNam(varn)
1090                  res@gsnRightString     = unit(varn)
1091                  if (.not. isvar("xs")) then
1092                     if (parameter(63) .NE. "x0") then
1093                        res@trXMinF = stringtofloat(parameter(63))
1094                     else
1095                        res@trXMinF            = minipt
1096                     end if
1097                  else
1098                     res@trXMinF            = xs
1099                  end if
1100                  if (.not. isvar("xe")) then
1101                     if (parameter(65) .NE. "xdim") then
1102                        res@trXMaxF = stringtofloat(parameter(65))
1103                     else
1104                        res@trXMaxF            = maxipt
1105                     end if
1106                  else
1107                     res@trXMaxF            = xe 
1108                  end if
1109                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1110               end if 
1111            end if
1112            if (vNam(varn) .EQ. "vpt") then
1113               minivpt=min(data(varn,:,:))
1114               maxivpt=max(data(varn,:,:))
1115               if (over .EQ. 1) then
1116                  res@xyDashPattern  = 1
1117                  plot_vpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1118               else
1119                  res@gsnLeftString      = vNam(varn)
1120                  res@gsnRightString     = unit(varn)
1121                  if (.not. isvar("xs")) then
1122                     if (parameter(63) .NE. "x0") then
1123                        res@trXMinF = stringtofloat(parameter(63))
1124                     else
1125                        res@trXMinF            = minivpt
1126                     end if
1127                  else
1128                     res@trXMinF            = xs
1129                  end if
1130                  if (.not. isvar("xe")) then
1131                     if (parameter(65) .NE. "xdim") then
1132                        res@trXMaxF = stringtofloat(parameter(65))
1133                     else
1134                        res@trXMaxF            = maxivpt
1135                     end if
1136                  else
1137                     res@trXMaxF            = xe 
1138                  end if
1139                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1140               end if
1141            end if
1142            if (vNam(varn) .EQ. "lpt") then
1143               minilpt=min(data(varn,:,:))
1144               maxilpt=max(data(varn,:,:))
1145               if (over .EQ. 1) then
1146                  res@xyDashPattern  = 2
1147                  plot_lpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1148               else
1149                  res@gsnLeftString      = vNam(varn)
1150                  res@gsnRightString     = unit(varn)
1151                  if (.not. isvar("xs")) then
1152                     if (parameter(63) .NE. "x0") then
1153                        res@trXMinF = stringtofloat(parameter(63))
1154                     else
1155                        res@trXMinF            = minilpt
1156                     end if
1157                  else
1158                     res@trXMinF            = xs
1159                  end if
1160                  if (.not. isvar("xe")) then
1161                     if (parameter(65) .NE. "xdim") then
1162                        res@trXMaxF = stringtofloat(parameter(65))
1163                     else
1164                        res@trXMaxF            = maxilpt
1165                     end if
1166                  else
1167                     res@trXMaxF            = xe 
1168                  end if
1169                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1170               end if
1171            end if
1172
1173            if (vNam(varn) .EQ. "q") then
1174               miniq=min(data(varn,:,:))
1175               maxiq=max(data(varn,:,:))
1176               if (over .EQ. 1) then
1177                  res@xyDashPattern  = 0
1178                  plot_q = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1179               else
1180                  res@gsnLeftString      = vNam(varn)
1181                  res@gsnRightString     = unit(varn)
1182                  if (.not. isvar("xs")) then
1183                     if (parameter(63) .NE. "x0") then
1184                        res@trXMinF = stringtofloat(parameter(63))
1185                     else
1186                        res@trXMinF            = minilq
1187                     end if
1188                  else
1189                     res@trXMinF            = xs
1190                  end if
1191                  if (.not. isvar("xe")) then
1192                     if (parameter(65) .NE. "xdim") then
1193                        res@trXMaxF = stringtofloat(parameter(65))
1194                     else
1195                        res@trXMaxF            = maxilq
1196                     end if
1197                  else
1198                     res@trXMaxF            = xe 
1199                  end if
1200                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1201               end if
1202            end if
1203            if (vNam(varn) .EQ. "qv") then
1204               miniqv=min(data(varn,:,:))
1205               maxiqv=max(data(varn,:,:))
1206               if (over .EQ. 1) then
1207                  res@xyDashPattern  = 1
1208                  plot_qv = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1209               else
1210                  res@gsnLeftString      = vNam(varn)
1211                  res@gsnRightString     = unit(varn)
1212                  if (.not. isvar("xs")) then
1213                     if (parameter(63) .NE. "x0") then
1214                        res@trXMinF = stringtofloat(parameter(63))
1215                     else
1216                        res@trXMinF            = minilqv
1217                     end if
1218                  else
1219                     res@trXMinF            = xs
1220                  end if
1221                  if (.not. isvar("xe")) then
1222                     if (parameter(65) .NE. "xdim") then
1223                        res@trXMaxF = stringtofloat(parameter(65))
1224                     else
1225                        res@trXMaxF            = maxilqv
1226                     end if
1227                  else
1228                     res@trXMaxF            = xe 
1229                  end if
1230                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1231               end if
1232            end if
1233            if (vNam(varn) .EQ. "ql") then
1234               miniql=min(data(varn,:,:))
1235               maxiql=max(data(varn,:,:))
1236               if (over .EQ. 1) then
1237                  res@xyDashPattern  = 2
1238                  plot_ql = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1239               else
1240                  res@gsnLeftString      = vNam(varn)
1241                  res@gsnRightString     = unit(varn)
1242                  if (.not. isvar("xs")) then
1243                     if (parameter(63) .NE. "x0") then
1244                        res@trXMinF = stringtofloat(parameter(63))
1245                     else
1246                        res@trXMinF            = miniql
1247                     end if
1248                  else
1249                     res@trXMinF            = xs
1250                  end if
1251                  if (.not. isvar("xe")) then
1252                     if (parameter(65) .NE. "xdim") then
1253                        res@trXMaxF = stringtofloat(parameter(65))
1254                     else
1255                        res@trXMaxF            = maxiql
1256                     end if
1257                  else
1258                     res@trXMaxF            = xe 
1259                  end if
1260                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1261               end if
1262            end if
1263
1264            if (vNam(varn) .EQ. "e") then
1265               minie=min(data(varn,:,:))
1266               maxie=max(data(varn,:,:))
1267               if (over .EQ. 1) then
1268                  res@xyDashPattern  = 0
1269                  plot_e = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1270               else
1271                  res@gsnLeftString      = vNam(varn)
1272                  res@gsnRightString     = unit(varn)
1273                  if (.not. isvar("xs")) then
1274                     if (parameter(63) .NE. "x0") then
1275                        res@trXMinF = stringtofloat(parameter(63))
1276                     else
1277                        res@trXMinF            = minie
1278                     end if
1279                  else
1280                     res@trXMinF            = xs
1281                  end if
1282                  if (.not. isvar("xe")) then
1283                     if (parameter(65) .NE. "xdim") then
1284                        res@trXMaxF = stringtofloat(parameter(65))
1285                     else
1286                        res@trXMaxF            = maxie
1287                     end if
1288                  else
1289                     res@trXMaxF            = xe 
1290                  end if
1291                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1292               end if
1293            end if
1294            if (vNam(varn) .EQ. "es") then
1295               minies=min(data(varn,:,:))
1296               maxies=max(data(varn,:,:))
1297               if (over .EQ. 1) then
1298                  res@xyDashPattern  = 1
1299                  plot_es = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1300               else
1301                  res@gsnLeftString      = vNam(varn)
1302                  res@gsnRightString     = unit(varn)
1303                  if (.not. isvar("xs")) then
1304                     if (parameter(63) .NE. "x0") then
1305                        res@trXMinF = stringtofloat(parameter(63))
1306                     else
1307                        res@trXMinF            = minies
1308                     end if
1309                  else
1310                     res@trXMinF            = xs
1311                  end if
1312                  if (.not. isvar("xe")) then
1313                     if (parameter(65) .NE. "xdim") then
1314                        res@trXMaxF = stringtofloat(parameter(65))
1315                     else
1316                        res@trXMaxF            = maxies
1317                     end if
1318                  else
1319                     res@trXMaxF            = xe 
1320                  end if
1321                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1322               end if
1323            end if
1324
1325            if (vNam(varn) .EQ. "km") then
1326               minikm=min(data(varn,:,:))
1327               maxikm=max(data(varn,:,:))
1328               if (over .EQ. 1) then
1329                  res@xyDashPattern  = 0
1330                  plot_km = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1331               else
1332                  res@gsnLeftString      = vNam(varn)
1333                  res@gsnRightString     = unit(varn)
1334                  if (.not. isvar("xs")) then
1335                     if (parameter(63) .NE. "x0") then
1336                        res@trXMinF = stringtofloat(parameter(63))
1337                     else
1338                        res@trXMinF            = minikm
1339                     end if
1340                  else
1341                     res@trXMinF            = xs
1342                  end if
1343                  if (.not. isvar("xe")) then
1344                     if (parameter(65) .NE. "xdim") then
1345                        res@trXMaxF = stringtofloat(parameter(65))
1346                     else
1347                        res@trXMaxF            = maxikm
1348                     end if
1349                  else
1350                     res@trXMaxF            = xe 
1351                  end if
1352                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1353               end if
1354            end if
1355            if (vNam(varn) .EQ. "kh") then
1356               minikh=min(data(varn,:,:))
1357               maxikh=max(data(varn,:,:))
1358               if (over .EQ. 1) then
1359                  res@xyDashPattern  = 1
1360                  plot_kh = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1361               else
1362                  res@gsnLeftString      = vNam(varn)
1363                  res@gsnRightString     = unit(varn)
1364                  if (.not. isvar("xs")) then
1365                     if (parameter(63) .NE. "x0") then
1366                        res@trXMinF = stringtofloat(parameter(63))
1367                     else
1368                        res@trXMinF            = minikh
1369                     end if
1370                  else
1371                     res@trXMinF            = xs
1372                  end if
1373                  if (.not. isvar("xe")) then
1374                     if (parameter(65) .NE. "xdim") then
1375                        res@trXMaxF = stringtofloat(parameter(65))
1376                     else
1377                        res@trXMaxF            = maxikh
1378                     end if
1379                  else
1380                     res@trXMaxF            = xe 
1381                  end if
1382                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
1383               end if
1384            end if
1385
1386            if (vNam(varn) .EQ. "wpup") then
1387               miniwpup=min(data(varn,:,:))
1388               maxiwpup=max(data(varn,:,:))
1389               if (over .EQ. 1) then
1390                  res@xyDashPattern  = 0
1391                  plot_wpup = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1392               else
1393                  res@gsnLeftString      = vNam(varn)
1394                  res@gsnRightString     = unit(varn)
1395                  if (.not. isvar("xs")) then
1396                     if (parameter(63) .NE. "x0") then
1397                        res@trXMinF = stringtofloat(parameter(63))
1398                     else
1399                        res@trXMinF            = miniwpup
1400                     end if
1401                  else
1402                     res@trXMinF            = xs
1403                  end if
1404                  if (.not. isvar("xe")) then
1405                     if (parameter(65) .NE. "xdim") then
1406                        res@trXMaxF = stringtofloat(parameter(65))
1407                     else
1408                        res@trXMaxF            = maxiwpup
1409                     end if
1410                  else
1411                     res@trXMaxF            = xe 
1412                  end if
1413                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1414               end if
1415            end if
1416            if (vNam(varn) .EQ. "wsus") then
1417               miniwsus=min(data(varn,:,:))
1418               maxiwsus=max(data(varn,:,:))
1419               if (over .EQ. 1) then
1420                  res@xyDashPattern  = 1
1421                  plot_wsus = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1422               else
1423                  res@gsnLeftString      = vNam(varn)
1424                  res@gsnRightString     = unit(varn)
1425                  if (.not. isvar("xs")) then
1426                     if (parameter(63) .NE. "x0") then
1427                        res@trXMinF = stringtofloat(parameter(63))
1428                     else
1429                        res@trXMinF            = miniwsus
1430                     end if
1431                  else
1432                     res@trXMinF            = xs
1433                  end if
1434                  if (.not. isvar("xe")) then
1435                     if (parameter(65) .NE. "xdim") then
1436                        res@trXMaxF = stringtofloat(parameter(65))
1437                     else
1438                        res@trXMaxF            = maxiwsus
1439                     end if
1440                  else
1441                     res@trXMaxF            = xe 
1442                  end if
1443                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1444               end if
1445            end if
1446            if (vNam(varn) .EQ. "wu") then
1447               miniwu=min(data(varn,:,:))
1448               maxiwu=max(data(varn,:,:))
1449               if (over .EQ. 1) then
1450                  res@xyDashPattern  = 2
1451                  plot_wu = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1452               else
1453                  res@gsnLeftString      = vNam(varn)
1454                  res@gsnRightString     = unit(varn)
1455                  if (.not. isvar("xs")) then
1456                     if (parameter(63) .NE. "x0") then
1457                        res@trXMinF = stringtofloat(parameter(63))
1458                     else
1459                        res@trXMinF            = miniwu
1460                     end if
1461                  else
1462                     res@trXMinF            = xs
1463                  end if
1464                  if (.not. isvar("xe")) then
1465                     if (parameter(65) .NE. "xdim") then
1466                        res@trXMaxF = stringtofloat(parameter(65))
1467                     else
1468                        res@trXMaxF            = maxiwu
1469                     end if
1470                  else
1471                     res@trXMaxF            = xe 
1472                  end if
1473                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1474               end if
1475            end if
1476
1477            if (vNam(varn) .EQ. "wpvp") then
1478               miniwpvp=min(data(varn,:,:))
1479               maxiwpvp=max(data(varn,:,:))
1480               if (over .EQ. 1) then
1481                  res@xyDashPattern  = 0
1482                  plot_wpvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1483               else
1484                  res@gsnLeftString      = vNam(varn)
1485                  res@gsnRightString     = unit(varn)
1486                  if (.not. isvar("xs")) then
1487                     if (parameter(63) .NE. "x0") then
1488                        res@trXMinF = stringtofloat(parameter(63))
1489                     else
1490                        res@trXMinF            = miniwpvp
1491                     end if
1492                  else
1493                     res@trXMinF            = xs
1494                  end if
1495                  if (.not. isvar("xe")) then
1496                     if (parameter(65) .NE. "xdim") then
1497                        res@trXMaxF = stringtofloat(parameter(65))
1498                     else
1499                        res@trXMaxF            = maxiwpvp
1500                     end if
1501                  else
1502                     res@trXMaxF            = xe 
1503                  end if
1504                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1505               end if
1506            end if
1507            if (vNam(varn) .EQ. "wsvs") then
1508               miniwsvs=min(data(varn,:,:))
1509               maxiwsvs=max(data(varn,:,:))
1510               if (over .EQ. 1) then
1511                  res@xyDashPattern  = 1
1512                  plot_wsvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1513               else
1514                  res@gsnLeftString      = vNam(varn)
1515                  res@gsnRightString     = unit(varn)
1516                  if (.not. isvar("xs")) then
1517                     if (parameter(63) .NE. "x0") then
1518                        res@trXMinF = stringtofloat(parameter(63))
1519                     else
1520                        res@trXMinF            = miniwsvs
1521                     end if
1522                  else
1523                     res@trXMinF            = xs
1524                  end if
1525                  if (.not. isvar("xe")) then
1526                     if (parameter(65) .NE. "xdim") then
1527                        res@trXMaxF = stringtofloat(parameter(65))
1528                     else
1529                        res@trXMaxF            = maxiwsvs
1530                     end if
1531                  else
1532                     res@trXMaxF            = xe 
1533                  end if
1534                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1535               end if
1536            end if
1537            if (vNam(varn) .EQ. "wv") then
1538               miniwv=min(data(varn,:,:))
1539               maxiwv=max(data(varn,:,:))
1540               if (over .EQ. 1) then
1541                  res@xyDashPattern  = 2
1542                  plot_wv = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1543               else
1544                  res@gsnLeftString      = vNam(varn)
1545                  res@gsnRightString     = unit(varn)
1546                  if (.not. isvar("xs")) then
1547                     if (parameter(63) .NE. "x0") then
1548                        res@trXMinF = stringtofloat(parameter(63))
1549                     else
1550                        res@trXMinF            = miniwv
1551                     end if
1552                  else
1553                     res@trXMinF            = xs
1554                  end if
1555                  if (.not. isvar("xe")) then
1556                     if (parameter(65) .NE. "xdim") then
1557                        res@trXMaxF = stringtofloat(parameter(65))
1558                     else
1559                        res@trXMaxF            = maxiwv
1560                     end if
1561                  else
1562                     res@trXMaxF            = xe 
1563                  end if
1564                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1565               end if
1566            end if
1567
1568            if (vNam(varn) .EQ. "wpptp") then
1569               miniwpptp=min(data(varn,:,:))
1570               maxiwpptp=max(data(varn,:,:))
1571               if (over .EQ. 1) then
1572                  res@xyDashPattern  = 0
1573                  plot_wpptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1574               else
1575                  res@gsnLeftString      = vNam(varn)
1576                  res@gsnRightString     = unit(varn)
1577                  if (.not. isvar("xs")) then
1578                     if (parameter(63) .NE. "x0") then
1579                        res@trXMinF = stringtofloat(parameter(63))
1580                     else
1581                        res@trXMinF            = miniwpptp
1582                     end if
1583                  else
1584                     res@trXMinF            = xs
1585                  end if
1586                  if (.not. isvar("xe")) then
1587                     if (parameter(65) .NE. "xdim") then
1588                        res@trXMaxF = stringtofloat(parameter(65))
1589                     else
1590                        res@trXMaxF            = maxiwpptp
1591                     end if
1592                  else
1593                     res@trXMaxF            = xe 
1594                  end if
1595                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1596               end if
1597            end if
1598            if (vNam(varn) .EQ. "wspts") then
1599               miniwspts=min(data(varn,:,:))
1600               maxiwspts=max(data(varn,:,:))
1601               if (over .EQ. 1) then
1602                  res@xyDashPattern  = 1
1603                  plot_wspts = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1604               else
1605                  res@gsnLeftString      = vNam(varn)
1606                  res@gsnRightString     = unit(varn)
1607                  if (.not. isvar("xs")) then
1608                     if (parameter(63) .NE. "x0") then
1609                        res@trXMinF = stringtofloat(parameter(63))
1610                     else
1611                        res@trXMinF            = miniwspts
1612                     end if
1613                  else
1614                     res@trXMinF            = xs
1615                  end if
1616                  if (.not. isvar("xe")) then
1617                     if (parameter(65) .NE. "xdim") then
1618                        res@trXMaxF = stringtofloat(parameter(65))
1619                     else
1620                        res@trXMaxF            = maxiwspts
1621                     end if
1622                  else
1623                     res@trXMaxF            = xe 
1624                  end if
1625                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1626               end if
1627            end if
1628            if (vNam(varn) .EQ. "wpt") then
1629               miniwpt=min(data(varn,:,:))
1630               maxiwpt=max(data(varn,:,:))
1631               if (over .EQ. 1) then
1632                  res@xyDashPattern  = 2
1633                  plot_wpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1634               else
1635                  res@gsnLeftString      = vNam(varn)
1636                  res@gsnRightString     = unit(varn)
1637                  if (.not. isvar("xs")) then
1638                     if (parameter(63) .NE. "x0") then
1639                        res@trXMinF = stringtofloat(parameter(63))
1640                     else
1641                        res@trXMinF            = miniwpt
1642                     end if
1643                  else
1644                     res@trXMinF            = xs
1645                  end if
1646                  if (.not. isvar("xe")) then
1647                     if (parameter(65) .NE. "xdim") then
1648                        res@trXMaxF = stringtofloat(parameter(65))
1649                     else
1650                        res@trXMaxF            = maxiwpt
1651                     end if
1652                  else
1653                     res@trXMaxF            = xe 
1654                  end if
1655                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1656               end if
1657            end if
1658
1659            if (vNam(varn) .EQ. "wsptsBC") then
1660               miniwsptsBC=min(data(varn,:,:))
1661               maxiwsptsBC=max(data(varn,:,:))
1662               if (over .EQ. 1) then
1663                  res@xyDashPattern  = 0
1664                  plot_wsptsBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1665               else
1666                  res@gsnLeftString      = vNam(varn)
1667                  res@gsnRightString     = unit(varn)
1668                  if (.not. isvar("xs")) then
1669                     if (parameter(63) .NE. "x0") then
1670                        res@trXMinF = stringtofloat(parameter(63))
1671                     else
1672                        res@trXMinF            = miniwsptsBC
1673                     end if
1674                  else
1675                     res@trXMinF            = xs
1676                  end if
1677                  if (.not. isvar("xe")) then
1678                     if (parameter(65) .NE. "xdim") then
1679                        res@trXMaxF = stringtofloat(parameter(65))
1680                     else
1681                        res@trXMaxF            = maxiwsptsBC
1682                     end if
1683                  else
1684                     res@trXMaxF            = xe 
1685                  end if
1686                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1687               end if
1688            end if             
1689            if (vNam(varn) .EQ. "wptBC") then
1690               miniwptBC=min(data(varn,:,:))
1691               maxiwptBC=max(data(varn,:,:))
1692               if (over .EQ. 1) then
1693                  res@xyDashPattern  = 1
1694                  plot_wptBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1695               else
1696                  res@gsnLeftString      = vNam(varn)
1697                  res@gsnRightString     = unit(varn)
1698                  if (.not. isvar("xs")) then
1699                     if (parameter(63) .NE. "x0") then
1700                        res@trXMinF = stringtofloat(parameter(63))
1701                     else
1702                        res@trXMinF            = miniwptBC
1703                     end if
1704                  else
1705                     res@trXMinF            = xs
1706                  end if
1707                  if (.not. isvar("xe")) then
1708                     if (parameter(65) .NE. "xdim") then
1709                        res@trXMaxF = stringtofloat(parameter(65))
1710                     else
1711                        res@trXMaxF            = maxiwptBC
1712                     end if
1713                  else
1714                     res@trXMaxF            = xe 
1715                  end if
1716                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1717               end if
1718            end if
1719
1720            if (vNam(varn) .EQ. "wpvptp") then
1721               miniwpvptp=min(data(varn,:,:))
1722               maxiwpvptp=max(data(varn,:,:))
1723               if (over .EQ. 1) then
1724                  res@xyDashPattern  = 0
1725                  plot_wpvptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1726               else
1727                  res@gsnLeftString      = vNam(varn)
1728                  res@gsnRightString     = unit(varn)
1729                  if (.not. isvar("xs")) then
1730                     if (parameter(63) .NE. "x0") then
1731                        res@trXMinF = stringtofloat(parameter(63))
1732                     else
1733                        res@trXMinF            = miniwpvptp
1734                     end if
1735                  else
1736                     res@trXMinF            = xs
1737                  end if
1738                  if (.not. isvar("xe")) then
1739                     if (parameter(65) .NE. "xdim") then
1740                        res@trXMaxF = stringtofloat(parameter(65))
1741                     else
1742                        res@trXMaxF            = maxiwpvptp
1743                     end if
1744                  else
1745                     res@trXMaxF            = xe 
1746                  end if
1747                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1748               end if
1749            end if
1750            if (vNam(varn) .EQ. "wsvpts") then
1751               miniwsvpts=min(data(varn,:,:))
1752               maxiwsvpts=max(data(varn,:,:))
1753               if (over .EQ. 1) then
1754                  res@xyDashPattern  = 1
1755                  plot_wsvpts = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1756               else
1757                  res@gsnLeftString      = vNam(varn)
1758                  res@gsnRightString     = unit(varn)
1759                  if (.not. isvar("xs")) then
1760                     if (parameter(63) .NE. "x0") then
1761                        res@trXMinF = stringtofloat(parameter(63))
1762                     else
1763                        res@trXMinF            = miniwsvpts
1764                     end if
1765                  else
1766                     res@trXMinF            = xs
1767                  end if
1768                  if (.not. isvar("xe")) then
1769                     if (parameter(65) .NE. "xdim") then
1770                        res@trXMaxF = stringtofloat(parameter(65))
1771                     else
1772                        res@trXMaxF            = maxiwsvpts
1773                     end if
1774                  else
1775                     res@trXMaxF            = xe 
1776                  end if
1777                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1778               end if
1779            end if
1780            if (vNam(varn) .EQ. "wvpt") then
1781               miniwvpt=min(data(varn,:,:))
1782               maxiwvpt=max(data(varn,:,:))
1783               if (over .EQ. 1) then
1784                  res@xyDashPattern  = 2
1785                  plot_wvpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1786               else
1787                  res@gsnLeftString      = vNam(varn)
1788                  res@gsnRightString     = unit(varn)
1789                  if (.not. isvar("xs")) then
1790                     if (parameter(63) .NE. "x0") then
1791                        res@trXMinF = stringtofloat(parameter(63))
1792                     else
1793                        res@trXMinF            = miniwvpt
1794                     end if
1795                  else
1796                     res@trXMinF            = xs
1797                  end if
1798                  if (.not. isvar("xe")) then
1799                     if (parameter(65) .NE. "xdim") then
1800                        res@trXMaxF = stringtofloat(parameter(65))
1801                     else
1802                        res@trXMaxF            = maxiwvpt
1803                     end if
1804                  else
1805                     res@trXMaxF            = xe 
1806                  end if
1807                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1808               end if
1809            end if
1810
1811            if (vNam(varn) .EQ. "wpqp") then
1812               miniwpqp=min(data(varn,:,:))
1813               maxiwpqp=max(data(varn,:,:))
1814               if (over .EQ. 1) then
1815                  res@xyDashPattern  = 0
1816                  plot_wpqp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1817               else
1818                  res@gsnLeftString      = vNam(varn)
1819                  res@gsnRightString     = unit(varn)
1820                  if (.not. isvar("xs")) then
1821                     if (parameter(63) .NE. "x0") then
1822                        res@trXMinF = stringtofloat(parameter(63))
1823                     else
1824                        res@trXMinF            = miniwpqp
1825                     end if
1826                  else
1827                     res@trXMinF            = xs
1828                  end if
1829                  if (.not. isvar("xe")) then
1830                     if (parameter(65) .NE. "xdim") then
1831                        res@trXMaxF = stringtofloat(parameter(65))
1832                     else
1833                        res@trXMaxF            = maxiwpqp
1834                     end if
1835                  else
1836                     res@trXMaxF            = xe 
1837                  end if
1838                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1839               end if
1840            end if
1841            if (vNam(varn) .EQ. "wsqs") then
1842               miniwsqs=min(data(varn,:,:))
1843               maxiwsqs=max(data(varn,:,:))
1844               if (over .EQ. 1) then
1845                  res@xyDashPattern  = 1
1846                  plot_wsqs = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1847               else
1848                  res@gsnLeftString      = vNam(varn)
1849                  res@gsnRightString     = unit(varn)
1850                  if (.not. isvar("xs")) then
1851                     if (parameter(63) .NE. "x0") then
1852                        res@trXMinF = stringtofloat(parameter(63))
1853                     else
1854                        res@trXMinF            = miniwsqs
1855                     end if
1856                  else
1857                     res@trXMinF            = xs
1858                  end if
1859                  if (.not. isvar("xe")) then
1860                     if (parameter(65) .NE. "xdim") then
1861                        res@trXMaxF = stringtofloat(parameter(65))
1862                     else
1863                        res@trXMaxF            = maxiwsqs
1864                     end if
1865                  else
1866                     res@trXMaxF            = xe 
1867                  end if
1868                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1869               end if
1870            end if
1871            if (vNam(varn) .EQ. "wq") then
1872               miniwq=min(data(varn,:,:))
1873               maxiwq=max(data(varn,:,:))
1874               if (over .EQ. 1) then
1875                  res@xyDashPattern  = 2
1876                  plot_wq = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1877               else
1878                  res@gsnLeftString      = vNam(varn)
1879                  res@gsnRightString     = unit(varn)
1880                  if (.not. isvar("xs")) then
1881                     if (parameter(63) .NE. "x0") then
1882                        res@trXMinF = stringtofloat(parameter(63))
1883                     else
1884                        res@trXMinF            = miniwq
1885                     end if
1886                  else
1887                     res@trXMinF            = xs
1888                  end if
1889                  if (.not. isvar("xe")) then
1890                     if (parameter(65) .NE. "xdim") then
1891                        res@trXMaxF = stringtofloat(parameter(65))
1892                     else
1893                        res@trXMaxF            = maxiwq
1894                     end if
1895                  else
1896                     res@trXMaxF            = xe 
1897                  end if
1898                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1899               end if
1900            end if
1901
1902            if (vNam(varn) .EQ. "wpqvp") then
1903               miniwpqvp=min(data(varn,:,:))
1904               maxiwpqvp=max(data(varn,:,:))
1905               if (over .EQ. 1) then
1906                  res@xyDashPattern  = 0
1907                  plot_wpqvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1908               else
1909                  res@gsnLeftString      = vNam(varn)
1910                  res@gsnRightString     = unit(varn)
1911                  if (.not. isvar("xs")) then
1912                     if (parameter(63) .NE. "x0") then
1913                        res@trXMinF = stringtofloat(parameter(63))
1914                     else
1915                        res@trXMinF            = miniwpqvp
1916                     end if
1917                  else
1918                     res@trXMinF            = xs
1919                  end if
1920                  if (.not. isvar("xe")) then
1921                     if (parameter(65) .NE. "xdim") then
1922                        res@trXMaxF = stringtofloat(parameter(65))
1923                     else
1924                        res@trXMaxF            = maxiwpqvp
1925                     end if
1926                  else
1927                     res@trXMaxF            = xe 
1928                  end if
1929                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1930               end if
1931            end if
1932            if (vNam(varn) .EQ. "wsqvs") then
1933               miniwsqvs=min(data(varn,:,:))
1934               maxiwsqvs=max(data(varn,:,:))
1935               if (over .EQ. 1) then
1936                  res@xyDashPattern  = 1
1937                  plot_wsqvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1938               else
1939                  res@gsnLeftString      = vNam(varn)
1940                  res@gsnRightString     = unit(varn)
1941                  if (.not. isvar("xs")) then
1942                     if (parameter(63) .NE. "x0") then
1943                        res@trXMinF = stringtofloat(parameter(63))
1944                     else
1945                        res@trXMinF            = miniwsqvs
1946                     end if
1947                  else
1948                     res@trXMinF            = xs
1949                  end if
1950                  if (.not. isvar("xe")) then
1951                     if (parameter(65) .NE. "xdim") then
1952                        res@trXMaxF = stringtofloat(parameter(65))
1953                     else
1954                        res@trXMaxF            = maxiwsqvs
1955                     end if
1956                  else
1957                     res@trXMaxF            = xe 
1958                  end if
1959                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1960               end if
1961            end if
1962            if (vNam(varn) .EQ. "wqv") then
1963               miniwqv=min(data(varn,:,:))
1964               maxiwqv=max(data(varn,:,:))
1965               if (over .EQ. 1) then
1966                  res@xyDashPattern  = 2
1967                  plot_wqv = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1968               else
1969                  res@gsnLeftString      = vNam(varn)
1970                  res@gsnRightString     = unit(varn)
1971                  if (.not. isvar("xs")) then
1972                     if (parameter(63) .NE. "x0") then
1973                        res@trXMinF = stringtofloat(parameter(63))
1974                     else
1975                        res@trXMinF            = miniwqv
1976                     end if
1977                  else
1978                     res@trXMinF            = xs
1979                  end if
1980                  if (.not. isvar("xe")) then
1981                     if (parameter(65) .NE. "xdim") then
1982                        res@trXMaxF = stringtofloat(parameter(65))
1983                     else
1984                        res@trXMaxF            = maxiwqv
1985                     end if
1986                  else
1987                     res@trXMaxF            = xe 
1988                  end if
1989                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1990               end if
1991            end if
1992
1993            if (vNam(varn) .EQ. "wpsp") then
1994               miniwpsp=min(data(varn,:,:))
1995               maxiwpsp=max(data(varn,:,:))
1996               if (over .EQ. 1) then
1997                  res@xyDashPattern  = 0
1998                  plot_wpsp = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
1999               else
2000                  res@gsnLeftString      = vNam(varn)
2001                  res@gsnRightString     = unit(varn)
2002                  if (.not. isvar("xs")) then
2003                     if (parameter(63) .NE. "x0") then
2004                        res@trXMinF = stringtofloat(parameter(63))
2005                     else
2006                        res@trXMinF            = miniwpsp
2007                     end if
2008                  else
2009                     res@trXMinF            = xs
2010                  end if
2011                  if (.not. isvar("xe")) then
2012                     if (parameter(65) .NE. "xdim") then
2013                        res@trXMaxF = stringtofloat(parameter(65))
2014                     else
2015                        res@trXMaxF            = maxiwpsp
2016                     end if
2017                  else
2018                     res@trXMaxF            = xe 
2019                  end if
2020                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2021               end if
2022            end if
2023            if (vNam(varn) .EQ. "wsss") then
2024               miniwsss=min(data(varn,:,:))
2025               maxiwsss=max(data(varn,:,:))
2026               if (over .EQ. 1) then
2027                  res@xyDashPattern  = 1
2028                  plot_wsss = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2029               else
2030                  res@gsnLeftString      = vNam(varn)
2031                  res@gsnRightString     = unit(varn)
2032                  if (.not. isvar("xs")) then
2033                     if (parameter(63) .NE. "x0") then
2034                        res@trXMinF = stringtofloat(parameter(63))
2035                     else
2036                        res@trXMinF            = miniwsss
2037                     end if
2038                  else
2039                     res@trXMinF            = xs
2040                  end if
2041                  if (.not. isvar("xe")) then
2042                     if (parameter(65) .NE. "xdim") then
2043                        res@trXMaxF = stringtofloat(parameter(65))
2044                     else
2045                        res@trXMaxF            = maxiwsss
2046                     end if
2047                  else
2048                     res@trXMaxF            = xe 
2049                  end if
2050                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2051               end if
2052            end if
2053            if (vNam(varn) .EQ. "ws") then
2054               miniws=min(data(varn,:,:))
2055               maxiws=max(data(varn,:,:))
2056               if (over .EQ. 1) then
2057                  res@xyDashPattern  = 2
2058                  plot_ws = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2059               else
2060                  res@gsnLeftString      = vNam(varn)
2061                  res@gsnRightString     = unit(varn)
2062                  if (.not. isvar("xs")) then
2063                     if (parameter(63) .NE. "x0") then
2064                        res@trXMinF = stringtofloat(parameter(63))
2065                     else
2066                        res@trXMinF            = miniws
2067                     end if
2068                  else
2069                     res@trXMinF            = xs
2070                  end if
2071                  if (.not. isvar("xe")) then
2072                     if (parameter(65) .NE. "xdim") then
2073                        res@trXMaxF = stringtofloat(parameter(65))
2074                     else
2075                        res@trXMaxF            = maxiws
2076                     end if
2077                  else
2078                     res@trXMaxF            = xe 
2079                  end if
2080                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2081               end if
2082            end if
2083
2084            if (vNam(varn) .EQ. "wpsap") then
2085               miniwpsap=min(data(varn,:,:))
2086               maxiwpsap=max(data(varn,:,:))
2087               if (over .EQ. 1) then
2088                  res@xyDashPattern  = 0
2089                  plot_wpsap = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2090               else
2091                  res@gsnLeftString      = vNam(varn)
2092                  res@gsnRightString     = unit(varn)
2093                  if (.not. isvar("xs")) then
2094                     if (parameter(63) .NE. "x0") then
2095                        res@trXMinF = stringtofloat(parameter(63))
2096                     else
2097                        res@trXMinF            = miniwpsap
2098                     end if
2099                  else
2100                     res@trXMinF            = xs
2101                  end if
2102                  if (.not. isvar("xe")) then
2103                     if (parameter(65) .NE. "xdim") then
2104                        res@trXMaxF = stringtofloat(parameter(65))
2105                     else
2106                        res@trXMaxF            = maxiwpsap
2107                     end if
2108                  else
2109                     res@trXMaxF            = xe 
2110                  end if
2111                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2112               end if
2113            end if
2114            if (vNam(varn) .EQ. "wssas") then
2115               miniwssas=min(data(varn,:,:))
2116               maxiwssas=max(data(varn,:,:))
2117               if (over .EQ. 1) then
2118                  res@xyDashPattern  = 1
2119                  plot_wssas = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2120               else
2121                  res@gsnLeftString      = vNam(varn)
2122                  res@gsnRightString     = unit(varn)
2123                  if (.not. isvar("xs")) then
2124                     if (parameter(63) .NE. "x0") then
2125                        res@trXMinF = stringtofloat(parameter(63))
2126                     else
2127                        res@trXMinF            = miniwssas
2128                     end if
2129                  else
2130                     res@trXMinF            = xs
2131                  end if
2132                  if (.not. isvar("xe")) then
2133                     if (parameter(65) .NE. "xdim") then
2134                        res@trXMaxF = stringtofloat(parameter(65))
2135                     else
2136                        res@trXMaxF            = maxiwssas
2137                     end if
2138                  else
2139                     res@trXMaxF            = xe 
2140                  end if
2141                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2142               end if
2143            end if
2144            if (vNam(varn) .EQ. "wsa") then
2145               miniwsa=min(data(varn,:,:))
2146               maxiwsa=max(data(varn,:,:))
2147               if (over .EQ. 1) then
2148                  res@xyDashPattern  = 2
2149                  plot_wsa = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2150               else
2151                  res@gsnLeftString      = vNam(varn)
2152                  res@gsnRightString     = unit(varn)
2153                  if (.not. isvar("xs")) then
2154                     if (parameter(63) .NE. "x0") then
2155                        res@trXMinF = stringtofloat(parameter(63))
2156                     else
2157                        res@trXMinF            = miniwsa
2158                     end if
2159                  else
2160                     res@trXMinF            = xs
2161                  end if
2162                  if (.not. isvar("xe")) then
2163                     if (parameter(65) .NE. "xdim") then
2164                        res@trXMaxF = stringtofloat(parameter(65))
2165                     else
2166                        res@trXMaxF            = maxiwsa
2167                     end if
2168                  else
2169                     res@trXMaxF            = xe 
2170                  end if
2171                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2172               end if
2173            end if
2174
2175            if (vNam(varn) .EQ. "us2") then
2176               minius2=min(data(varn,:,:))
2177               maxius2=max(data(varn,:,:))
2178               if (over .EQ. 1) then
2179                  res@xyDashPattern  = 0
2180                  plot_us2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2181               else
2182                  res@gsnLeftString      = vNam(varn)
2183                  res@gsnRightString     = unit(varn)
2184                  if (.not. isvar("xs")) then
2185                     if (parameter(63) .NE. "x0") then
2186                        res@trXMinF = stringtofloat(parameter(63))
2187                     else
2188                        res@trXMinF            = minius2
2189                     end if
2190                  else
2191                     res@trXMinF            = xs
2192                  end if
2193                  if (.not. isvar("xe")) then
2194                     if (parameter(65) .NE. "xdim") then
2195                        res@trXMaxF = stringtofloat(parameter(65))
2196                     else
2197                        res@trXMaxF            = maxius2
2198                     end if
2199                  else
2200                     res@trXMaxF            = xe 
2201                  end if
2202                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2203               end if
2204            end if
2205            if (vNam(varn) .EQ. "vs2") then
2206               minivs2=min(data(varn,:,:))
2207               maxivs2=max(data(varn,:,:))
2208               if (over .EQ. 1) then
2209                  res@xyDashPattern  = 1
2210                  plot_vs2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2211               else
2212                  res@gsnLeftString      = vNam(varn)
2213                  res@gsnRightString     = unit(varn)
2214                  if (.not. isvar("xs")) then
2215                     if (parameter(63) .NE. "x0") then
2216                        res@trXMinF = stringtofloat(parameter(63))
2217                     else
2218                        res@trXMinF            = minivs2
2219                     end if
2220                  else
2221                     res@trXMinF            = xs
2222                  end if
2223                  if (.not. isvar("xe")) then
2224                     if (parameter(65) .NE. "xdim") then
2225                        res@trXMaxF = stringtofloat(parameter(65))
2226                     else
2227                        res@trXMaxF            = maxivs2
2228                     end if
2229                  else
2230                     res@trXMaxF            = xe 
2231                  end if
2232                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2233               end if
2234            end if
2235            if (vNam(varn) .EQ. "ws2") then
2236               miniws2=min(data(varn,:,:))
2237               maxiws2=max(data(varn,:,:))
2238               if (over .EQ. 1) then
2239                  res@xyDashPattern  = 2
2240                  plot_ws2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2241               else
2242                  res@gsnLeftString      = vNam(varn)
2243                  res@gsnRightString     = unit(varn)
2244                  if (.not. isvar("xs")) then
2245                     if (parameter(63) .NE. "x0") then
2246                        res@trXMinF = stringtofloat(parameter(63))
2247                     else
2248                        res@trXMinF            = miniws2
2249                     end if
2250                  else
2251                     res@trXMinF            = xs
2252                  end if
2253                  if (.not. isvar("xe")) then
2254                     if (parameter(65) .NE. "xdim") then
2255                        res@trXMaxF = stringtofloat(parameter(65))
2256                     else
2257                        res@trXMaxF            = maxiws2
2258                     end if
2259                  else
2260                     res@trXMaxF            = xe 
2261                  end if
2262                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2263               end if
2264            end if
2265
2266            if (vNam(varn) .EQ. "wsususodz") then
2267               miniwsususodz=min(data(varn,:,:))
2268               maxiwsususodz=max(data(varn,:,:))
2269               if (over .EQ. 1) then
2270                  res@xyDashPattern  = 0
2271                  plot_wsususodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2272               else
2273                  res@gsnLeftString      = vNam(varn)
2274                  res@gsnRightString     = unit(varn)
2275                  if (.not. isvar("xs")) then
2276                     if (parameter(63) .NE. "x0") then
2277                        res@trXMinF = stringtofloat(parameter(63))
2278                     else
2279                        res@trXMinF            = miniwsususodz
2280                     end if
2281                  else
2282                     res@trXMinF            = xs
2283                  end if
2284                  if (.not. isvar("xe")) then
2285                     if (parameter(65) .NE. "xdim") then
2286                        res@trXMaxF = stringtofloat(parameter(65))
2287                     else
2288                        res@trXMaxF            = maxiwsususodz
2289                     end if
2290                  else
2291                     res@trXMaxF            = xe 
2292                  end if
2293                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2294               end if 
2295            end if
2296            if (vNam(varn) .EQ. "wspsodz") then
2297               miniwspsodz=min(data(varn,:,:))
2298               maxiwspsodz=max(data(varn,:,:))
2299               if (over .EQ. 1) then
2300                  res@xyDashPattern  = 1
2301                  plot_wspsodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2302               else
2303                  res@gsnLeftString      = vNam(varn)
2304                  res@gsnRightString     = unit(varn)
2305                  if (.not. isvar("xs")) then
2306                     if (parameter(63) .NE. "x0") then
2307                        res@trXMinF = stringtofloat(parameter(63))
2308                     else
2309                        res@trXMinF            = miniwspsodz
2310                     end if
2311                  else
2312                     res@trXMinF            = xs
2313                  end if
2314                  if (.not. isvar("xe")) then
2315                     if (parameter(65) .NE. "xdim") then
2316                        res@trXMaxF = stringtofloat(parameter(65))
2317                     else
2318                        res@trXMaxF            = maxiwspsodz
2319                     end if
2320                  else
2321                     res@trXMaxF            = xe 
2322                  end if
2323                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2324               end if
2325            end if
2326            if (vNam(varn) .EQ. "wpeodz") then
2327               miniwpeodz=min(data(varn,:,:))
2328               maxiwpeodz=max(data(varn,:,:))
2329               if (over .EQ. 1) then
2330                  res@xyDashPattern  = 2
2331                  plot_wpeodz = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2332               else
2333                  res@gsnLeftString      = vNam(varn)
2334                  res@gsnRightString     = unit(varn)
2335                  if (.not. isvar("xs")) then
2336                     if (parameter(63) .NE. "x0") then
2337                        res@trXMinF = stringtofloat(parameter(63))
2338                     else
2339                        res@trXMinF            = miniwpeodz
2340                     end if
2341                  else
2342                     res@trXMinF            = xs
2343                  end if
2344                  if (.not. isvar("xe")) then
2345                     if (parameter(65) .NE. "xdim") then
2346                        res@trXMaxF = stringtofloat(parameter(65))
2347                     else
2348                        res@trXMaxF            = maxiwpeodz
2349                     end if
2350                  else
2351                     res@trXMaxF            = xe 
2352                  end if
2353                  plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2354               end if
2355            end if
2356         if (over .EQ. 0) then
2357            n=n+1 
2358         end if     
2359         varn=varn+1
2360      end if     
2361      delete(temp)
2362   end do
2363
2364   if (count_var .EQ. 0) then
2365      print(" ")
2366      print("Please select a variable 'var=' or use the default value")
2367      print(" ")
2368      print("Your selection '"+var+"' does not exist on the input file")
2369      print(" ")
2370      exit
2371   end if
2372
2373   if (over .EQ. 1 ) then
2374
2375      overlay(plot_u,plot_v)
2376      overlay(plot_u,plot_w)
2377      u=0
2378      overlay(plot_pt,plot_vpt)
2379      overlay(plot_pt,plot_lpt)
2380      pt=0
2381      overlay(plot_q,plot_qv)
2382      overlay(plot_q,plot_ql)
2383      q=0
2384      overlay(plot_e,plot_es)
2385      e=0
2386      overlay(plot_km,plot_kh)
2387      km=0
2388      overlay(plot_wpup,plot_wsus)
2389      overlay(plot_wpup,plot_wu)
2390      wpup=0
2391      overlay(plot_wpvp,plot_wsvs)
2392      overlay(plot_wpvp,plot_wv)
2393      wpvp=0
2394      overlay(plot_wpptp,plot_wspts)
2395      overlay(plot_wpptp,plot_wpt)
2396      wpptp=0
2397      overlay(plot_wsptsBC,plot_wptBC)
2398      wsptsBC=0
2399      overlay(plot_wpvptp,plot_wsvpts)
2400      overlay(plot_wpvptp,plot_wvpt)
2401      wpvptp=0
2402      overlay(plot_wpqp,plot_wsqs)
2403      overlay(plot_wpqp,plot_wq)
2404      wpqp=0
2405      overlay(plot_wpqvp,plot_wsqvs)
2406      overlay(plot_wpqvp,plot_wqv)
2407      wpqvp=0
2408      overlay(plot_wpsp,plot_wsss)
2409      overlay(plot_wpsp,plot_ws)
2410      wpsp=0
2411      overlay(plot_wpsap,plot_wssas)
2412      overlay(plot_wpsap,plot_wsa)
2413      wpsap=0
2414      overlay(plot_us2,plot_vs2)
2415      overlay(plot_us2,plot_ws2)
2416      us2=0
2417      overlay(plot_wsususodz,plot_wspsodz)
2418      overlay(plot_wsususodz,plot_wpeodz)
2419      wsususodz=0
2420
2421   end if
2422
2423   if (over .EQ. 1) then
2424   
2425      do varn = 0,dim-1   
2426     
2427         if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then
2428            check = False
2429         else
2430            if (.not. isvar("var")) then
2431               check = True
2432               if (parameter(21) .NE. "variables") then
2433                  var=parameter(21)
2434                  check = isStrSubset( var,","+vNam(varn)+"," )
2435               end if
2436            else         
2437               check = isStrSubset( var,","+vNam(varn)+"," )
2438            end if
2439         end if 
2440
2441         if (check)
2442           
2443            z = f->$vNam(varn+1)$
2444           
2445            res@gsnLeftString      = vNam(varn)
2446            res@gsnRightString     = unit(varn)
2447            res@trYMinF            = min_z
2448            res@trYMaxF            = max_z 
2449            if (.not. isvar("xs")) then
2450               if (parameter(63) .NE. "x0") then
2451                  res@trXMinF = stringtofloat(parameter(63))
2452               else
2453                  res@trXMinF = min(data(varn,:,:))
2454               end if
2455            else
2456               res@trXMinF = xs
2457            end if
2458            if (.not. isvar("xe")) then
2459               if (parameter(65) .NE. "xdim") then
2460                  res@trXMaxF = stringtofloat(parameter(65))
2461               else
2462                  res@trXMaxF = max(data(varn,:,:))
2463            end if
2464            else
2465               res@trXMaxF = xe 
2466            end if
2467            plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 
2468           
2469            if (vNam(varn) .EQ. "u" .OR. vNam(varn) .EQ. "v" .OR. vNam(varn) .EQ. "w") then
2470               if (u .EQ. 0) then
2471                  res@gsnLeftString      = "u, v and w"
2472                  res@gsnRightString     = unit(varn)
2473                  if (.not. isvar("xs")) then
2474                     if (parameter(63) .NE. "x0") then
2475                        res@trXMinF = stringtofloat(parameter(63))
2476                     else
2477                        res@trXMinF = min((/miniu,miniv,miniw/))
2478                     end if
2479                  else
2480                     res@trXMinF = xs
2481                  end if
2482                  if (.not. isvar("xe")) then
2483                     if (parameter(65) .NE. "xdim") then
2484                        res@trXMaxF = stringtofloat(parameter(65))
2485                     else
2486                        res@trXMaxF = max((/maxiu,maxiv,maxiw/))
2487                     end if
2488                  else
2489                     res@trXMaxF = xe 
2490                  end if
2491                  if (vNam(varn) .EQ. "w") then
2492                     plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2493                  else
2494                     plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2495                  end if
2496
2497                  ; ***************************************************
2498                  ; legend for combined plot
2499                  ; ***************************************************
2500     
2501                  lgres                    = True
2502                  lgMonoDashIndex          = False
2503                  lgres@lgLabelFont        = "helvetica"   
2504                  lgres@lgLabelFontHeightF = .1           
2505                  lgres@vpWidthF           = 0.12           
2506                  lgres@vpHeightF          = 0.1         
2507                  lgres@lgDashIndexes      = (/0,1,2/)
2508                  lbid = gsn_create_legend(wks,3,(/"u","v","w"/),lgres)       
2509
2510                  amres = True
2511                  amres@amParallelPosF   = 0.65             
2512                  amres@amOrthogonalPosF = -0.2           
2513                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2514                  overlay(plot(n),plot_u)
2515                  u=1
2516               else
2517                  varn=varn+1
2518                  continue
2519               end if       
2520            end if 
2521     
2522            if (vNam(varn) .EQ. "pt" .OR. vNam(varn) .EQ. "vpt" .OR. vNam(varn) .EQ. "lpt") then
2523               if (pt .EQ. 0) then
2524                  res@gsnLeftString      = "pt, vpt and lpt"
2525                  res@gsnRightString     = unit(varn)
2526                  if (.not. isvar("xs")) then
2527                     if (parameter(63) .NE. "x0") then
2528                        res@trXMinF = stringtofloat(parameter(63))
2529                     else
2530                        res@trXMinF = min((/minipt,minivpt,minilpt/))
2531                     end if
2532                  else
2533                     res@trXMinF = xs
2534                  end if
2535                  if (.not. isvar("xe")) then
2536                     if (parameter(65) .NE. "xdim") then
2537                        res@trXMaxF = stringtofloat(parameter(65))
2538                     else
2539                        res@trXMaxF = max((/maxipt,maxivpt,maxilpt/))
2540                     end if
2541                  else
2542                     res@trXMaxF = xe 
2543                  end if
2544
2545                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2546
2547                  ; ***************************************************
2548                  ; legend for combined plot
2549                  ; ***************************************************
2550     
2551                  lgres                    = True
2552                  lgMonoDashIndex          = False
2553                  lgres@lgLabelFont        = "helvetica"   
2554                  lgres@lgLabelFontHeightF = .1           
2555                  lgres@vpWidthF           = 0.12           
2556                  lgres@vpHeightF          = 0.1         
2557                  lgres@lgDashIndexes      = (/0,1,2/)
2558                  lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres)       
2559
2560                  amres = True
2561                  amres@amParallelPosF   = 0.65             
2562                  amres@amOrthogonalPosF = -0.2           
2563                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2564                  overlay(plot(n),plot_pt)
2565                  pt=1
2566               else
2567                  varn=varn+1
2568                  continue       
2569               end if
2570            end if           
2571            if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then
2572               if (q .EQ. 0) then
2573                  res@gsnLeftString      = "q, qv and ql"
2574                  res@gsnRightString     = unit(varn)
2575                  if (.not. isvar("xs")) then
2576                     if (parameter(63) .NE. "x0") then
2577                        res@trXMinF = stringtofloat(parameter(63))
2578                     else
2579                        res@trXMinF = min((/miniq,miniqv,miniql/))
2580                     end if
2581                  else
2582                     res@trXMinF = xs
2583                  end if
2584                  if (.not. isvar("xe")) then
2585                     if (parameter(65) .NE. "xdim") then
2586                        res@trXMaxF = stringtofloat(parameter(65))
2587                     else
2588                        res@trXMaxF = max((/maxiq,maxiqv,maxiql/))
2589                     end if
2590                  else
2591                     res@trXMaxF = xe 
2592                  end if
2593
2594                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2595
2596                  ; ***************************************************
2597                  ; legend for combined plot
2598                  ; ***************************************************
2599     
2600                  lgres                    = True
2601                  lgMonoDashIndex          = False
2602                  lgres@lgLabelFont        = "helvetica"   
2603                  lgres@lgLabelFontHeightF = .1           
2604                  lgres@vpWidthF           = 0.12           
2605                  lgres@vpHeightF          = 0.1         
2606                  lgres@lgDashIndexes      = (/0,1,2/)
2607                  lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres)       
2608
2609                  amres = True
2610                  amres@amParallelPosF   = 0.65             
2611                  amres@amOrthogonalPosF = -0.2           
2612                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2613                  overlay(plot(n),plot_q)
2614                  q=1
2615               else
2616                  varn=varn+1
2617                  continue   
2618               end if
2619            end if   
2620           
2621            if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then
2622               if (e .EQ. 0) then
2623                  res@gsnLeftString      = "e and es"
2624                  res@gsnRightString     = unit(varn)
2625                  if (.not. isvar("xs")) then
2626                     if (parameter(63) .NE. "x0") then
2627                        res@trXMinF = stringtofloat(parameter(63))
2628                     else
2629                        res@trXMinF = min((/minie,minies/))
2630                     end if
2631                  else
2632                     res@trXMinF = xs
2633                  end if
2634                  if (.not. isvar("xe")) then
2635                     if (parameter(65) .NE. "xdim") then
2636                        res@trXMaxF = stringtofloat(parameter(65))
2637                     else
2638                        res@trXMaxF = max((/maxie,maxies/))
2639                     end if
2640                  else
2641                     res@trXMaxF = xe 
2642                  end if
2643
2644                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2645
2646                  ; ***************************************************
2647                  ; legend for combined plot
2648                  ; ***************************************************
2649     
2650                  lgres                    = True
2651                  lgMonoDashIndex          = False
2652                  lgres@lgLabelFont        = "helvetica"   
2653                  lgres@lgLabelFontHeightF = .1           
2654                  lgres@vpWidthF           = 0.12           
2655                  lgres@vpHeightF          = 0.1         
2656                  lgres@lgDashIndexes      = (/0,1,2/)
2657                  lbid = gsn_create_legend(wks,3,(/"e","es"/),lgres)       
2658
2659                  amres = True
2660                  amres@amParallelPosF   = 0.65             
2661                  amres@amOrthogonalPosF = -0.2           
2662                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2663                  overlay(plot(n),plot_e)
2664                  e=1
2665               else
2666                  varn=varn+1
2667                  continue   
2668               end if
2669            end if           
2670            if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then
2671               if (km .EQ. 0) then
2672                  res@gsnLeftString      = "km and kh"
2673                  res@gsnRightString     = unit(varn)
2674                  if (.not. isvar("xs")) then
2675                     if (parameter(63) .NE. "x0") then
2676                        res@trXMinF = stringtofloat(parameter(63))
2677                     else
2678                        res@trXMinF = min((/minikm,minikh/))
2679                     end if
2680                  else
2681                     res@trXMinF = xs
2682                  end if
2683                  if (.not. isvar("xe")) then
2684                     if (parameter(65) .NE. "xdim") then
2685                        res@trXMaxF = stringtofloat(parameter(65))
2686                     else
2687                        res@trXMaxF = max((/maxikm,maxikh/))
2688                     end if
2689                  else
2690                     res@trXMaxF = xe 
2691                  end if
2692
2693                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_u,res)
2694
2695                  ; ***************************************************
2696                  ; legend for combined plot
2697                  ; ***************************************************
2698     
2699                  lgres                    = True
2700                  lgMonoDashIndex          = False
2701                  lgres@lgLabelFont        = "helvetica"   
2702                  lgres@lgLabelFontHeightF = .1           
2703                  lgres@vpWidthF           = 0.12           
2704                  lgres@vpHeightF          = 0.1         
2705                  lgres@lgDashIndexes      = (/0,1,2/)
2706                  lbid = gsn_create_legend(wks,3,(/"km","kh"/),lgres)       
2707
2708                  amres = True
2709                  amres@amParallelPosF   = 0.65             
2710                  amres@amOrthogonalPosF = -0.2           
2711                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2712                  overlay(plot(n),plot_km)
2713                  km=1
2714               else
2715                  varn=varn+1
2716                  continue   
2717               end if
2718            end if           
2719           
2720            if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then
2721               if (wpup .EQ. 0) then
2722                  res@gsnLeftString      = "wpup, wsus and wu"
2723                  res@gsnRightString     = unit(varn)
2724                  if (.not. isvar("xs")) then
2725                     if (parameter(63) .NE. "x0") then
2726                        res@trXMinF = stringtofloat(parameter(63))
2727                     else
2728                        res@trXMinF = min((/miniwpup,miniwsus,miniwu/))
2729                     end if
2730                  else
2731                     res@trXMinF = xs
2732                  end if
2733                  if (.not. isvar("xe")) then
2734                     if (parameter(65) .NE. "xdim") then
2735                        res@trXMaxF = stringtofloat(parameter(65))
2736                     else
2737                        res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/))
2738                     end if
2739                  else
2740                     res@trXMaxF = xe 
2741                  end if
2742
2743                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2744
2745                  ; ***************************************************
2746                  ; legend for combined plot
2747                  ; ***************************************************
2748     
2749                  lgres                    = True
2750                  lgMonoDashIndex          = False
2751                  lgres@lgLabelFont        = "helvetica"   
2752                  lgres@lgLabelFontHeightF = .1           
2753                  lgres@vpWidthF           = 0.12           
2754                  lgres@vpHeightF          = 0.1         
2755                  lgres@lgDashIndexes      = (/0,1,2/)
2756                  lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres)       
2757
2758                  amres = True
2759                  amres@amParallelPosF   = 0.65             
2760                  amres@amOrthogonalPosF = -0.2           
2761                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2762                  overlay(plot(n),plot_wpup)
2763                  wpup=1
2764               else
2765                  varn=varn+1
2766                  continue   
2767               end if
2768            end if
2769            if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then
2770               if (wpvp .EQ. 0) then
2771                  res@gsnLeftString      = "wpvp, wsus and wv"
2772                  res@gsnRightString     = unit(varn)
2773                  if (.not. isvar("xs")) then
2774                     if (parameter(63) .NE. "x0") then
2775                        res@trXMinF = stringtofloat(parameter(63))
2776                     else
2777                        res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/))
2778                     end if
2779                  else
2780                     res@trXMinF = xs
2781                  end if
2782                  if (.not. isvar("xe")) then
2783                     if (parameter(65) .NE. "xdim") then
2784                        res@trXMaxF = stringtofloat(parameter(65))
2785                     else
2786                        res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/))
2787                     end if
2788                  else
2789                     res@trXMaxF = xe 
2790                  end if
2791
2792                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2793
2794                  ; ***************************************************
2795                  ; legend for combined plot
2796                  ; ***************************************************
2797     
2798                  lgres                    = True
2799                  lgMonoDashIndex          = False
2800                  lgres@lgLabelFont        = "helvetica"   
2801                  lgres@lgLabelFontHeightF = .1           
2802                  lgres@vpWidthF           = 0.12           
2803                  lgres@vpHeightF          = 0.1         
2804                  lgres@lgDashIndexes      = (/0,1,2/)
2805                  lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres)       
2806
2807                  amres = True
2808                  amres@amParallelPosF   = 0.65             
2809                  amres@amOrthogonalPosF = -0.2           
2810                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2811                  overlay(plot(n),plot_wpvp)
2812                  wpup=0
2813               else
2814                  varn=varn+1
2815                  continue   
2816               end if
2817            end if
2818            if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then
2819               if (wpptp .EQ. 0) then
2820                  res@gsnLeftString      = "wpptp, wspts and wv"
2821                  res@gsnRightString     = unit(varn)
2822                  if (.not. isvar("xs")) then
2823                     if (parameter(63) .NE. "x0") then
2824                        res@trXMinF = stringtofloat(parameter(63))
2825                     else
2826                        res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/))
2827                     end if
2828                  else
2829                     res@trXMinF = xs
2830                  end if
2831                  if (.not. isvar("xe")) then
2832                     if (parameter(65) .NE. "xdim") then
2833                        res@trXMaxF = stringtofloat(parameter(65))
2834                     else
2835                        res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/))
2836                     end if
2837                  else
2838                     res@trXMaxF = xe 
2839                  end if
2840
2841                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2842
2843                  ; ***************************************************
2844                  ; legend for combined plot
2845                  ; ***************************************************
2846     
2847                  lgres                    = True
2848                  lgMonoDashIndex          = False
2849                  lgres@lgLabelFont        = "helvetica"   
2850                  lgres@lgLabelFontHeightF = .1           
2851                  lgres@vpWidthF           = 0.12           
2852                  lgres@vpHeightF          = 0.1         
2853                  lgres@lgDashIndexes      = (/0,1,2/)
2854                  lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres)       
2855
2856                  amres = True
2857                  amres@amParallelPosF   = 0.65             
2858                  amres@amOrthogonalPosF = -0.2           
2859                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2860                  overlay(plot(n),plot_wpptp)
2861                  wpptp=1
2862               else
2863                  varn=varn+1
2864                  continue   
2865               end if
2866            end if
2867            if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then
2868               if (wsptsBC .EQ. 0) then
2869                  res@gsnLeftString      = "wsptsBC and wptBC"
2870                  res@gsnRightString     = unit(varn)
2871                  if (.not. isvar("xs")) then
2872                     if (parameter(63) .NE. "x0") then
2873                        res@trXMinF = stringtofloat(parameter(63))
2874                     else
2875                        res@trXMinF = min((/miniwsptsBC,miniwptBC/))
2876                     end if
2877                  else
2878                     res@trXMinF = xs
2879                  end if
2880                  if (.not. isvar("xe")) then
2881                     if (parameter(65) .NE. "xdim") then
2882                        res@trXMaxF = stringtofloat(parameter(65))
2883                     else
2884                        res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/))
2885                     end if
2886                  else
2887                     res@trXMaxF = xe 
2888                  end if
2889
2890                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2891
2892                  ; ***************************************************
2893                  ; legend for combined plot
2894                  ; ***************************************************
2895     
2896                  lgres                    = True
2897                  lgMonoDashIndex          = False
2898                  lgres@lgLabelFont        = "helvetica"   
2899                  lgres@lgLabelFontHeightF = .1           
2900                  lgres@vpWidthF           = 0.12           
2901                  lgres@vpHeightF          = 0.1         
2902                  lgres@lgDashIndexes      = (/0,1,2/)
2903                  lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres)       
2904
2905                  amres = True
2906                  amres@amParallelPosF   = 0.65             
2907                  amres@amOrthogonalPosF = -0.2           
2908                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2909                  overlay(plot(n),plot_wsptsBC)
2910                  wsptsBC=1
2911               else
2912                  varn=varn+1
2913                  continue   
2914               end if 
2915            end if             
2916            if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then
2917               if (wpvptp .EQ. 0) then
2918                  res@gsnLeftString      = "wpvptp, wsvpts and wv"
2919                  res@gsnRightString     = unit(varn)
2920                  if (.not. isvar("xs")) then
2921                     if (parameter(63) .NE. "x0") then
2922                        res@trXMinF = stringtofloat(parameter(63))
2923                     else
2924                        res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/))
2925                     end if
2926                  else
2927                     res@trXMinF = xs
2928                  end if
2929                  if (.not. isvar("xe")) then
2930                     if (parameter(65) .NE. "xdim") then
2931                        res@trXMaxF = stringtofloat(parameter(65))
2932                     else
2933                        res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/))
2934                     end if
2935                  else
2936                     res@trXMaxF = xe 
2937                  end if
2938
2939                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2940
2941                  ; ***************************************************
2942                  ; legend for combined plot
2943                  ; ***************************************************
2944     
2945                  lgres                    = True
2946                  lgMonoDashIndex          = False
2947                  lgres@lgLabelFont        = "helvetica"   
2948                  lgres@lgLabelFontHeightF = .1           
2949                  lgres@vpWidthF           = 0.12           
2950                  lgres@vpHeightF          = 0.1         
2951                  lgres@lgDashIndexes      = (/0,1,2/)
2952                  lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres)       
2953
2954                  amres = True
2955                  amres@amParallelPosF   = 0.65             
2956                  amres@amOrthogonalPosF = -0.2           
2957                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
2958                  overlay(plot(n),plot_wpvptp)
2959                  wpvptp=1
2960               else
2961                  varn=varn+1
2962                  continue   
2963               end if
2964            end if
2965            if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then
2966               if (wpqp .EQ. 0) then
2967                  res@gsnLeftString      = "wpqp, wsqs and wq"
2968                  res@gsnRightString     = unit(varn)
2969                  if (.not. isvar("xs")) then
2970                     if (parameter(63) .NE. "x0") then
2971                        res@trXMinF = stringtofloat(parameter(63))
2972                     else
2973                        res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/))
2974                     end if
2975                  else
2976                     res@trXMinF = xs
2977                  end if
2978                  if (.not. isvar("xe")) then
2979                     if (parameter(65) .NE. "xdim") then
2980                        res@trXMaxF = stringtofloat(parameter(65))
2981                     else
2982                        res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/))
2983                     end if
2984                  else
2985                     res@trXMaxF = xe 
2986                  end if
2987
2988                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
2989
2990                  ; ***************************************************
2991                  ; legend for combined plot
2992                  ; ***************************************************
2993     
2994                  lgres                    = True
2995                  lgMonoDashIndex          = False
2996                  lgres@lgLabelFont        = "helvetica"   
2997                  lgres@lgLabelFontHeightF = .1           
2998                  lgres@vpWidthF           = 0.12           
2999                  lgres@vpHeightF          = 0.1         
3000                  lgres@lgDashIndexes      = (/0,1,2/)
3001                  lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres)       
3002
3003                  amres = True
3004                  amres@amParallelPosF   = 0.65             
3005                  amres@amOrthogonalPosF = -0.2           
3006                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
3007                  overlay(plot(n),plot_wpqp)
3008                  wpqp=1
3009               else
3010                  varn=varn+1
3011                  continue   
3012               end if
3013            end if
3014            if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then
3015               if (wpqvp .EQ. 0) then
3016                  res@gsnLeftString      = "wpqvp, wsqvs and wqv"
3017                  res@gsnRightString     = unit(varn)
3018                  if (.not. isvar("xs")) then
3019                     if (parameter(63) .NE. "x0") then
3020                        res@trXMinF = stringtofloat(parameter(63))
3021                     else
3022                        res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/))
3023                     end if
3024                  else
3025                     res@trXMinF = xs
3026                  end if
3027                  if (.not. isvar("xe")) then
3028                     if (parameter(65) .NE. "xdim") then
3029                        res@trXMaxF = stringtofloat(parameter(65))
3030                     else
3031                        res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/))
3032                     end if
3033                  else
3034                     res@trXMaxF = xe 
3035                  end if
3036
3037                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3038
3039                  ; ***************************************************
3040                  ; legend for combined plot
3041                  ; ***************************************************
3042     
3043                  lgres                    = True
3044                  lgMonoDashIndex          = False
3045                  lgres@lgLabelFont        = "helvetica"   
3046                  lgres@lgLabelFontHeightF = .1           
3047                  lgres@vpWidthF           = 0.12           
3048                  lgres@vpHeightF          = 0.1         
3049                  lgres@lgDashIndexes      = (/0,1,2/)
3050                  lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres)       
3051
3052                  amres = True
3053                  amres@amParallelPosF   = 0.65             
3054                  amres@amOrthogonalPosF = -0.2           
3055                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
3056                  overlay(plot(n),plot_wpqvp)
3057                  wpqvp=1
3058               else
3059                  varn=varn+1
3060                  continue   
3061               end if
3062            end if
3063            if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then
3064               if (wpsp .EQ. 0) then
3065                  res@gsnLeftString      = "wpsp, wsss and ws"
3066                  res@gsnRightString     = unit(varn)
3067                  if (.not. isvar("xs")) then
3068                     if (parameter(63) .NE. "x0") then
3069                        res@trXMinF = stringtofloat(parameter(63))
3070                     else
3071                        res@trXMinF = min((/miniwpsp,miniwsss,miniws/))
3072                     end if
3073                  else
3074                     res@trXMinF = xs
3075                  end if
3076                  if (.not. isvar("xe")) then
3077                     if (parameter(65) .NE. "xdim") then
3078                        res@trXMaxF = stringtofloat(parameter(65))
3079                     else
3080                        res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/))
3081                     end if
3082                  else
3083                     res@trXMaxF = xe 
3084                  end if
3085
3086                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3087
3088                  ; ***************************************************
3089                  ; legend for combined plot
3090                  ; ***************************************************
3091     
3092                  lgres                    = True
3093                  lgMonoDashIndex          = False
3094                  lgres@lgLabelFont        = "helvetica"   
3095                  lgres@lgLabelFontHeightF = .1           
3096                  lgres@vpWidthF           = 0.12           
3097                  lgres@vpHeightF          = 0.1         
3098                  lgres@lgDashIndexes      = (/0,1,2/)
3099                  lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres)       
3100
3101                  amres = True
3102                  amres@amParallelPosF   = 0.65             
3103                  amres@amOrthogonalPosF = -0.2           
3104                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
3105                  overlay(plot(n),plot_wpsp)
3106                  wpsp=1
3107               else
3108                  varn=varn+1
3109                  continue   
3110               end if
3111            end if
3112            if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then
3113               if (wpsap .EQ. 0) then
3114                  res@gsnLeftString      = "wpsap, wssas and wsa"
3115                  res@gsnRightString     = unit(varn)
3116                  if (.not. isvar("xs")) then
3117                     if (parameter(63) .NE. "x0") then
3118                        res@trXMinF = stringtofloat(parameter(63))
3119                     else
3120                        res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/))
3121                     end if
3122                  else
3123                     res@trXMinF = xs
3124                  end if
3125                  if (.not. isvar("xe")) then
3126                     if (parameter(65) .NE. "xdim") then
3127                        res@trXMaxF = stringtofloat(parameter(65))
3128                     else
3129                        res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/))
3130                     end if
3131                  else
3132                     res@trXMaxF = xe 
3133                  end if
3134
3135                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3136
3137                  ; ***************************************************
3138                  ; legend for combined plot
3139                  ; ***************************************************
3140     
3141                  lgres                    = True
3142                  lgMonoDashIndex          = False
3143                  lgres@lgLabelFont        = "helvetica"   
3144                  lgres@lgLabelFontHeightF = .1           
3145                  lgres@vpWidthF           = 0.12           
3146                  lgres@vpHeightF          = 0.1         
3147                  lgres@lgDashIndexes      = (/0,1,2/)
3148                  lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres)       
3149
3150                  amres = True
3151                  amres@amParallelPosF   = 0.65             
3152                  amres@amOrthogonalPosF = -0.2           
3153                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
3154                  overlay(plot(n),plot_wpsap)
3155                  wpsap=1
3156               else
3157                  varn=varn+1
3158                  continue   
3159               end if
3160            end if
3161         
3162            if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then
3163               if (us2 .EQ. 0) then
3164                  res@gsnLeftString      = "us2, vs2 and ws2"
3165                  res@gsnRightString     = unit(varn)
3166                  if (.not. isvar("xs")) then
3167                     if (parameter(63) .NE. "x0") then
3168                        res@trXMinF = stringtofloat(parameter(63))
3169                     else
3170                        res@trXMinF = min((/minius2,minivs2,miniws2/))
3171                     end if
3172                  else
3173                     res@trXMinF = xs
3174                  end if
3175                  if (.not. isvar("xe")) then
3176                     if (parameter(65) .NE. "xdim") then
3177                        res@trXMaxF = stringtofloat(parameter(65))
3178                     else
3179                        res@trXMaxF = max((/maxius2,maxivs2,maxiws2/))
3180                     end if
3181                  else
3182                     res@trXMaxF = xe 
3183                  end if
3184
3185                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3186
3187                  ; ***************************************************
3188                  ; legend for combined plot
3189                  ; ***************************************************
3190     
3191                  lgres                    = True
3192                  lgMonoDashIndex          = False
3193                  lgres@lgLabelFont        = "helvetica"   
3194                  lgres@lgLabelFontHeightF = .1           
3195                  lgres@vpWidthF           = 0.12           
3196                  lgres@vpHeightF          = 0.1         
3197                  lgres@lgDashIndexes      = (/0,1,2/)
3198                  lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres)       
3199
3200                  amres = True
3201                  amres@amParallelPosF   = 0.65             
3202                  amres@amOrthogonalPosF = -0.2           
3203                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
3204                  overlay(plot(n),plot_us2)
3205                  us2=1
3206               else
3207                  varn=varn+1
3208                  continue   
3209               end if
3210            end if
3211           
3212            if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then
3213               if (wsususodz .EQ. 0) then
3214                  res@gsnLeftString      = "wsususodz, wspsodz and ws2"
3215                  res@gsnRightString     = unit(varn)
3216                  if (.not. isvar("xs")) then
3217                     if (parameter(63) .NE. "x0") then
3218                        res@trXMinF = stringtofloat(parameter(63))
3219                     else
3220                        res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/))
3221                     end if
3222                  else
3223                     res@trXMinF = xs
3224                  end if
3225                  if (.not. isvar("xe")) then
3226                     if (parameter(65) .NE. "xdim") then
3227                        res@trXMaxF = stringtofloat(parameter(65))
3228                     else
3229                        res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/))
3230                     end if
3231                  else
3232                     res@trXMaxF = xe 
3233                  end if
3234
3235                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z_w,res)
3236
3237                  ; ***************************************************
3238                  ; legend for combined plot
3239                  ; ***************************************************
3240     
3241                  lgres                    = True
3242                  lgMonoDashIndex          = False
3243                  lgres@lgLabelFont        = "helvetica"   
3244                  lgres@lgLabelFontHeightF = .1           
3245                  lgres@vpWidthF           = 0.12           
3246                  lgres@vpHeightF          = 0.1         
3247                  lgres@lgDashIndexes      = (/0,1,2/)
3248                  lbid = gsn_create_legend(wks,3,(/"wsususodz","wspsodz","wpeodz"/),lgres)       
3249
3250                  amres = True
3251                  amres@amParallelPosF   = 0.65             
3252                  amres@amOrthogonalPosF = -0.2           
3253                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
3254                  overlay(plot(n),plot_wsususodz)
3255                  wsususodz=1
3256               else
3257                  varn=varn+1
3258                  continue   
3259               end if
3260            end if     
3261            n=n+1   
3262            varn=varn+1
3263         end if
3264      end do
3265   end if     
3266
3267   if (combine .EQ. 1) then
3268     
3269      check=True     
3270      n_o=0
3271      do varn = 0,dim-1
3272 
3273         if ( isStrSubset( vNam(varn), "NORM") .or. isStrSubset( vNam(varn), "time") )
3274            check = False
3275         else
3276            if (.not. isvar("var")) then
3277               check = True
3278               if (parameter(21) .NE. "variables") then
3279                  var=parameter(21)
3280                  check = isStrSubset( var,","+vNam(varn)+"," )
3281               end if
3282            else         
3283               check = isStrSubset( var,","+vNam(varn)+"," )
3284            end if
3285         end if
3286     
3287         if(check) then
3288
3289            z = f->$vNam(varn+1)$
3290            com=isStrSubset( c_var,","+vNam(varn)+"," )
3291     
3292            if (com)
3293                             
3294               if (n_o .EQ. 1) then
3295                  res@xyDashPattern  = 1
3296                                 
3297               else           
3298                  if (n_o .EQ. 2) then
3299                     res@xyDashPattern  = 2
3300                  else
3301                     res@xyDashPattern  = 0
3302                     res@gsnLeftString  = c_var
3303                     res@gsnRightString = unit(varn)
3304                     if (.not. isvar("xs")) then
3305                        if (parameter(63) .NE. "x0") then
3306                           res@trXMinF = stringtofloat(parameter(63))
3307                        else
3308                           res@trXMinF = min(mini)
3309                        end if
3310                     else
3311                        res@trXMinF = xs
3312                     end if
3313                     if (.not. isvar("xe")) then
3314                        if (parameter(65) .NE. "xdim") then
3315                           res@trXMaxF = stringtofloat(parameter(65))
3316                        else
3317                           res@trXMaxF = max(maxi)
3318                        end if
3319                     else
3320                        res@trXMaxF = xe 
3321                     end if
3322                  end if
3323               end if
3324               label(n_o)=vNam(varn)
3325               color_o(n_o)=237
3326               plot_o(n_o)=gsn_csm_xy(wks,data(varn,:,:),z,res)
3327               n_o=n_o+1
3328            end if         
3329            varn=varn+1   
3330         end if
3331      end do
3332   
3333      if(number_comb .EQ. 2)then
3334         overlay(plot_o(0),plot_o(1))
3335      end if
3336      if(number_comb .EQ. 3)then
3337         overlay(plot_o(0),plot_o(1))
3338         overlay(plot_o(0),plot_o(2))
3339      end if
3340
3341      ; ***************************************************
3342      ; legend for combined plot
3343      ; ***************************************************
3344     
3345      lgres                    = True
3346      lgMonoDashIndex          = False
3347      lgres@lgDashIndexes      = (/0,1,2/)
3348      lgres@lgLabelFont        = "helvetica"   
3349      lgres@lgLabelFontHeightF = .1           
3350      lgres@vpWidthF           = 0.12           
3351      lgres@vpHeightF          = 0.1           
3352 
3353      lbid = gsn_create_legend(wks,number_comb,label,lgres)       
3354
3355      amres = True
3356      amres@amParallelPosF   = 0.65                 
3357      amres@amOrthogonalPosF = -0.2           
3358      annoid1 = gsn_add_annotation(plot_o(0),lbid,amres)
3359   
3360      plot(0) = plot_o(0)
3361
3362   end if
3363
3364   ; ***************************************************
3365   ; merge plots onto one page
3366   ; ***************************************************
3367
3368   do m=0,n-1
3369      plot_(m)=plot(n-1-m)
3370   end do
3371
3372   if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
3373      gsn_panel(wks,plot_(0:n-1),(/n,1/),resP)
3374   else   
3375      do i = 0,n-1, no_lines*no_columns
3376         if( (i+no_lines*no_columns) .gt. (n-1)) then
3377            gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP)
3378         else
3379           gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
3380         end if
3381      end do
3382   end if
3383
3384   print(" ")
3385   print("Output to: " + file_out +"."+ format_out)
3386   print(" ")
3387
3388end
Note: See TracBrowser for help on using the repository browser.