Changeset 250 for palm/trunk/SCRIPTS/NCL


Ignore:
Timestamp:
Mar 2, 2009 12:18:17 PM (13 years ago)
Author:
letzel
Message:
  • NCL scripts in trunk/SCRIPTS/NCL updated
Location:
palm/trunk/SCRIPTS/NCL
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/profiles.ncl

    r236 r250  
    441441            else
    442442               if(vNam(varn) .EQ. "vpt" .OR. isStrSubset(vNam(varn), "vpt_"))then
    443                   z_u = f->$vNam(varn+1)$
     443                  z_u = f_att->$vNam(varn+1)$
    444444                  break
    445445               else
     
    812812   
    813813   delta_t=t_all(nt-1)/nt
    814    
    815    
    816    print(nt)
    817    print(t_all)
    818    
    819    
    820    ;exit
    821814
    822815   ; ****************************************************       
     
    946939   res@pmLegendOrthogonalPosF  = -1.0
    947940   res@pmLegendWidthF          = 0.12
    948    res@pmLegendHeightF         = 0.04*(end_time_step-start_time_step+1)
     941   res@pmLegendHeightF         = 0.05*(end_time_step-start_time_step+1)
    949942   res@lgLabelFontHeightF     = font_size
    950943   res@lgTitleString      = "Time [h]"
     
    957950   res@tiXAxisString      = " "
    958951   if ( black .eq. 0 ) then 
    959       res@xyLineColors = ispan(2,237,235/np)
     952      res@xyLineColors = -(ispan(-237,-2,235/np))
    960953   end if
    961954   if (norm_z .EQ. 1)then
     
    972965      res@xyMonoDashPattern       = True
    973966   else
     967      res@xyMonoDashPattern       = False
    974968      if (no_files .GT. 1)
    975969         res@xyMonoDashPattern       = True 
     
    979973      end if       
    980974   end if
    981    
     975
    982976   res@tmXBMinorPerMajor = 4
    983977   res@tmYLMinorPerMajor = 4
     
    14041398   n_o=0
    14051399   count_var=0
     1400
    14061401   res@xyDashPattern = 1*nof
    1407 
     1402   
    14081403   do varn = 0,dim-1
    14091404
     
    15791574            end if         
    15801575            plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 
    1581          end if       
    1582          
     1576         end if 
     1577
    15831578            if (vNam(varn) .EQ. "u") then
    15841579               miniu=min(data(varn,:,:))
     
    16081603               maxiv=max(data(varn,:,:))
    16091604               if (over .EQ. 1) then
     1605                  res@xyMonoDashPattern = True
    16101606                  res@xyDashPattern  = 1
    16111607                  plot_v = gsn_csm_xy(wks,data(varn,:,:),z,res)
     
    16311627               maxiw=max(data(varn,:,:))
    16321628               if (over .EQ. 1) then
    1633                   res@xyDashPattern  = 2
     1629                  res@xyDashPattern = 2
    16341630                  plot_w = gsn_csm_xy(wks,data(varn,:,:),z,res)
    16351631               else
     
    28022798            if (prof3d .EQ. 0) then
    28032799               if (log_z .EQ. 1) then
    2804                   z = f->$vNam(varn+1)$(1:dimz-1)
    2805                else
    2806                   z = f->$vNam(varn+1)$               
     2800                  z = f_att->$vNam(varn+1)$(1:dimz-1)
     2801               else
     2802                  z = f_att->$vNam(varn+1)$               
    28072803               end if
    28082804            else
     
    28352831            res@trYMinF            = min_z
    28362832            res@trYMaxF            = max_z 
     2833            res@xyDashPattern = 0
     2834
    28372835            if (xs .EQ. -1) then
    28382836               res@trXMinF = min(data(varn,:,:))
     
    28622860                     res@trXMaxF = xe 
    28632861                  end if 
     2862                  if (vNam(varn) .EQ. "v")then
     2863                     res@xyDashPattern = 1
     2864                  end if
     2865                  if (vNam(varn) .EQ. "w")then
     2866                     res@xyDashPattern = 2
     2867                  end if
    28642868                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    28652869
     
    28722876                  lgres@lgLabelFont        = "helvetica"   
    28732877                  lgres@lgLabelFontHeightF = font_size_legend           
    2874                   lgres@vpWidthF           = 0.12           
    2875                   lgres@vpHeightF          = 0.1         
     2878                  lgres@vpWidthF           = 0.06           
     2879                  lgres@vpHeightF          = 0.12         
    28762880                  lgres@lgDashIndexes      = (/0,1,2/)
    28772881                  lbid = gsn_create_legend(wks,3,(/"u","v","w"/),lgres)       
     
    29062910                     res@trXMaxF = xe 
    29072911                  end if
    2908 
     2912                  if (vNam(varn) .EQ. "vpt")then
     2913                     res@xyDashPattern = 1
     2914                  end if
     2915                  if (vNam(varn) .EQ. "lpt")then
     2916                     res@xyDashPattern = 2
     2917                  end if
     2918                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     2919
     2920                  ; ***************************************************
     2921                  ; legend for combined plot
     2922                  ; ***************************************************
     2923     
     2924                  lgres                    = True
     2925                  lgMonoDashIndex          = False
     2926                  lgres@lgLabelFont        = "helvetica"   
     2927                  lgres@lgLabelFontHeightF = font_size_legend         
     2928                  lgres@vpWidthF           = 0.07           
     2929                  lgres@vpHeightF          = 0.12         
     2930                  lgres@lgDashIndexes      = (/0,1,2/)
     2931                  lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres)       
     2932
     2933                  amres = True
     2934                  amres@amParallelPosF   = 0.65             
     2935                  amres@amOrthogonalPosF = -0.2           
     2936                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     2937                  overlay(plot(n),plot_pt)
     2938                  pt=1
     2939               else
     2940                  if (prof3d .EQ. 0)then
     2941                     varn=varn+1
     2942                  end if
     2943                  continue       
     2944               end if
     2945            end if           
     2946            if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then
     2947               if (q .EQ. 0) then
     2948                  res@gsnLeftString      = "q, qv and ql"
     2949                  res@tiXAxisString      = "["+unit(varn)+"]"
     2950                  res@gsnRightString     = " "
     2951                  if (xs .EQ. -1) then
     2952                     res@trXMinF = min((/miniq,miniqv,miniql/))
     2953                  else
     2954                     res@trXMinF = xs
     2955                  end if
     2956                  if (xe .EQ. -1) then
     2957                     res@trXMaxF = max((/maxiq,maxiqv,maxiql/))
     2958                  else
     2959                     res@trXMaxF = xe 
     2960                  end if
     2961                  if (vNam(varn) .EQ. "qv")then
     2962                     res@xyDashPattern = 1
     2963                  end if
     2964                  if (vNam(varn) .EQ. "ql")then
     2965                     res@xyDashPattern = 2
     2966                  end if
     2967                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     2968
     2969                  ; ***************************************************
     2970                  ; legend for combined plot
     2971                  ; ***************************************************
     2972     
     2973                  lgres                    = True
     2974                  lgMonoDashIndex          = False
     2975                  lgres@lgLabelFont        = "helvetica"   
     2976                  lgres@lgLabelFontHeightF = font_size_legend           
     2977                  lgres@vpWidthF           = 0.06         
     2978                  lgres@vpHeightF          = 0.12         
     2979                  lgres@lgDashIndexes      = (/0,1,2/)
     2980                  lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres)       
     2981
     2982                  amres = True
     2983                  amres@amParallelPosF   = 0.65             
     2984                  amres@amOrthogonalPosF = -0.2           
     2985                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     2986                  overlay(plot(n),plot_q)
     2987                  q=1
     2988               else
     2989                  if (prof3d .EQ. 0)then
     2990                     varn=varn+1
     2991                  end if
     2992                  continue   
     2993               end if
     2994            end if   
     2995           
     2996            if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then
     2997               if (e .EQ. 0) then
     2998                  res@gsnLeftString      = "e and es"
     2999                  res@tiXAxisString      = "["+unit(varn)+"]"
     3000                  res@gsnRightString     = " "
     3001                  if (xs .EQ. -1) then
     3002                     res@trXMinF = min((/minie,minies/))
     3003                  else
     3004                     res@trXMinF = xs
     3005                  end if
     3006                  if (xe .EQ. -1) then
     3007                     res@trXMaxF = max((/maxie,maxies/))
     3008                  else
     3009                     res@trXMaxF = xe 
     3010                  end if
     3011                  if (vNam(varn) .EQ. "es")then
     3012                     res@xyDashPattern = 1
     3013                  end if
     3014                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3015
     3016                  ; ***************************************************
     3017                  ; legend for combined plot
     3018                  ; ***************************************************
     3019     
     3020                  lgres                    = True
     3021                  lgMonoDashIndex          = False
     3022                  lgres@lgLabelFont        = "helvetica"   
     3023                  lgres@lgLabelFontHeightF = font_size_legend           
     3024                  lgres@vpWidthF           = 0.06           
     3025                  lgres@vpHeightF          = 0.08         
     3026                  lgres@lgDashIndexes      = (/0,1,2/)
     3027                  lbid = gsn_create_legend(wks,2,(/"e","es"/),lgres)       
     3028
     3029                  amres = True
     3030                  amres@amParallelPosF   = 0.65             
     3031                  amres@amOrthogonalPosF = -0.2           
     3032                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3033                  overlay(plot(n),plot_e)
     3034                  e=1
     3035               else
     3036                  if (prof3d .EQ. 0)then
     3037                     varn=varn+1
     3038                  end if
     3039                  continue   
     3040               end if
     3041            end if           
     3042            if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then
     3043               if (km .EQ. 0) then
     3044                  res@gsnLeftString      = "km and kh"
     3045                  res@tiXAxisString      = "["+unit(varn)+"]"
     3046                  res@gsnRightString     = " "
     3047                  if (xs .EQ. -1) then
     3048                     res@trXMinF = min((/minikm,minikh/))
     3049                  else
     3050                     res@trXMinF = xs
     3051                  end if
     3052                  if (xe .EQ. -1) then
     3053                     res@trXMaxF = max((/maxikm,maxikh/))
     3054                  else
     3055                     res@trXMaxF = xe 
     3056                  end if
     3057                  if (vNam(varn) .EQ. "kh")then
     3058                     res@xyDashPattern = 1
     3059                  end if
     3060                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3061
     3062                  ; ***************************************************
     3063                  ; legend for combined plot
     3064                  ; ***************************************************
     3065     
     3066                  lgres                    = True
     3067                  lgMonoDashIndex          = False
     3068                  lgres@lgLabelFont        = "helvetica"   
     3069                  lgres@lgLabelFontHeightF = font_size_legend           
     3070                  lgres@vpWidthF           = 0.06           
     3071                  lgres@vpHeightF          = 0.08         
     3072                  lgres@lgDashIndexes      = (/0,1,2/)
     3073                  lbid = gsn_create_legend(wks,2,(/"km","kh"/),lgres)       
     3074
     3075                  amres = True
     3076                  amres@amParallelPosF   = 0.65             
     3077                  amres@amOrthogonalPosF = -0.2           
     3078                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3079                  overlay(plot(n),plot_km)
     3080                  km=1
     3081               else
     3082                  if (prof3d .EQ. 0)then
     3083                     varn=varn+1
     3084                  end if
     3085                  continue   
     3086               end if
     3087            end if           
     3088           
     3089            if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then
     3090               if (wpup .EQ. 0) then
     3091                  res@gsnLeftString      = "wpup, wsus and wu"
     3092                  res@tiXAxisString      = "["+unit(varn)+"]"
     3093                  res@gsnRightString     = " "
     3094                  if (xs .EQ. -1) then
     3095                     res@trXMinF = min((/miniwpup,miniwsus,miniwu/))
     3096                  else
     3097                     res@trXMinF = xs
     3098                  end if
     3099                  if (xe .EQ. -1) then
     3100                     res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/))
     3101                  else
     3102                     res@trXMaxF = xe 
     3103                  end if 
     3104                  if (vNam(varn) .EQ. "wsus")then
     3105                     res@xyDashPattern = 1
     3106                  end if
     3107                  if (vNam(varn) .EQ. "wu")then
     3108                     res@xyDashPattern = 2
     3109                  end if
     3110                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3111
     3112                  ; ***************************************************
     3113                  ; legend for combined plot
     3114                  ; ***************************************************
     3115     
     3116                  lgres                    = True
     3117                  lgMonoDashIndex          = False
     3118                  lgres@lgLabelFont        = "helvetica"   
     3119                  lgres@lgLabelFontHeightF = font_size_legend           
     3120                  lgres@vpWidthF           = 0.08           
     3121                  lgres@vpHeightF          = 0.12         
     3122                  lgres@lgDashIndexes      = (/0,1,2/)
     3123                  lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres)       
     3124
     3125                  amres = True
     3126                  amres@amParallelPosF   = 0.65             
     3127                  amres@amOrthogonalPosF = -0.2           
     3128                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3129                  overlay(plot(n),plot_wpup)
     3130                  wpup=1
     3131               else
     3132                  if (prof3d .EQ. 0)then
     3133                     varn=varn+1
     3134                  end if
     3135                  continue   
     3136               end if
     3137            end if
     3138            if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then
     3139               if (wpvp .EQ. 0) then
     3140                  res@gsnLeftString      = "wpvp, wsvs and wv"
     3141                  res@tiXAxisString      = "["+unit(varn)+"]"
     3142                  res@gsnRightString     = " "
     3143                  if (xs .EQ. -1) then
     3144                     res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/))
     3145                  else
     3146                     res@trXMinF = xs
     3147                  end if
     3148                  if (xe .EQ. -1) then
     3149                     res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/))
     3150                  else
     3151                     res@trXMaxF = xe 
     3152                  end if
     3153                  if (vNam(varn) .EQ. "wsvs")then
     3154                     res@xyDashPattern = 1
     3155                  end if
     3156                  if (vNam(varn) .EQ. "wv")then
     3157                     res@xyDashPattern = 2
     3158                  end if
     3159                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3160
     3161                  ; ***************************************************
     3162                  ; legend for combined plot
     3163                  ; ***************************************************
     3164     
     3165                  lgres                    = True
     3166                  lgMonoDashIndex          = False
     3167                  lgres@lgLabelFont        = "helvetica"   
     3168                  lgres@lgLabelFontHeightF = font_size_legend           
     3169                  lgres@vpWidthF           = 0.08         
     3170                  lgres@vpHeightF          = 0.12         
     3171                  lgres@lgDashIndexes      = (/0,1,2/)
     3172                  lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres)       
     3173
     3174                  amres = True
     3175                  amres@amParallelPosF   = 0.65             
     3176                  amres@amOrthogonalPosF = -0.2           
     3177                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3178                  overlay(plot(n),plot_wpvp)
     3179                  wpvp=1
     3180               else
     3181                  if (prof3d .EQ. 0)then
     3182                     varn=varn+1
     3183                  end if
     3184                  continue   
     3185               end if
     3186            end if
     3187            if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then
     3188               if (wpptp .EQ. 0) then
     3189                  res@gsnLeftString      = "wpptp, wspts and wpt"
     3190                  res@tiXAxisString      = "["+unit(varn)+"]"
     3191                  res@gsnRightString     = " "
     3192                  if (xs .EQ. -1) then
     3193                     res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/))
     3194                  else
     3195                     res@trXMinF = xs
     3196                  end if
     3197                  if (xe .EQ. -1) then
     3198                     res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/))
     3199                  else
     3200                     res@trXMaxF = xe 
     3201                  end if
     3202                  if (vNam(varn) .EQ. "wspts")then
     3203                     res@xyDashPattern = 1
     3204                  end if
     3205                  if (vNam(varn) .EQ. "wpt")then
     3206                     res@xyDashPattern = 2
     3207                  end if
     3208                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3209
     3210                  ; ***************************************************
     3211                  ; legend for combined plot
     3212                  ; ***************************************************
     3213     
     3214                  lgres                    = True
     3215                  lgMonoDashIndex          = False
     3216                  lgres@lgLabelFont        = "helvetica"   
     3217                  lgres@lgLabelFontHeightF = font_size_legend           
     3218                  lgres@vpWidthF           = 0.09           
     3219                  lgres@vpHeightF          = 0.12         
     3220                  lgres@lgDashIndexes      = (/0,1,2/)
     3221                  lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres)       
     3222
     3223                  amres = True
     3224                  amres@amParallelPosF   = 0.65             
     3225                  amres@amOrthogonalPosF = -0.2           
     3226                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3227                  overlay(plot(n),plot_wpptp)
     3228                  wpptp=1
     3229               else
     3230                  if (prof3d .EQ. 0)then
     3231                     varn=varn+1
     3232                  end if
     3233                  continue   
     3234               end if
     3235            end if
     3236            if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then
     3237               if (wsptsBC .EQ. 0) then
     3238                  res@gsnLeftString      = "wsptsBC and wptBC"
     3239                  res@tiXAxisString      = "["+unit(varn)+"]"
     3240                  res@gsnRightString     = " "
     3241                  if (xs .EQ. -1) then
     3242                     res@trXMinF = min((/miniwsptsBC,miniwptBC/))
     3243                  else
     3244                     res@trXMinF = xs
     3245                  end if
     3246                  if (xe .EQ. -1) then
     3247                     res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/))
     3248                  else
     3249                     res@trXMaxF = xe 
     3250                  end if
     3251                  if (vNam(varn) .EQ. "wptBC")then
     3252                     res@xyDashPattern = 1
     3253                  end if
     3254                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3255
     3256                  ; ***************************************************
     3257                  ; legend for combined plot
     3258                  ; ***************************************************
     3259     
     3260                  lgres                    = True
     3261                  lgMonoDashIndex          = False
     3262                  lgres@lgLabelFont        = "helvetica"   
     3263                  lgres@lgLabelFontHeightF = font_size_legend           
     3264                  lgres@vpWidthF           = 0.1           
     3265                  lgres@vpHeightF          = 0.12         
     3266                  lgres@lgDashIndexes      = (/0,1,2/)
     3267                  lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres)       
     3268
     3269                  amres = True
     3270                  amres@amParallelPosF   = 0.65             
     3271                  amres@amOrthogonalPosF = -0.2           
     3272                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3273                  overlay(plot(n),plot_wsptsBC)
     3274                  wsptsBC=1
     3275               else
     3276                  if (prof3d .EQ. 0)then
     3277                     varn=varn+1
     3278                  end if
     3279                  continue   
     3280               end if 
     3281            end if             
     3282            if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then
     3283               if (wpvptp .EQ. 0) then
     3284                  res@gsnLeftString      = "wpvptp, wsvpts and wvpt"
     3285                  res@tiXAxisString      = "["+unit(varn)+"]"
     3286                  res@gsnRightString     = " "
     3287                  if (xs .EQ. -1) then
     3288                     res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/))
     3289                  else
     3290                     res@trXMinF = xs
     3291                  end if
     3292                  if (xe .EQ. -1) then
     3293                     res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/))
     3294                  else
     3295                     res@trXMaxF = xe 
     3296                  end if
     3297                  if (vNam(varn) .EQ. "wsvpts")then
     3298                     res@xyDashPattern = 1
     3299                  end if
     3300                  if (vNam(varn) .EQ. "wvpt")then
     3301                     res@xyDashPattern = 2
     3302                  end if
     3303                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3304
     3305                  ; ***************************************************
     3306                  ; legend for combined plot
     3307                  ; ***************************************************
     3308     
     3309                  lgres                    = True
     3310                  lgMonoDashIndex          = False
     3311                  lgres@lgLabelFont        = "helvetica"   
     3312                  lgres@lgLabelFontHeightF = font_size_legend           
     3313                  lgres@vpWidthF           = 0.1           
     3314                  lgres@vpHeightF          = 0.12         
     3315                  lgres@lgDashIndexes      = (/0,1,2/)
     3316                  lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres)       
     3317
     3318                  amres = True
     3319                  amres@amParallelPosF   = 0.65             
     3320                  amres@amOrthogonalPosF = -0.2           
     3321                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3322                  overlay(plot(n),plot_wpvptp)
     3323                  wpvptp=1
     3324               else
     3325                  if (prof3d .EQ. 0)then
     3326                     varn=varn+1
     3327                  end if
     3328                  continue   
     3329               end if
     3330            end if
     3331            if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then
     3332               if (wpqp .EQ. 0) then
     3333                  res@gsnLeftString      = "wpqp, wsqs and wq"
     3334                  res@tiXAxisString      = "["+unit(varn)+"]"
     3335                  res@gsnRightString     = " "
     3336                  if (xs .EQ. -1) then
     3337                     res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/))
     3338                  else
     3339                     res@trXMinF = xs
     3340                  end if
     3341                  if (xe .EQ. -1) then
     3342                     res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/))
     3343                  else
     3344                     res@trXMaxF = xe 
     3345                  end if 
     3346                  if (vNam(varn) .EQ. "wsqs")then
     3347                     res@xyDashPattern = 1
     3348                  end if
     3349                  if (vNam(varn) .EQ. "wq")then
     3350                     res@xyDashPattern = 2
     3351                  end if
     3352                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3353
     3354                  ; ***************************************************
     3355                  ; legend for combined plot
     3356                  ; ***************************************************
     3357     
     3358                  lgres                    = True
     3359                  lgMonoDashIndex          = False
     3360                  lgres@lgLabelFont        = "helvetica"   
     3361                  lgres@lgLabelFontHeightF = font_size_legend           
     3362                  lgres@vpWidthF           = 0.08           
     3363                  lgres@vpHeightF          = 0.12         
     3364                  lgres@lgDashIndexes      = (/0,1,2/)
     3365                  lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres)       
     3366
     3367                  amres = True
     3368                  amres@amParallelPosF   = 0.65             
     3369                  amres@amOrthogonalPosF = -0.2           
     3370                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3371                  overlay(plot(n),plot_wpqp)
     3372                  wpqp=1
     3373               else
     3374                  if (prof3d .EQ. 0)then
     3375                     varn=varn+1
     3376                  end if
     3377                  continue   
     3378               end if
     3379            end if
     3380            if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then
     3381               if (wpqvp .EQ. 0) then
     3382                  res@gsnLeftString      = "wpqvp, wsqvs and wqv"
     3383                  res@tiXAxisString      = "["+unit(varn)+"]"
     3384                  res@gsnRightString     = " "
     3385                  if (xs .EQ. -1) then
     3386                     res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/))
     3387                  else
     3388                     res@trXMinF = xs
     3389                  end if
     3390                  if (xe .EQ. -1) then
     3391                     res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/))
     3392                  else
     3393                     res@trXMaxF = xe 
     3394                  end if
     3395                  if (vNam(varn) .EQ. "wsqvs")then
     3396                     res@xyDashPattern = 1
     3397                  end if
     3398                  if (vNam(varn) .EQ. "wqv")then
     3399                     res@xyDashPattern = 2
     3400                  end if
     3401                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3402
     3403                  ; ***************************************************
     3404                  ; legend for combined plot
     3405                  ; ***************************************************
     3406     
     3407                  lgres                    = True
     3408                  lgMonoDashIndex          = False
     3409                  lgres@lgLabelFont        = "helvetica"   
     3410                  lgres@lgLabelFontHeightF = font_size_legend           
     3411                  lgres@vpWidthF           = 0.09           
     3412                  lgres@vpHeightF          = 0.12         
     3413                  lgres@lgDashIndexes      = (/0,1,2/)
     3414                  lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres)       
     3415
     3416                  amres = True
     3417                  amres@amParallelPosF   = 0.65             
     3418                  amres@amOrthogonalPosF = -0.2           
     3419                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3420                  overlay(plot(n),plot_wpqvp)
     3421                  wpqvp=1
     3422               else
     3423                  if (prof3d .EQ. 0)then
     3424                     varn=varn+1
     3425                  end if
     3426                  continue   
     3427               end if
     3428            end if
     3429            if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then
     3430               if (wpsp .EQ. 0) then
     3431                  res@gsnLeftString      = "wpsp, wsss and ws"
     3432                  res@tiXAxisString      = "["+unit(varn)+"]"
     3433                  res@gsnRightString     = " "
     3434                  if (xs .EQ. -1) then
     3435                     res@trXMinF = min((/miniwpsp,miniwsss,miniws/))
     3436                  else
     3437                     res@trXMinF = xs
     3438                  end if
     3439                  if (xe .EQ. -1) then
     3440                     res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/))
     3441                  else
     3442                     res@trXMaxF = xe 
     3443                  end if
     3444                  if (vNam(varn) .EQ. "wsss")then
     3445                     res@xyDashPattern = 1
     3446                  end if
     3447                  if (vNam(varn) .EQ. "ws")then
     3448                     res@xyDashPattern = 2
     3449                  end if
     3450                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3451
     3452                  ; ***************************************************
     3453                  ; legend for combined plot
     3454                  ; ***************************************************
     3455     
     3456                  lgres                    = True
     3457                  lgMonoDashIndex          = False
     3458                  lgres@lgLabelFont        = "helvetica"   
     3459                  lgres@lgLabelFontHeightF = font_size_legend           
     3460                  lgres@vpWidthF           = 0.08           
     3461                  lgres@vpHeightF          = 0.12         
     3462                  lgres@lgDashIndexes      = (/0,1,2/)
     3463                  lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres)       
     3464
     3465                  amres = True
     3466                  amres@amParallelPosF   = 0.65             
     3467                  amres@amOrthogonalPosF = -0.2           
     3468                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3469                  overlay(plot(n),plot_wpsp)
     3470                  wpsp=1
     3471               else
     3472                  if (prof3d .EQ. 0)then
     3473                     varn=varn+1
     3474                  end if
     3475                  continue   
     3476               end if
     3477            end if
     3478            if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then
     3479               if (wpsap .EQ. 0) then
     3480                  res@gsnLeftString      = "wpsap, wssas and wsa"
     3481                  res@tiXAxisString      = "["+unit(varn)+"]"
     3482                  res@gsnRightString     = " "
     3483                  if (xs .EQ. -1) then
     3484                     res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/))
     3485                  else
     3486                     res@trXMinF = xs
     3487                  end if
     3488                  if (xe .EQ. -1) then
     3489                     res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/))
     3490                  else
     3491                     res@trXMaxF = xe 
     3492                  end if
     3493                  if (vNam(varn) .EQ. "wssas")then
     3494                     res@xyDashPattern = 1
     3495                  end if
     3496                  if (vNam(varn) .EQ. "wsa")then
     3497                     res@xyDashPattern = 2
     3498                  end if
     3499                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3500
     3501                  ; ***************************************************
     3502                  ; legend for combined plot
     3503                  ; ***************************************************
     3504     
     3505                  lgres                    = True
     3506                  lgMonoDashIndex          = False
     3507                  lgres@lgLabelFont        = "helvetica"   
     3508                  lgres@lgLabelFontHeightF = font_size_legend           
     3509                  lgres@vpWidthF           = 0.09           
     3510                  lgres@vpHeightF          = 0.12         
     3511                  lgres@lgDashIndexes      = (/0,1,2/)
     3512                  lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres)       
     3513
     3514                  amres = True
     3515                  amres@amParallelPosF   = 0.65             
     3516                  amres@amOrthogonalPosF = -0.2           
     3517                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3518                  overlay(plot(n),plot_wpsap)
     3519                  wpsap=1
     3520               else
     3521                  if (prof3d .EQ. 0)then
     3522                     varn=varn+1
     3523                  end if
     3524                  continue   
     3525               end if
     3526            end if
     3527         
     3528            if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then
     3529               if (us2 .EQ. 0) then
     3530                  res@gsnLeftString      = "us2, vs2 and ws2"
     3531                  res@tiXAxisString      = "["+unit(varn)+"]"
     3532                  res@gsnRightString     = " "
     3533                  if (xs .EQ. -1) then
     3534                     res@trXMinF = min((/minius2,minivs2,miniws2/))
     3535                  else
     3536                     res@trXMinF = xs
     3537                  end if
     3538                  if (xe .EQ. -1) then
     3539                     res@trXMaxF = max((/maxius2,maxivs2,maxiws2/))
     3540                  else
     3541                     res@trXMaxF = xe 
     3542                  end if
     3543                  if (vNam(varn) .EQ. "vs2")then
     3544                     res@xyDashPattern = 1
     3545                  end if
     3546                  if (vNam(varn) .EQ. "ws2")then
     3547                     res@xyDashPattern = 2
     3548                  end if
     3549                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
     3550
     3551                  ; ***************************************************
     3552                  ; legend for combined plot
     3553                  ; ***************************************************
     3554     
     3555                  lgres                    = True
     3556                  lgMonoDashIndex          = False
     3557                  lgres@lgLabelFont        = "helvetica"   
     3558                  lgres@lgLabelFontHeightF = font_size_legend           
     3559                  lgres@vpWidthF           = 0.07           
     3560                  lgres@vpHeightF          = 0.12         
     3561                  lgres@lgDashIndexes      = (/0,1,2/)
     3562                  lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres)       
     3563
     3564                  amres = True
     3565                  amres@amParallelPosF   = 0.65             
     3566                  amres@amOrthogonalPosF = -0.2           
     3567                  annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
     3568                  overlay(plot(n),plot_us2)
     3569                  us2=1
     3570               else
     3571                  if (prof3d .EQ. 0)then
     3572                     varn=varn+1
     3573                  end if
     3574                  continue   
     3575               end if
     3576            end if
     3577           
     3578            if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then
     3579               if (wsususodz .EQ. 0) then
     3580                  res@gsnLeftString      = "wsususodz, wspsodz and wpeodz"
     3581                  res@tiXAxisString      = "["+unit(varn)+"]"
     3582                  res@gsnRightString     = " "
     3583                  if (xs .EQ. -1) then
     3584                     res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/))
     3585                  else
     3586                     res@trXMinF = xs
     3587                  end if
     3588                  if (xe .EQ. -1) then
     3589                     res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/))
     3590                  else
     3591                     res@trXMaxF = xe 
     3592                  end if
     3593                  if (vNam(varn) .EQ. "wspsodz")then
     3594                     res@xyDashPattern = 1
     3595                  end if
     3596                  if (vNam(varn) .EQ. "wpeodz")then
     3597                     res@xyDashPattern = 2
     3598                  end if
    29093599                  plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    29103600
     
    29183608                  lgres@lgLabelFontHeightF = font_size_legend           
    29193609                  lgres@vpWidthF           = 0.12           
    2920                   lgres@vpHeightF          = 0.1         
    2921                   lgres@lgDashIndexes      = (/0,1,2/)
    2922                   lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres)       
    2923 
    2924                   amres = True
    2925                   amres@amParallelPosF   = 0.65             
    2926                   amres@amOrthogonalPosF = -0.2           
    2927                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    2928                   overlay(plot(n),plot_pt)
    2929                   pt=1
    2930                else
    2931                   if (prof3d .EQ. 0)then
    2932                      varn=varn+1
    2933                   end if
    2934                   continue       
    2935                end if
    2936             end if           
    2937             if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then
    2938                if (q .EQ. 0) then
    2939                   res@gsnLeftString      = "q, qv and ql"
    2940                   res@tiXAxisString      = "["+unit(varn)+"]"
    2941                   res@gsnRightString     = " "
    2942                   if (xs .EQ. -1) then
    2943                      res@trXMinF = min((/miniq,miniqv,miniql/))
    2944                   else
    2945                      res@trXMinF = xs
    2946                   end if
    2947                   if (xe .EQ. -1) then
    2948                      res@trXMaxF = max((/maxiq,maxiqv,maxiql/))
    2949                   else
    2950                      res@trXMaxF = xe 
    2951                   end if
    2952 
    2953                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    2954 
    2955                   ; ***************************************************
    2956                   ; legend for combined plot
    2957                   ; ***************************************************
    2958      
    2959                   lgres                    = True
    2960                   lgMonoDashIndex          = False
    2961                   lgres@lgLabelFont        = "helvetica"   
    2962                   lgres@lgLabelFontHeightF = font_size_legend           
    2963                   lgres@vpWidthF           = 0.12           
    2964                   lgres@vpHeightF          = 0.1         
    2965                   lgres@lgDashIndexes      = (/0,1,2/)
    2966                   lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres)       
    2967 
    2968                   amres = True
    2969                   amres@amParallelPosF   = 0.65             
    2970                   amres@amOrthogonalPosF = -0.2           
    2971                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    2972                   overlay(plot(n),plot_q)
    2973                   q=1
    2974                else
    2975                   if (prof3d .EQ. 0)then
    2976                      varn=varn+1
    2977                   end if
    2978                   continue   
    2979                end if
    2980             end if   
    2981            
    2982             if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then
    2983                if (e .EQ. 0) then
    2984                   res@gsnLeftString      = "e and es"
    2985                   res@tiXAxisString      = "["+unit(varn)+"]"
    2986                   res@gsnRightString     = " "
    2987                   if (xs .EQ. -1) then
    2988                      res@trXMinF = min((/minie,minies/))
    2989                   else
    2990                      res@trXMinF = xs
    2991                   end if
    2992                   if (xe .EQ. -1) then
    2993                      res@trXMaxF = max((/maxie,maxies/))
    2994                   else
    2995                      res@trXMaxF = xe 
    2996                   end if
    2997                  
    2998                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    2999 
    3000                   ; ***************************************************
    3001                   ; legend for combined plot
    3002                   ; ***************************************************
    3003      
    3004                   lgres                    = True
    3005                   lgMonoDashIndex          = False
    3006                   lgres@lgLabelFont        = "helvetica"   
    3007                   lgres@lgLabelFontHeightF = font_size_legend           
    3008                   lgres@vpWidthF           = 0.12           
    3009                   lgres@vpHeightF          = 0.1         
    3010                   lgres@lgDashIndexes      = (/0,1,2/)
    3011                   lbid = gsn_create_legend(wks,3,(/"e","es"/),lgres)       
    3012 
    3013                   amres = True
    3014                   amres@amParallelPosF   = 0.65             
    3015                   amres@amOrthogonalPosF = -0.2           
    3016                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3017                   overlay(plot(n),plot_e)
    3018                   e=1
    3019                else
    3020                   if (prof3d .EQ. 0)then
    3021                      varn=varn+1
    3022                   end if
    3023                   continue   
    3024                end if
    3025             end if           
    3026             if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then
    3027                if (km .EQ. 0) then
    3028                   res@gsnLeftString      = "km and kh"
    3029                   res@tiXAxisString      = "["+unit(varn)+"]"
    3030                   res@gsnRightString     = " "
    3031                   if (xs .EQ. -1) then
    3032                      res@trXMinF = min((/minikm,minikh/))
    3033                   else
    3034                      res@trXMinF = xs
    3035                   end if
    3036                   if (xe .EQ. -1) then
    3037                      res@trXMaxF = max((/maxikm,maxikh/))
    3038                   else
    3039                      res@trXMaxF = xe 
    3040                   end if
    3041 
    3042                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3043 
    3044                   ; ***************************************************
    3045                   ; legend for combined plot
    3046                   ; ***************************************************
    3047      
    3048                   lgres                    = True
    3049                   lgMonoDashIndex          = False
    3050                   lgres@lgLabelFont        = "helvetica"   
    3051                   lgres@lgLabelFontHeightF = font_size_legend           
    3052                   lgres@vpWidthF           = 0.12           
    3053                   lgres@vpHeightF          = 0.1         
    3054                   lgres@lgDashIndexes      = (/0,1,2/)
    3055                   lbid = gsn_create_legend(wks,3,(/"km","kh"/),lgres)       
    3056 
    3057                   amres = True
    3058                   amres@amParallelPosF   = 0.65             
    3059                   amres@amOrthogonalPosF = -0.2           
    3060                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3061                   overlay(plot(n),plot_km)
    3062                   km=1
    3063                else
    3064                   if (prof3d .EQ. 0)then
    3065                      varn=varn+1
    3066                   end if
    3067                   continue   
    3068                end if
    3069             end if           
    3070            
    3071             if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then
    3072                if (wpup .EQ. 0) then
    3073                   res@gsnLeftString      = "wpup, wsus and wu"
    3074                   res@tiXAxisString      = "["+unit(varn)+"]"
    3075                   res@gsnRightString     = " "
    3076                   if (xs .EQ. -1) then
    3077                      res@trXMinF = min((/miniwpup,miniwsus,miniwu/))
    3078                   else
    3079                      res@trXMinF = xs
    3080                   end if
    3081                   if (xe .EQ. -1) then
    3082                      res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/))
    3083                   else
    3084                      res@trXMaxF = xe 
    3085                   end if 
    3086 
    3087                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3088 
    3089                   ; ***************************************************
    3090                   ; legend for combined plot
    3091                   ; ***************************************************
    3092      
    3093                   lgres                    = True
    3094                   lgMonoDashIndex          = False
    3095                   lgres@lgLabelFont        = "helvetica"   
    3096                   lgres@lgLabelFontHeightF = font_size_legend           
    3097                   lgres@vpWidthF           = 0.12           
    3098                   lgres@vpHeightF          = 0.1         
    3099                   lgres@lgDashIndexes      = (/0,1,2/)
    3100                   lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres)       
    3101 
    3102                   amres = True
    3103                   amres@amParallelPosF   = 0.65             
    3104                   amres@amOrthogonalPosF = -0.2           
    3105                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3106                   overlay(plot(n),plot_wpup)
    3107                   wpup=1
    3108                else
    3109                   if (prof3d .EQ. 0)then
    3110                      varn=varn+1
    3111                   end if
    3112                   continue   
    3113                end if
    3114             end if
    3115             if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then
    3116                if (wpvp .EQ. 0) then
    3117                   res@gsnLeftString      = "wpvp, wsus and wv"
    3118                   res@tiXAxisString      = "["+unit(varn)+"]"
    3119                   res@gsnRightString     = " "
    3120                   if (xs .EQ. -1) then
    3121                      res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/))
    3122                   else
    3123                      res@trXMinF = xs
    3124                   end if
    3125                   if (xe .EQ. -1) then
    3126                      res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/))
    3127                   else
    3128                      res@trXMaxF = xe 
    3129                   end if
    3130                  
    3131                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3132 
    3133                   ; ***************************************************
    3134                   ; legend for combined plot
    3135                   ; ***************************************************
    3136      
    3137                   lgres                    = True
    3138                   lgMonoDashIndex          = False
    3139                   lgres@lgLabelFont        = "helvetica"   
    3140                   lgres@lgLabelFontHeightF = font_size_legend           
    3141                   lgres@vpWidthF           = 0.12           
    3142                   lgres@vpHeightF          = 0.1         
    3143                   lgres@lgDashIndexes      = (/0,1,2/)
    3144                   lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres)       
    3145 
    3146                   amres = True
    3147                   amres@amParallelPosF   = 0.65             
    3148                   amres@amOrthogonalPosF = -0.2           
    3149                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3150                   overlay(plot(n),plot_wpvp)
    3151                   wpup=0
    3152                else
    3153                   if (prof3d .EQ. 0)then
    3154                      varn=varn+1
    3155                   end if
    3156                   continue   
    3157                end if
    3158             end if
    3159             if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then
    3160                if (wpptp .EQ. 0) then
    3161                   res@gsnLeftString      = "wpptp, wspts and wv"
    3162                   res@tiXAxisString      = "["+unit(varn)+"]"
    3163                   res@gsnRightString     = " "
    3164                   if (xs .EQ. -1) then
    3165                      res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/))
    3166                   else
    3167                      res@trXMinF = xs
    3168                   end if
    3169                   if (xe .EQ. -1) then
    3170                      res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/))
    3171                   else
    3172                      res@trXMaxF = xe 
    3173                   end if
    3174 
    3175                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3176 
    3177                   ; ***************************************************
    3178                   ; legend for combined plot
    3179                   ; ***************************************************
    3180      
    3181                   lgres                    = True
    3182                   lgMonoDashIndex          = False
    3183                   lgres@lgLabelFont        = "helvetica"   
    3184                   lgres@lgLabelFontHeightF = font_size_legend           
    3185                   lgres@vpWidthF           = 0.12           
    3186                   lgres@vpHeightF          = 0.1         
    3187                   lgres@lgDashIndexes      = (/0,1,2/)
    3188                   lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres)       
    3189 
    3190                   amres = True
    3191                   amres@amParallelPosF   = 0.65             
    3192                   amres@amOrthogonalPosF = -0.2           
    3193                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3194                   overlay(plot(n),plot_wpptp)
    3195                   wpptp=1
    3196                else
    3197                   if (prof3d .EQ. 0)then
    3198                      varn=varn+1
    3199                   end if
    3200                   continue   
    3201                end if
    3202             end if
    3203             if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then
    3204                if (wsptsBC .EQ. 0) then
    3205                   res@gsnLeftString      = "wsptsBC and wptBC"
    3206                   res@tiXAxisString      = "["+unit(varn)+"]"
    3207                   res@gsnRightString     = " "
    3208                   if (xs .EQ. -1) then
    3209                      res@trXMinF = min((/miniwsptsBC,miniwptBC/))
    3210                   else
    3211                      res@trXMinF = xs
    3212                   end if
    3213                   if (xe .EQ. -1) then
    3214                      res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/))
    3215                   else
    3216                      res@trXMaxF = xe 
    3217                   end if
    3218 
    3219                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3220 
    3221                   ; ***************************************************
    3222                   ; legend for combined plot
    3223                   ; ***************************************************
    3224      
    3225                   lgres                    = True
    3226                   lgMonoDashIndex          = False
    3227                   lgres@lgLabelFont        = "helvetica"   
    3228                   lgres@lgLabelFontHeightF = font_size_legend           
    3229                   lgres@vpWidthF           = 0.12           
    3230                   lgres@vpHeightF          = 0.1         
    3231                   lgres@lgDashIndexes      = (/0,1,2/)
    3232                   lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres)       
    3233 
    3234                   amres = True
    3235                   amres@amParallelPosF   = 0.65             
    3236                   amres@amOrthogonalPosF = -0.2           
    3237                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3238                   overlay(plot(n),plot_wsptsBC)
    3239                   wsptsBC=1
    3240                else
    3241                   if (prof3d .EQ. 0)then
    3242                      varn=varn+1
    3243                   end if
    3244                   continue   
    3245                end if 
    3246             end if             
    3247             if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then
    3248                if (wpvptp .EQ. 0) then
    3249                   res@gsnLeftString      = "wpvptp, wsvpts and wv"
    3250                   res@tiXAxisString      = "["+unit(varn)+"]"
    3251                   res@gsnRightString     = " "
    3252                   if (xs .EQ. -1) then
    3253                      res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/))
    3254                   else
    3255                      res@trXMinF = xs
    3256                   end if
    3257                   if (xe .EQ. -1) then
    3258                      res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/))
    3259                   else
    3260                      res@trXMaxF = xe 
    3261                   end if
    3262                  
    3263                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3264 
    3265                   ; ***************************************************
    3266                   ; legend for combined plot
    3267                   ; ***************************************************
    3268      
    3269                   lgres                    = True
    3270                   lgMonoDashIndex          = False
    3271                   lgres@lgLabelFont        = "helvetica"   
    3272                   lgres@lgLabelFontHeightF = font_size_legend           
    3273                   lgres@vpWidthF           = 0.12           
    3274                   lgres@vpHeightF          = 0.1         
    3275                   lgres@lgDashIndexes      = (/0,1,2/)
    3276                   lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres)       
    3277 
    3278                   amres = True
    3279                   amres@amParallelPosF   = 0.65             
    3280                   amres@amOrthogonalPosF = -0.2           
    3281                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3282                   overlay(plot(n),plot_wpvptp)
    3283                   wpvptp=1
    3284                else
    3285                   if (prof3d .EQ. 0)then
    3286                      varn=varn+1
    3287                   end if
    3288                   continue   
    3289                end if
    3290             end if
    3291             if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then
    3292                if (wpqp .EQ. 0) then
    3293                   res@gsnLeftString      = "wpqp, wsqs and wq"
    3294                   res@tiXAxisString      = "["+unit(varn)+"]"
    3295                   res@gsnRightString     = " "
    3296                   if (xs .EQ. -1) then
    3297                      res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/))
    3298                   else
    3299                      res@trXMinF = xs
    3300                   end if
    3301                   if (xe .EQ. -1) then
    3302                      res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/))
    3303                   else
    3304                      res@trXMaxF = xe 
    3305                   end if 
    3306                  
    3307                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3308 
    3309                   ; ***************************************************
    3310                   ; legend for combined plot
    3311                   ; ***************************************************
    3312      
    3313                   lgres                    = True
    3314                   lgMonoDashIndex          = False
    3315                   lgres@lgLabelFont        = "helvetica"   
    3316                   lgres@lgLabelFontHeightF = font_size_legend           
    3317                   lgres@vpWidthF           = 0.12           
    3318                   lgres@vpHeightF          = 0.1         
    3319                   lgres@lgDashIndexes      = (/0,1,2/)
    3320                   lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres)       
    3321 
    3322                   amres = True
    3323                   amres@amParallelPosF   = 0.65             
    3324                   amres@amOrthogonalPosF = -0.2           
    3325                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3326                   overlay(plot(n),plot_wpqp)
    3327                   wpqp=1
    3328                else
    3329                   if (prof3d .EQ. 0)then
    3330                      varn=varn+1
    3331                   end if
    3332                   continue   
    3333                end if
    3334             end if
    3335             if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then
    3336                if (wpqvp .EQ. 0) then
    3337                   res@gsnLeftString      = "wpqvp, wsqvs and wqv"
    3338                   res@tiXAxisString      = "["+unit(varn)+"]"
    3339                   res@gsnRightString     = " "
    3340                   if (xs .EQ. -1) then
    3341                      res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/))
    3342                   else
    3343                      res@trXMinF = xs
    3344                   end if
    3345                   if (xe .EQ. -1) then
    3346                      res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/))
    3347                   else
    3348                      res@trXMaxF = xe 
    3349                   end if
    3350                  
    3351                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3352 
    3353                   ; ***************************************************
    3354                   ; legend for combined plot
    3355                   ; ***************************************************
    3356      
    3357                   lgres                    = True
    3358                   lgMonoDashIndex          = False
    3359                   lgres@lgLabelFont        = "helvetica"   
    3360                   lgres@lgLabelFontHeightF = font_size_legend           
    3361                   lgres@vpWidthF           = 0.12           
    3362                   lgres@vpHeightF          = 0.1         
    3363                   lgres@lgDashIndexes      = (/0,1,2/)
    3364                   lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres)       
    3365 
    3366                   amres = True
    3367                   amres@amParallelPosF   = 0.65             
    3368                   amres@amOrthogonalPosF = -0.2           
    3369                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3370                   overlay(plot(n),plot_wpqvp)
    3371                   wpqvp=1
    3372                else
    3373                   if (prof3d .EQ. 0)then
    3374                      varn=varn+1
    3375                   end if
    3376                   continue   
    3377                end if
    3378             end if
    3379             if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then
    3380                if (wpsp .EQ. 0) then
    3381                   res@gsnLeftString      = "wpsp, wsss and ws"
    3382                   res@tiXAxisString      = "["+unit(varn)+"]"
    3383                   res@gsnRightString     = " "
    3384                   if (xs .EQ. -1) then
    3385                      res@trXMinF = min((/miniwpsp,miniwsss,miniws/))
    3386                   else
    3387                      res@trXMinF = xs
    3388                   end if
    3389                   if (xe .EQ. -1) then
    3390                      res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/))
    3391                   else
    3392                      res@trXMaxF = xe 
    3393                   end if
    3394 
    3395                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3396 
    3397                   ; ***************************************************
    3398                   ; legend for combined plot
    3399                   ; ***************************************************
    3400      
    3401                   lgres                    = True
    3402                   lgMonoDashIndex          = False
    3403                   lgres@lgLabelFont        = "helvetica"   
    3404                   lgres@lgLabelFontHeightF = font_size_legend           
    3405                   lgres@vpWidthF           = 0.12           
    3406                   lgres@vpHeightF          = 0.1         
    3407                   lgres@lgDashIndexes      = (/0,1,2/)
    3408                   lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres)       
    3409 
    3410                   amres = True
    3411                   amres@amParallelPosF   = 0.65             
    3412                   amres@amOrthogonalPosF = -0.2           
    3413                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3414                   overlay(plot(n),plot_wpsp)
    3415                   wpsp=1
    3416                else
    3417                   if (prof3d .EQ. 0)then
    3418                      varn=varn+1
    3419                   end if
    3420                   continue   
    3421                end if
    3422             end if
    3423             if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then
    3424                if (wpsap .EQ. 0) then
    3425                   res@gsnLeftString      = "wpsap, wssas and wsa"
    3426                   res@tiXAxisString      = "["+unit(varn)+"]"
    3427                   res@gsnRightString     = " "
    3428                   if (xs .EQ. -1) then
    3429                      res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/))
    3430                   else
    3431                      res@trXMinF = xs
    3432                   end if
    3433                   if (xe .EQ. -1) then
    3434                      res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/))
    3435                   else
    3436                      res@trXMaxF = xe 
    3437                   end if
    3438 
    3439                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3440 
    3441                   ; ***************************************************
    3442                   ; legend for combined plot
    3443                   ; ***************************************************
    3444      
    3445                   lgres                    = True
    3446                   lgMonoDashIndex          = False
    3447                   lgres@lgLabelFont        = "helvetica"   
    3448                   lgres@lgLabelFontHeightF = font_size_legend           
    3449                   lgres@vpWidthF           = 0.12           
    3450                   lgres@vpHeightF          = 0.1         
    3451                   lgres@lgDashIndexes      = (/0,1,2/)
    3452                   lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres)       
    3453 
    3454                   amres = True
    3455                   amres@amParallelPosF   = 0.65             
    3456                   amres@amOrthogonalPosF = -0.2           
    3457                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3458                   overlay(plot(n),plot_wpsap)
    3459                   wpsap=1
    3460                else
    3461                   if (prof3d .EQ. 0)then
    3462                      varn=varn+1
    3463                   end if
    3464                   continue   
    3465                end if
    3466             end if
    3467          
    3468             if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then
    3469                if (us2 .EQ. 0) then
    3470                   res@gsnLeftString      = "us2, vs2 and ws2"
    3471                   res@tiXAxisString      = "["+unit(varn)+"]"
    3472                   res@gsnRightString     = " "
    3473                   if (xs .EQ. -1) then
    3474                      res@trXMinF = min((/minius2,minivs2,miniws2/))
    3475                   else
    3476                      res@trXMinF = xs
    3477                   end if
    3478                   if (xe .EQ. -1) then
    3479                      res@trXMaxF = max((/maxius2,maxivs2,maxiws2/))
    3480                   else
    3481                      res@trXMaxF = xe 
    3482                   end if
    3483                  
    3484                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3485 
    3486                   ; ***************************************************
    3487                   ; legend for combined plot
    3488                   ; ***************************************************
    3489      
    3490                   lgres                    = True
    3491                   lgMonoDashIndex          = False
    3492                   lgres@lgLabelFont        = "helvetica"   
    3493                   lgres@lgLabelFontHeightF = font_size_legend           
    3494                   lgres@vpWidthF           = 0.12           
    3495                   lgres@vpHeightF          = 0.1         
    3496                   lgres@lgDashIndexes      = (/0,1,2/)
    3497                   lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres)       
    3498 
    3499                   amres = True
    3500                   amres@amParallelPosF   = 0.65             
    3501                   amres@amOrthogonalPosF = -0.2           
    3502                   annoid1 = gsn_add_annotation(plot(n),lbid,amres)                 
    3503                   overlay(plot(n),plot_us2)
    3504                   us2=1
    3505                else
    3506                   if (prof3d .EQ. 0)then
    3507                      varn=varn+1
    3508                   end if
    3509                   continue   
    3510                end if
    3511             end if
    3512            
    3513             if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then
    3514                if (wsususodz .EQ. 0) then
    3515                   res@gsnLeftString      = "wsususodz, wspsodz and ws2"
    3516                   res@tiXAxisString      = "["+unit(varn)+"]"
    3517                   res@gsnRightString     = " "
    3518                   if (xs .EQ. -1) then
    3519                      res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/))
    3520                   else
    3521                      res@trXMinF = xs
    3522                   end if
    3523                   if (xe .EQ. -1) then
    3524                      res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/))
    3525                   else
    3526                      res@trXMaxF = xe 
    3527                   end if
    3528                  
    3529                   plot(n) =  gsn_csm_xy(wks,data(varn,:,:),z,res)
    3530 
    3531                   ; ***************************************************
    3532                   ; legend for combined plot
    3533                   ; ***************************************************
    3534      
    3535                   lgres                    = True
    3536                   lgMonoDashIndex          = False
    3537                   lgres@lgLabelFont        = "helvetica"   
    3538                   lgres@lgLabelFontHeightF = font_size_legend           
    3539                   lgres@vpWidthF           = 0.12           
    3540                   lgres@vpHeightF          = 0.1         
     3610                  lgres@vpHeightF          = 0.12         
    35413611                  lgres@lgDashIndexes      = (/0,1,2/)
    35423612                  lbid = gsn_create_legend(wks,3,(/"wsususodz","wspsodz","wpeodz"/),lgres)       
     
    35903660            if (prof3d .EQ. 0) then
    35913661               if (log_z .EQ. 1) then
    3592                   z = f->$vNam(varn+1)$(1:dimz-1)
    3593                else
    3594                   z = f->$vNam(varn+1)$               
     3662                  z = f_att->$vNam(varn+1)$(1:dimz-1)
     3663               else
     3664                  z = f_att->$vNam(varn+1)$               
    35953665               end if
    35963666            else
     
    37183788   if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
    37193789      gsn_panel(wks,plot_,(/n,1/),resP)
     3790      print(" ")
     3791      print("Outputs to .eps or .epsi have only one frame")
     3792      print(" ")
    37203793   else   
    37213794      do i = 0,n-1, no_rows*no_columns
     
    37233796            gsn_panel(wks,plot_(i:n-1),(/no_rows,no_columns/),resP)
    37243797         else
    3725            gsn_panel(wks,plot_(i:i+no_rows*no_columns-1),(/no_rows,no_columns/),resP)
     3798            gsn_panel(wks,plot_(i:i+no_rows*no_columns-1),(/no_rows,no_columns/),resP)
    37263799         end if
    37273800      end do
  • palm/trunk/SCRIPTS/NCL/spectra.ncl

    r218 r250  
    639639   if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
    640640      gsn_panel(wks,plot,(/n,1/),resP)
     641      print(" ")
     642      print("Outputs to .eps or .epsi have only one frame")
     643      print(" ")
    641644   else   
    642645      do i = 0,n-1, no_rows*no_columns
  • palm/trunk/SCRIPTS/NCL/timeseries.ncl

    r218 r250  
    12431243   if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
    12441244      gsn_panel(wks_ps,plot_ps(1:n),(/n,1/),resP)
     1245      print(" ")
     1246      print("Outputs to .eps or .epsi have only one frame")
     1247      print(" ")
    12451248   else
    12461249      do np = 1,n,no_rows*no_columns   
Note: See TracChangeset for help on using the changeset viewer.