gphys 1.2.2.1 → 1.4.3

data/bin/gpview

@@ -66,6 +66,16 @@ if you want to know more description, please read
 :--map_axis [uxc,uyc,rot]:
   (for 10<=itr<=33) set mapping parameters which are fed in DCL.umpcnt. 
 
+:--map_radius <radius>:
+  (for itr>=20) set clipping radius (degree) around the tangential point. 
+  Deafault=90.
+
+:--xcoord [xcoord]
+  name of x-coordinate (use for associate coordinates)
+
+:--ycoord [ycoord]
+  name of y-coordinate (use for associate coordinates)
+
 :--title:
   set title of figure
 
@@ -96,9 +106,16 @@ if you want to know more description, please read
 :--mean <dim>:
   mean along axis <dim>. 
 
+:--stddev <dim>:
+  standard deviation along axis <dim>.
+
 :--eddy <dim>:
   deviation from mean along axis <dim>. 
 
+:--diff:
+  calculate the difference of two variables.
+  you must provide two gturls whose size and dimension are same.  
+
 :--m, --map <map_type>: 
   plot map. itr number must be set.  this option is neglect if itr
   number is 1-4.  abailable map type is coast_world, border_world,
@@ -118,6 +135,11 @@ if you want to know more description, please read
               in dclext_datetime_ax.rb)
     "ymd" (=> like "h" but for y-m-d type using DCL.uc[xy]acl)
 
+:-xx:yyyy=value:
+  specify Dennou Club Library internal variable. 
+  'xx' is a package name prefix of DCL (eg. sg=>sgpack).
+  'yy' is a variable name of DCL.
+  'value' is a value of the DCL variable.
 
 == LINE/MARK OPTIONS
 
@@ -138,8 +160,11 @@ if you want to know more description, please read
   set min/max value for line plot. min or max must be set.
 
 :--overplot
-  set number of lines on each figure
+  set number of lines on each figure.
 
+:--overplot_color, --Opc
+  multiple lines are distinguished with color rather than line type.
+  
 == CONTOUR/TONE OPTIONS
 
 :--shade:
@@ -183,8 +208,13 @@ if you want to know more description, please read
 :--patterns:           
   set each patterns for tone plot. 
 
-:--tonf:
-  use tonf instead of tone.
+:--tone [a|e|f|b|c]:
+  set tone subroutine:
+    a (=> tone routine is selected automatically depending on the datasize)
+    e (=> DCL.uetone is used)
+    f (=> DCL.uetonf is used)
+    b (=> DCL.uetonb is used)
+    c (=> DCL.uetonc is used)
 
 :--nocolorbar:         
   do not draw color bar
@@ -192,6 +222,9 @@ if you want to know more description, please read
 :--nozero:             
   do not draw zero contour
 
+:--udsfmt [strings]:
+  change contour label format. see UDCNTR/DCL manual for the format.
+
 == EXAMPLES
 
 For a single GPhys variable, 
@@ -211,6 +244,9 @@ For a single GPhys variable,
     % gpview --operation log10 data.nc@temp,lon=0
     % gpview --time_ax false data.nc@temp,lon=0,lat=0
 
+Example to set DCL internal variable,
+    % gpview -ud:lmsg=false --itr 5 --exch data.nc@temp,lon=0
+
 For multiple GPhys variables, 
     % gpview data1.nc@temp data2.nc@temp
     % gpview data1.nc@temp,lon=130:150 data2.nc@temp,lon=150:170
@@ -256,6 +292,19 @@ For multiple GPhys variables,
   2011/11/04  S Takehiro (option --time_ax implemented)
   2011/11/20  S Takehiro (option --eddy implemented)
   2012/02/19  S Takehiro (description for gturl format updated)
+  2013/01/29  T Tanigawa && S Takehiro (calculation of aspect ratio 
+                                        of viewport improved)
+  2013/01/29  S Nishizawa && S Takehiro (sequential operation implemented)
+  2014/02/21  Y Naito && S Takehiro (implement DCL internal variable setting,
+                                     add option --udsfmt)
+  2014/02/21  S Noda && S Takehiro (add option --stddev)
+  2014/02/23  S Takehiro (option --tonf is changed to --tone)
+  2014/03/05  S Takehiro (add option --map_radius)
+  2014/03/10  Y Kawai && S Takehiro (setting labels for contours)
+  2014/04/06  S Takehiro (option --map_radius bug fix)
+  2014/06/27  S Nishizawa && S Takehiro (add options --xcoord, --ycoord)
+  2014/11/28  H Kashimura && S Takehiro (add option --diff, --tone bug fix)
+  2014/11/28  Y O Takahashi && S Takehiro (add option --overplot_color, --Opc)
 
 =end
 #################################################
@@ -315,7 +364,7 @@ def each_along_dims(gphys, loopdim)
 end
 
 
-def draw_setup(gp)
+def draw_setup
 
   # set missing value
   DCLExt.gl_set_params('lmiss'=>true)
@@ -347,34 +396,20 @@ def draw_setup(gp)
                                     'min' => ( min_srange  || min_range ),
                                     'max' => ( max_srange  || max_range )
                                  )
+
   if ( $OPT_clevels || $OPT_levels )
+    $OPT_label=($OPT_clevels || $OPT_levels).split(',')
     $OPT_clevels=($OPT_clevels || $OPT_levels).split(',').map!{|v| v.to_f }
   end
 
   if ( $OPT_slevels || $OPT_levels )
     $OPT_slevels=($OPT_slevels||$OPT_levels).split(',').map!{|v| v.to_f }
-    p $OPT_slevels
   end
 
   if ( $OPT_patterns )
     $OPT_patterns=$OPT_patterns.split(',').map!{|v| v.to_f }
   end
 
-  # judge draw kind
-  gp_rank = gp.rank
-  gp_rank = gp_rank - 1 if ( $OPT_animate || $OPT_anim )
-  if ($OPT_mark)
-    draw_flag = "mark"
-  elsif ($OPT_line || gp_rank == 1)
-    draw_flag = "line"
-  elsif (!$OPT_line && gp_rank >= 2) && !$OPT_noshade && $OPT_nocont 
-    draw_flag = "nocont"
-  elsif (!$OPT_line && gp_rank >= 2) && $OPT_noshade && !$OPT_nocont 
-    draw_flag = "noshade"
-  elsif (!$OPT_line && gp_rank >= 2) && !$OPT_noshade && !$OPT_nocont
-    draw_flag = "full"
-  end  
-
   # similar projection
   if ($OPT_similar)
     if /([\d\-.]*),([\d\-.]*),([\d\-.]*)/ =~ $OPT_similar
@@ -399,6 +434,13 @@ def draw_setup(gp)
     GGraph.set_fig('map_axis'=>map_axis)
   end
 
+  # clipping parameter
+  if ( $OPT_map_radius )
+    map_radius=$OPT_map_radius.to_f
+  else
+    map_radius=90.0
+  end
+  GGraph.set_fig('map_radius'=>map_radius)
 
   # map
   if ( $OPT_m || $OPT_map)
@@ -413,8 +455,6 @@ def draw_setup(gp)
     GGraph.set_axes('time_ax'=>$OPT_time_ax)
   end
 
-  return draw_flag
-
 end
 
 def draw(gp, draw_flag)
@@ -433,6 +473,17 @@ def draw(gp, draw_flag)
                   "min" => __split_range($OPT_range)[0],
                   "max" => __split_range($OPT_range)[1]
                   )
+    elsif ( $OPT_overplot_color || $OPT_Opc )
+      GGraph.line(gp, 
+                  false,
+                  "title"=>$OPT_title,
+                  "index"=>($OPT_index.to_i+$Overplot*10),
+                  "type" =>($OPT_type ||1),
+                  "exchange"=>$OPT_exch,
+                  "annotate"=>$annotate,
+                  "min" => __split_range($OPT_range)[0],
+                  "max" => __split_range($OPT_range)[1]
+                  )
     else
       GGraph.line(gp, 
                   false,
@@ -489,13 +540,21 @@ def draw(gp, draw_flag)
                 "transpose"=>$OPT_exch,
                 "levels"=>$OPT_slevels,
                 "patterns"=>$OPT_patterns,
-                "tonf"=>$OPT_tonf
+                "auto"=>$auto,
+                "tonf"=>$tonf,
+                "tonb"=>$tonb,
+                "tonc"=>$tonc,
+                "xcoord"=>$OPT_xcoord,
+                "ycoord"=>$OPT_ycoord
                 )
     GGraph.contour(gp, 
 		   false,
 		   "transpose"=>$OPT_exch,
 		   "levels"=>$OPT_clevels,
-		   "nozero"=>$OPT_nozero
+		   "nozero"=>$OPT_nozero,
+                   "label"=>$OPT_label,
+                   "xcoord"=>$OPT_xcoord,
+                   "ycoord"=>$OPT_ycoord
 		   )
   when "nocont"
     GGraph.tone(gp,
@@ -505,7 +564,12 @@ def draw(gp, draw_flag)
                 "transpose"=>$OPT_exch,
                 "levels"=>$OPT_slevels,
                 "patterns"=>$OPT_patterns,
-                "tonf"=>$OPT_tonf
+                "auto"=>$auto,
+                "tonf"=>$tonf,
+                "tonb"=>$tonb,
+                "tonc"=>$tonc,
+                "xcoord"=>$OPT_xcoord,
+                "ycoord"=>$OPT_ycoord
                 )
   when "noshade"
     mj = DCL.udpget('indxmj')
@@ -513,11 +577,13 @@ def draw(gp, draw_flag)
     GGraph.contour(gp, 
                    true, 
                    "title"    =>$OPT_title,
-                   "label"    =>true,
                    "annotate"=>$annotate,
 		   "transpose"=>$OPT_exch,
                    "levels"=>$OPT_clevels, 
-                   "nozero"=>$OPT_nozero
+                   "nozero"=>$OPT_nozero,
+                   "label"=>$OPT_label,
+                   "xcoord"=>$OPT_xcoord,
+                   "ycoord"=>$OPT_ycoord
                    )
   end
 
@@ -529,8 +595,6 @@ def draw(gp, draw_flag)
                       )
   end
 
-
-
 end
 
 
@@ -548,7 +612,7 @@ def set_vpsize( default_vp, aspect=2.0 )
   # center grid of viewport
   cen_of_vp = [ x0 + hlength/2.0, y0 + vlength/2.0  ] 
 
-  if aspect <= 2.0
+  if aspect <= hlength / vlength
     hlength = vlength * aspect
     x0 = cen_of_vp[0] - hlength/2.0
     x1 = cen_of_vp[0] + hlength/2.0
@@ -584,10 +648,60 @@ def __split_range(range)
   return min, max
 end
 
+def check_dclopts
+
+  indices = []
+    ARGV.each_index{|i|
+    if /^-(.)(.):(.*)=(.*)/ =~ ARGV[i]
+      indices << i
+    end
+    }
+
+  dclopts = []
+    indices.reverse_each{|i|
+    dclopts << ARGV[i]
+    ARGV.slice!(i)
+    }
+
+  dclopts.each{|opt|
+    pkg = opt[1..2]
+    name = opt[4..-1].split('=')[0]
+    value = opt[4..-1].split('=')[1]
+    dcl_set_params(pkg,name,value)
+  }
+  
+end
+
+def dcl_set_params(pkg,name,value)
+
+  set = 'stx'
+  case name
+  when /^c/i
+    eval( "DCL.#{pkg}c#{set}(name,value.to_s)" )
+  when /^l/i
+    if /(.*)(T|t)(.*)/ =~ value
+      eval( "DCL.#{pkg}l#{set}(name,true)" )
+    else
+    if /(.*)(F|f)(.*)/ =~ value
+      eval( "DCL.#{pkg}l#{set}(name,false)" )
+    else
+      raise "value of logical parameter must include 't' or 'f'"
+    end
+    end
+  when /^[i-n]/i
+    eval( "DCL.#{pkg}i#{set}(name,value.to_i)" )
+  else
+    eval( "DCL.#{pkg}r#{set}(name,value.to_f)" )
+  end
+  
+end
 
 #####################################################
 ###++++++           Main Routine            ++++++###
 
+## options for DCL
+check_dclopts
+
 ## parse options
 parser = GetoptLong.new
 parser.set_options(
@@ -598,6 +712,9 @@ parser.set_options(
                    ['--itr',                      GetoptLong::REQUIRED_ARGUMENT],
                    ['--similar',                  GetoptLong::REQUIRED_ARGUMENT],
                    ['--map_axis',                 GetoptLong::REQUIRED_ARGUMENT],
+                   ['--map_radius',               GetoptLong::REQUIRED_ARGUMENT],
+                   ['--xcoord',                   GetoptLong::REQUIRED_ARGUMENT],
+                   ['--ycoord',                   GetoptLong::REQUIRED_ARGUMENT],
                    ['--title',                    GetoptLong::REQUIRED_ARGUMENT],
                    ['--aspect',                   GetoptLong::REQUIRED_ARGUMENT],
                    ['--anim',                     GetoptLong::REQUIRED_ARGUMENT],
@@ -613,7 +730,10 @@ parser.set_options(
                    ['--reverse',                  GetoptLong::NO_ARGUMENT],
                    ['--Gr',                       GetoptLong::NO_ARGUMENT],
                    ['--mean',                     GetoptLong::REQUIRED_ARGUMENT],
+                   ['--stddev',                   GetoptLong::REQUIRED_ARGUMENT
+],
                    ['--eddy',                     GetoptLong::REQUIRED_ARGUMENT],
+                   ['--diff',                     GetoptLong::NO_ARGUMENT],                   
                    ['--map',                      GetoptLong::REQUIRED_ARGUMENT],
                    ['--m',                        GetoptLong::NO_ARGUMENT],
                    ['--operation',                GetoptLong::REQUIRED_ARGUMENT],
@@ -624,6 +744,8 @@ parser.set_options(
                    ['--index',                    GetoptLong::REQUIRED_ARGUMENT],
                    ['--type',                     GetoptLong::REQUIRED_ARGUMENT],
                    ['--overplot',                 GetoptLong::REQUIRED_ARGUMENT],
+                   ['--overplot_color',           GetoptLong::NO_ARGUMENT],
+                   ['--Opc',                      GetoptLong::NO_ARGUMENT],
                    ### tone or cont option ###
                    ['--nocont',                   GetoptLong::NO_ARGUMENT],
                    ['--noshade',                  GetoptLong::NO_ARGUMENT],
@@ -638,7 +760,8 @@ parser.set_options(
                    ['--clevels',                  GetoptLong::REQUIRED_ARGUMENT],
                    ['--slevels',                  GetoptLong::REQUIRED_ARGUMENT],
                    ['--patterns',                 GetoptLong::REQUIRED_ARGUMENT],
-                   ['--tonf',                     GetoptLong::NO_ARGUMENT],
+                   ['--tone',                     GetoptLong::REQUIRED_ARGUMENT],
+                   ['--udsfmt',                   GetoptLong::REQUIRED_ARGUMENT],
                    ['--nocolorbar',               GetoptLong::NO_ARGUMENT],
                    ['--nozero',                   GetoptLong::NO_ARGUMENT],
                    ['--help',                     GetoptLong::NO_ARGUMENT]
@@ -669,6 +792,7 @@ if ($OPT_noannotate)
 else
   $annotate = true
 end
+
 $Overplot_max = ( $OPT_overplot.to_i || 1 )
 $Overplot = 1
 
@@ -678,12 +802,35 @@ else
   $colorbar = true
 end
 
+if ($OPT_tone) 
+  case $OPT_tone
+  when "a"
+    $auto, $tonf, $tonb, $tonc = true,false,false,false
+  when "e"
+    $auto, $tonf, $tonb, $tonc = false,false,false,false
+  when "f"
+    $auto, $tonf, $tonb, $tonc = false,true,false,false
+  when "b"
+    $auto, $tonf, $tonb, $tonc = false,false,true,false
+  when "c"
+    $auto, $tonf, $tonb, $tonc = false,false,false,true
+  else
+    raise "The value of option --tone should be 'a','e','f','b' or 'c'."
+  end
+else
+  $auto, $tonf, $tonb, $tonc = true,false,false,false
+end
+
 ## decide VIEWPORT
 $VIEWPORT = set_vpsize( VIEWPORT, ($OPT_aspect||2.0) )
 
 ## tune the size of axis parameters.
 DCL.uzfact(0.7)
 
+## set the format of contour labels.
+udsfmt = ($OPT_udsfmt||DCL.udqfmt())
+DCL.udsfmt(udsfmt)
+
 ## draw figure
 loopdim   = ( $OPT_animate || $OPT_anim ) 
 loopdim = loopdim.to_i if loopdim.to_i.to_s == loopdim
@@ -691,7 +838,7 @@ loopdim = loopdim.to_i if loopdim.to_i.to_s == loopdim
 ### set colormap
 DCL.sgscmn($OPT_clrmap||1)
 ## open work station
-DCL.gropn($OPT_wsn||1)
+DCL.gropn($OPT_wsn||4)
 
 ## open netcdf variables
 
@@ -703,38 +850,94 @@ while ARGV[0] do
   print "  Reading #{gturl}\n"
   ARGV.shift 
 
-## mean along any axis
+  ## for case of calculating difference of two gphys object
+  if ($OPT_diff)
+    raise "--diff option must be used with even numbers (2, 4, 6,...) of gturls" if ARGV[0] == nil
+    prev_gturl = gturl
+    gturl = ARGV[0]
+    gturl = gturl+'@'+$OPT_var if $OPT_var
+    print "  Reading #{gturl}\n"
+    gp = gp - GPhys::IO.open_gturl(gturl)
+    print "  Taking difference: #{prev_gturl} - #{gturl}\n"
+    gturl = prev_gturl+" - "+gturl
+    ARGV.shift 
+  end
+  
+  ## mean along any axis
   if ($OPT_mean)
-    dims = ($OPT_mean).split(/\s*,\s*/)
-    dims.each{|dim|
-      dim = dim.to_i if dim.to_i.to_s == dim
-      gp = gp.mean(dim)
-    }
+    dims_mean = ($OPT_mean).split(/\s*,\s*/)
+    dims_mean = dims_mean.map{|dim| dim.to_i.to_s == dim ? dim.to_i : dim}
   end
 
-## deviation from mean along any axis
-  if ($OPT_eddy)
-    dims = ($OPT_eddy).split(/\s*,\s*/)
-    dims.each{|dim|
-      dim = dim.to_i if dim.to_i.to_s == dim
-      gp = gp.eddy(dim)
-    }
+  ## standard deviation along any axis
+  if ($OPT_stddev)
+    dims_stddev = ($OPT_stddev).split(/\s*,\s*/)
+    dims_stddev = dims_stddev.map{|dim| dim.to_i.to_s == dim ? dim.to_i : dim}
   end
 
-## operation of a mathematical function
-  if ($OPT_operation)
-    eval "gp = gp.#{$OPT_operation}"
+  ## deviation from mean along any axis
+  if ($OPT_eddy)
+    dims_eddy = ($OPT_eddy).split(/\s*,\s*/)
+    dims_eddy = dims_eddy.map{|dim| dim.to_i.to_s == dim ? dim.to_i : dim}
   end
 
   GGraph.margin_info($0, gturl) if $annotate  # draw margin infomation
-  kind_of_fig = draw_setup(gp)                # determine figure kind
+  draw_setup                                  # determine figure kind
+
+
+  kind_of_fig = nil
+  proc = Proc.new do |g|
+      ## mean along any axis
+      if ($OPT_mean)
+        dims_mean.each{|dim|
+          g = g.mean(dim)
+        }
+      end
+
+      ## standard deviation along any axis
+      if ($OPT_stddev)
+        dims_stddev.each{|dim|
+          g = g.stddev(dim)
+        }
+      end
+
+      ## deviation from mean along any axis
+      if ($OPT_eddy)
+        dims_eddy.each{|dim|
+          g = g.eddy(dim)
+        }
+      end
+
+      ## operation of a mathematical function
+      if ($OPT_operation)
+        eval "g = g.#{$OPT_operation}"
+      end
+
+      # judge draw kind
+      unless kind_of_fig
+        gp_rank = g.rank
+        if ($OPT_mark)
+          kind_of_fig = "mark"
+        elsif ($OPT_line || gp_rank == 1)
+          kind_of_fig = "line"
+        elsif (!$OPT_line && gp_rank >= 2) && !$OPT_noshade && $OPT_nocont 
+          kind_of_fig = "nocont"
+        elsif (!$OPT_line && gp_rank >= 2) && $OPT_noshade && !$OPT_nocont 
+          kind_of_fig = "noshade"
+        elsif (!$OPT_line && gp_rank >= 2) && !$OPT_noshade && !$OPT_nocont
+          kind_of_fig = "full"
+        end
+      end  
+
+      draw(g, kind_of_fig)
+  end
 
   if loopdim           # animation
     each_along_dims(gp, loopdim){|gp_subset|
-      draw(gp_subset, kind_of_fig)
+      proc.call(gp_subset)
     }
   else
-    draw( gp, kind_of_fig )                # single figure
+    proc.call(gp)
   end
 
 end