gphys 1.5.5 → 1.5.7

data/bin/gpview

@@ -3,17 +3,17 @@
 =begin
 = NAME
 
-gpview - quick viewer for the values of a variable specified by a gtool4-type URL. 
+gpview - quick viewer for the values of a variable specified by a gtool4-type URL.
 
 
 (1) for 1-dimensional variable, make line plot.
 (2) for 2-dimensional variable, make contour/tone plot.
-(3) for 3/more-dimensional variable, make contour/tone plot, 
+(3) for 3/more-dimensional variable, make contour/tone plot,
     select first 2D. possible to make animation.
-    
+
 = USAGE
 
-The first command form is 
+The first command form is
 
     % gpview [options] gturl1 gturl2 gturl3 ...
 
@@ -21,11 +21,20 @@ where the format of the gturl is
 
     path[@|/]varname[,dimname=pos1[:pos2[:thinning_intv]][,dimname=...]]
 
-The second commnad form is 
+When the slice parameters are the same throughout the all Gphys
+variables, the following form can be used:
+
+     % gpview --slice <slice> [options] file1@var1 file2@var2 ...
+
+where the format of the slice is
+
+     dimname=pos1[:pos2[:thinning_intv]][,dimname=...]]
+
+The second commnad form is
 
     % gpview --var url [options] gtfile1 gtfile2 ...
 
-where the format of the url is 
+where the format of the url is
 
     varname[,dimname=pos1[:pos2[:thinning_intv]][,dimname=...]]
 
@@ -37,19 +46,22 @@ if you want to know more description, please read
 
 == GLOBAL OPTIONS
 
-:--help: 
-  print this message. 
+:--help:
+  print this message.
 
 :--var url:
   set the variable name and slicing parameters.
 
+:--slice url:
+  set the slicing parameters.
+
 :--wsn [1-4]:
   set work staion number. each number represent output device:
-  * for dcl-6, 
+  * for dcl-6,
      1 : X window
      2.: PDF file (named dcl.pdf).
-  * for dcl-5-?, 
-     1 : X window, 
+  * for dcl-5-?,
+     1 : X window,
      2 : PS file. (named dcl.ps)
      3 : Tektronix terminal
      4 : GTK Windows
@@ -59,19 +71,19 @@ if you want to know more description, please read
 
 :--itr [1-4,5-7,10-15,20-23,30-33]:
   set axis scale. default is 1.
-  1 : linear scale for x/y axis 
+  1 : linear scale for x/y axis
   2 : linear scale for x , log scale for y axis
   3 : log scale for x , linear scale for y axis
   4 : log scale for x/y axis
 
 :--similar [simfac,vxoff,vyoff]:
-  (for 5<=itr<=7) set similarity parameters which are fed in DCL.grssim. 
+  (for 5<=itr<=7) set similarity parameters which are fed in DCL.grssim.
 
 :--map_axis [uxc,uyc,rot]:
-  (for 10<=itr<=33) set mapping parameters which are fed in DCL.umpcnt. 
+  (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. 
+  (for itr>=20) set clipping radius (degree) around the tangential point.
   Deafault=90.
 
 :--xcoord [xcoord]
@@ -87,16 +99,16 @@ if you want to know more description, please read
   set aspect ratio of Viewport. default is 2.0.
 
 :--noannotate:
-  not draw annotations. 
+  not draw annotations.
 
 :--animate/anim <dim>:
-  plot animation along <dim>. <dim> must be name of dimension. 
+  plot animation along <dim>. <dim> must be name of dimension.
 
-:--alternate, --Ga: 
+:--alternate, --Ga:
   enable to backing store.
 
 :--nowait, --Gw:
-  not wait for any actions if animate 
+  not wait for any actions if animate
 
 :--delay <number>:
   wait <number> seconds for each new page drawing of animation
@@ -111,39 +123,42 @@ if you want to know more description, please read
   exchange(transpose) x/y axis.
 
 :--mean <dim>:
-  mean along axis <dim>. 
+  mean along axis <dim>.
+
+:--sum <dim>:
+  sum up along axis <dim>.
 
 :--stddev <dim>:
   standard deviation along axis <dim>.
 
 :--eddy <dim>:
-  deviation from mean along axis <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.  
+  you must provide two gturls whose size and dimension are same.
 
-:--m, --map <map_type>: 
+:--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,
   plate_world, state_usa, coast_japan, pref_japan
 
-:--operation <math_func>: 
+:--operation <math_func>:
   operation of the specified math function on the data.
-  <math_func> should be a math function with one argument 
+  <math_func> should be a math function with one argument
   such as log10, sqrt, sin, etc.
 
-:--time_ax [nil|false|h|ymd]: 
-  specify type of calendar-type time axis: 
+:--time_ax [nil|false|h|ymd]:
+  specify type of calendar-type time axis:
     nil   (=> auto slection)
-    false (=> do not use the time axis even if 
+    false (=> do not use the time axis even if
               the units of the axis is a time one with since field)
     "h"   (=> like nil, but always use the hour-resolving datetime_ax method
               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. 
+  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.
@@ -171,7 +186,7 @@ if you want to know more description, please read
 
 :--overplot_color, --Opc
   multiple lines are distinguished with color rather than line type.
-  
+
 == CONTOUR/TONE OPTIONS
 
 :--shade:
@@ -192,28 +207,46 @@ if you want to know more description, please read
 :--srange:
   set min/max value for tone plot. this is more dominant than --interval/int
 
-:--interval,--int [num]: 
+:--interval,--int [num]:
   set interval value for contour/tone plot. set the number of lines if
   you set negative value.
 
 :--cint:
-  set interval value for contour plot. this is more dominant 
+  set interval value for contour plot. this is more dominant
   than --interval/int
 
-:--sint:               
+:--sint:
   set interval value for tone plot. this is more dominant than --interval/int.
-                           
-:--levels:             
-  set values of contour/tone levels. 
 
-:--clevels:            
-  set values of contour levels. 
+:--levels:
+  set values of contour/tone levels.
+
+:--clevels:
+  set values of contour levels.
+
+:--slevels:
+  set tone of contour levels.
+
+:--patterns:
+  set each patterns for tone plot.
 
-:--slevels:           
-  set tone of contour levels. 
+:--clog:
+  approximately log-scaled levels for contour plot.
 
-:--patterns:           
-  set each patterns for tone plot. 
+:--slog:
+  approximately log-scaled levels for tone plot.
+
+:--clog_cycle [1|2|3|]::
+  number of levels in one-order for log-countour plot (default 3).
+
+:--slog_cycle [1|2|3|]:
+  number of levels in one-order for log-tone plot (default 3). 
+
+:--log:
+  approximately log-scaled levels for tone and contour plot.
+
+:--log_cycle [1|2|3|]::
+  number of levels in one-order for log- tone and countour plot (default 3).
 
 :--tone [a|e|f|b|c]:
   set tone subroutine:
@@ -223,10 +256,15 @@ if you want to know more description, please read
     b (=> DCL.uetonb is used)
     c (=> DCL.uetonc is used)
 
-:--nocolorbar:         
+:--colorbar_pos [b|l]:
+  set position of color bar
+    b (=> bottom (default)
+    l (=> left
+
+:--nocolorbar:
   do not draw color bar
 
-:--nozero:             
+:--nozero:
   do not draw zero contour
 
 :--udsfmt [strings]:
@@ -234,12 +272,12 @@ if you want to know more description, please read
 
 = EXAMPLES
 
-For a single GPhys variable, 
+For a single GPhys variable,
 
     % gpview data.nc@temp
     % gpview data.nc@temp,lon=130:150,lat=0:90:2
-    % gpview --nocont data.nc@temp,lon=130:150,lat=0    
-    % gpview --noshade data.nc@temp,lon=130:150,lat=0    
+    % gpview --nocont data.nc@temp,lon=130:150,lat=0
+    % gpview --noshade data.nc@temp,lon=130:150,lat=0
     % gpview --mean lon data.nc@temp,lon=130:150,lat=0
     % gpview --exch data.nc@temp,lon=130:150,lat=0
     % gpview --animate lon data.nc@temp,lon=130:150
@@ -255,12 +293,13 @@ For a single GPhys variable,
 Example to set DCL internal variable,
     % gpview -ud:lmsg=false --itr 5 --exch data.nc@temp,lon=0
 
-For multiple GPhys variables, 
+For multiple GPhys variables,
     % gpview data1.nc@temp data2.nc@temp
-    % gpview data1.nc@temp,lon=130:150 data2.nc@temp,lon=150:170
-    % gpview --overplot=3 --anim t data.nc@temp,lon=0 data.nc@temp,lon=10
-    % gpview --var temp,lon=130:150 data?.nc 
-    % gpview --anim t --var temp,lon=130:150 data*.nc 
+    % gpview data1.nc@temp,lon=130:150 data2.nc@temp,lon=150:170 
+    % gpview --overplot=3 data.nc@temp,lon=0 data.nc@temp,lon=10
+    % gpview --slice lon=130,lat=0 data1.nc@temp data2.nc@temp --overplot 2
+    % gpview --var temp,lon=130:150 data?.nc
+    % gpview --anim t --var temp,lon=130:150 data*.nc
 
 =end
 #
@@ -271,9 +310,9 @@ For multiple GPhys variables,
 #   2005/01/09  D Tsukahara (add option --animate, smooth, alternate, index )
 #   2005/01/10  D Tsukahara (transpose axis with attr "positive:down" ,
 #                            abailable loopsense_flag. )
-#   2005/01/10  D Tsukahara (implement GGraph::color_bar, and margin_info, 
+#   2005/01/10  D Tsukahara (implement GGraph::color_bar, and margin_info,
 #                            which file name, date, and toolname. )
-#   2005/01/11  D Tsukahara ( 1. write document about OPTIONS. 
+#   2005/01/11  D Tsukahara ( 1. write document about OPTIONS.
 #                             2. add many options. more info please see document. )
 #   2005/01/23  S Takehiro  (add option --similar, map_axis)
 #   2005/02/09  D Tsukahara && M Nakano (add option --reverse, --Gr, --map)
@@ -287,7 +326,7 @@ For multiple GPhys variables,
 #   2005/08/07  S Takehiro (option --overplot implemented)
 #   2005/08/09  S Takehiro && T Horinouchi (add option --help and help function.
 #                                             URL information corrected. )
-#   2005/08/23  S Takehiro (multiple Gphys variables plotting 
+#   2005/08/23  S Takehiro (multiple Gphys variables plotting
 #                            and option --var implemented)
 #   2005/08/23  S Takehiro (common methods to gp* command moved to gpcommon.rb)
 #   2005/10/30  S Takehiro (option --animate recoverd)
@@ -301,7 +340,7 @@ 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 
+#   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,
@@ -318,6 +357,11 @@ For multiple GPhys variables,
 #   2015/09/05  Y O Takahashi && S Takehiro (add option --delay)
 #   2015/09/06  Y O Takahashi && S Takehiro (defalt wsn number changed)
 #   2017/06/07  S Takehiro (int option bug fix)
+#   2021/02/22  S Takehiro (added options --log, --log_scale, --clog, --clog_scale, --slog, --slog_scale)
+#   2021/05/19  S Takehiro (added options --sum)
+#   2021/05/25  S Takehiro (added options --slice)
+#   2021/07/16  S Takehiro (added options --colorbar_pos)
+#   2021/07/16  S Takehiro (DCL colormap setting fixed)
 #
 #################################################
 require "getoptlong"        # for option_parse
@@ -359,7 +403,7 @@ def each_along_dims(gphys, loopdim)
   if loopdim.is_a?(String)
     dimname = loopdim
   elsif
-    if loopdim < 0 
+    if loopdim < 0
       dimname = gphys.coord(gphys.rank + loopdim).name
     else
       dimname = gphys.coord(loopdim).name
@@ -369,7 +413,7 @@ def each_along_dims(gphys, loopdim)
   end
 
   loopdim_na = gphys.coord(dimname).val                      # get coord ary
-  loopdim_na = loopdim_na[-1..0] if $OPT_reverse || $OPT_Gr  # reverse  
+  loopdim_na = loopdim_na[-1..0] if $OPT_reverse || $OPT_Gr  # reverse
   loopdim_na.each { |x|
     yield( gphys.cut(dimname=>x) )
   }
@@ -387,12 +431,12 @@ def draw_setup
   DCL.sgpset('lfull', true)               # use full area in the window
   DCL.sgpset('lfprop',true)               # use proportional font
   DCL.uscset('cyspos', 'B' )              # move unit y axis
-  
+
   # viewport size
   GGraph.set_fig('viewport'=>$VIEWPORT)
   GGraph.set_fig( 'itr'=>($OPT_itr == nil) ? 1 : $OPT_itr.to_i )
-  GGraph.set_fig("xrev"=>"units:mb,units:hPa,units:millibar,positive:down",  
-                 "yrev"=>"units:mb,units:hPa,units:millibar,positive:down") 
+  GGraph.set_fig("xrev"=>"units:mb,units:hPa,units:millibar,positive:down",
+                 "yrev"=>"units:mb,units:hPa,units:millibar,positive:down")
 
   # set options
   min_range,  max_range  = __split_range($OPT_range)
@@ -409,7 +453,7 @@ def draw_setup
     tone_interval = nil
   end
   GGraph.set_linear_contour_options(
-                                    'int' => cont_interval, 
+                                    'int' => cont_interval,
                                     'min' => ( min_crange  || min_range ),
                                     'max' => ( max_crange  || max_range )
                                     )
@@ -432,6 +476,19 @@ def draw_setup
     $OPT_patterns=$OPT_patterns.split(',').map!{|v| v.to_f }
   end
 
+  if ( $OPT_log ) 
+    $OPT_slog=$OPT_log;  $OPT_clog=$OPT_log
+  end
+
+  $OPT_slog_cycle = 3 if ( ! $OPT_slog_cycle ) 
+  $OPT_clog_cycle = 3 if ( ! $OPT_clog_cycle ) 
+
+  if ( $OPT_log_cycle ) 
+    $OPT_slog_cycle=$OPT_log_cycle;  $OPT_clog_cycle=$OPT_log_cycle
+  end
+  $OPT_slog_cycle=$OPT_slog_cycle.to_i;  $OPT_clog_cycle=$OPT_clog_cycle.to_i
+
+
   # similar projection
   if ($OPT_similar)
     if /([\d\-.]*),([\d\-.]*),([\d\-.]*)/ =~ $OPT_similar
@@ -471,7 +528,7 @@ def draw_setup
     GGraph::set_map(map_type=>true)
   end
 
-  # time axis 
+  # time axis
   if ($OPT_time_ax)
     $OPT_time_ax = false if $OPT_time_ax == "false"
     GGraph.set_axes('time_ax'=>$OPT_time_ax)
@@ -484,8 +541,8 @@ def draw(gp, draw_flag)
   # draw hontai
   case draw_flag
   when "line"
-    if ( $Overplot == 1 ) 
-      GGraph.line(gp, 
+    if ( $Overplot == 1 )
+      GGraph.line(gp,
                   true,
                   "title"=>$OPT_title,
                   "index"=>($OPT_index||1),
@@ -496,7 +553,7 @@ def draw(gp, draw_flag)
                   "max" => __split_range($OPT_range)[1]
                   )
     elsif ( $OPT_overplot_color || $OPT_Opc )
-      GGraph.line(gp, 
+      GGraph.line(gp,
                   false,
                   "title"=>$OPT_title,
                   "index"=>($OPT_index.to_i+$Overplot*10),
@@ -507,7 +564,7 @@ def draw(gp, draw_flag)
                   "max" => __split_range($OPT_range)[1]
                   )
     else
-      GGraph.line(gp, 
+      GGraph.line(gp,
                   false,
                   "title"=>$OPT_title,
                   "index"=>($OPT_index||1),
@@ -521,12 +578,12 @@ def draw(gp, draw_flag)
     if ( $Overplot < $Overplot_max )
       $Overplot += 1
     else
-      $Overplot = 1 
+      $Overplot = 1
     end
 
   when "mark"
-    if ( $Overplot == 1 ) 
-      GGraph.mark(gp, 
+    if ( $Overplot == 1 )
+      GGraph.mark(gp,
                   true,
                   "title"=>$OPT_title,
                   "index"=>($OPT_index||1),
@@ -537,7 +594,7 @@ def draw(gp, draw_flag)
                   "max" => __split_range($OPT_range)[1]
                   )
     else
-      GGraph.mark(gp, 
+      GGraph.mark(gp,
                   false,
                   "title"=>$OPT_title,
                   "index"=>($OPT_index||1),
@@ -551,17 +608,19 @@ def draw(gp, draw_flag)
     if ( $Overplot < $Overplot_max )
       $Overplot += 1
     else
-      $Overplot = 1 
+      $Overplot = 1
     end
 
   when "full"
     GGraph.tone(gp,
-                true, 
+                true,
                 "title"=>$OPT_title,
                 "annotate"=>$annotate,
                 "transpose"=>$OPT_exch,
                 "levels"=>$OPT_slevels,
                 "patterns"=>$OPT_patterns,
+                "log"=>$OPT_slog,
+                "log_cycle"=>$OPT_slog_cycle,
                 "auto"=>$auto,
                 "tonf"=>$tonf,
                 "tonb"=>$tonb,
@@ -569,23 +628,27 @@ def draw(gp, draw_flag)
                 "xcoord"=>$OPT_xcoord,
                 "ycoord"=>$OPT_ycoord
                 )
-    GGraph.contour(gp, 
+    GGraph.contour(gp,
 		   false,
 		   "transpose"=>$OPT_exch,
 		   "levels"=>$OPT_clevels,
 		   "nozero"=>$OPT_nozero,
+                   "log"=>$OPT_clog,
+                   "log_cycle"=>$OPT_clog_cycle,
                    "label"=>$OPT_label,
                    "xcoord"=>$OPT_xcoord,
                    "ycoord"=>$OPT_ycoord
 		   )
   when "nocont"
     GGraph.tone(gp,
-                true, 
+                true,
                 "title"=>$OPT_title,
                 "annotate"=>$annotate,
                 "transpose"=>$OPT_exch,
                 "levels"=>$OPT_slevels,
                 "patterns"=>$OPT_patterns,
+                "log"=>$OPT_slog,
+                "log_cycle"=>$OPT_slog_cycle,
                 "auto"=>$auto,
                 "tonf"=>$tonf,
                 "tonb"=>$tonb,
@@ -596,13 +659,15 @@ def draw(gp, draw_flag)
   when "noshade"
     mj = DCL.udpget('indxmj')
     mn = DCL.udpget('indxmn')
-    GGraph.contour(gp, 
-                   true, 
+    GGraph.contour(gp,
+                   true,
                    "title"    =>$OPT_title,
                    "annotate"=>$annotate,
 		   "transpose"=>$OPT_exch,
-                   "levels"=>$OPT_clevels, 
+                   "levels"=>$OPT_clevels,
                    "nozero"=>$OPT_nozero,
+                   "log"=>$OPT_clog,
+                   "log_cycle"=>$OPT_clog_cycle,
                    "label"=>$OPT_label,
                    "xcoord"=>$OPT_xcoord,
                    "ycoord"=>$OPT_ycoord
@@ -611,12 +676,20 @@ def draw(gp, draw_flag)
 
   # color bar
   if ( ( draw_flag == "full") || ( draw_flag == "nocont") ) && $colorbar
-    GGraph::color_bar(
-                      "left"      => true,
-                      "landscape" => true
-                      )
+    if $OPT_colorbar_pos == "l" then
+      GGraph::color_bar(
+                      "left"      => false,
+                      "landscape" => false,
+                      "log" => $OPT_slog
+      )
+    else
+      GGraph::color_bar(
+                       "left"      => true,
+                       "landscape" => true,
+                       "log" => $OPT_slog
+                       )
+    end
   end
-
 end
 
 
@@ -632,17 +705,17 @@ def set_vpsize( default_vp, aspect=2.0 )
   hlength =   x1 - x0
   vlength =   y1 - y0
   # center grid of viewport
-  cen_of_vp = [ x0 + hlength/2.0, y0 + vlength/2.0  ] 
+  cen_of_vp = [ x0 + hlength/2.0, y0 + vlength/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
   else
-    vlength = hlength / aspect     
+    vlength = hlength / aspect
     y0 = cen_of_vp[1] - vlength/2.0
     y1 = cen_of_vp[1] + vlength/2.0
-  end  
+  end
 
   return [ x0, x1, y0, y1 ]
 
@@ -666,7 +739,7 @@ def __split_range(range)
   else
     raise "invalid range: variable subset specification error. split range with ':'\n\n"
   end
-  
+
   return min, max
 end
 
@@ -691,7 +764,7 @@ def check_dclopts
     value = opt[4..-1].split('=')[1]
     dcl_set_params(pkg,name,value)
   }
-  
+
 end
 
 def dcl_set_params(pkg,name,value)
@@ -715,7 +788,7 @@ def dcl_set_params(pkg,name,value)
   else
     eval( "DCL.#{pkg}r#{set}(name,value.to_f)" )
   end
-  
+
 end
 
 #####################################################
@@ -728,6 +801,7 @@ check_dclopts
 parser = GetoptLong.new
 parser.set_options(
                    ###    global option   ###
+                   ['--slice',                    GetoptLong::REQUIRED_ARGUMENT],
                    ['--var',                      GetoptLong::REQUIRED_ARGUMENT],
                    ['--wsn',                      GetoptLong::REQUIRED_ARGUMENT],
                    ['--clrmap',                   GetoptLong::REQUIRED_ARGUMENT],
@@ -753,10 +827,11 @@ parser.set_options(
                    ['--reverse',                  GetoptLong::NO_ARGUMENT],
                    ['--Gr',                       GetoptLong::NO_ARGUMENT],
                    ['--mean',                     GetoptLong::REQUIRED_ARGUMENT],
+                   ['--sum',                     GetoptLong::REQUIRED_ARGUMENT],
                    ['--stddev',                   GetoptLong::REQUIRED_ARGUMENT
 ],
                    ['--eddy',                     GetoptLong::REQUIRED_ARGUMENT],
-                   ['--diff',                     GetoptLong::NO_ARGUMENT],                   
+                   ['--diff',                     GetoptLong::NO_ARGUMENT],
                    ['--map',                      GetoptLong::REQUIRED_ARGUMENT],
                    ['--m',                        GetoptLong::NO_ARGUMENT],
                    ['--operation',                GetoptLong::REQUIRED_ARGUMENT],
@@ -783,8 +858,15 @@ parser.set_options(
                    ['--clevels',                  GetoptLong::REQUIRED_ARGUMENT],
                    ['--slevels',                  GetoptLong::REQUIRED_ARGUMENT],
                    ['--patterns',                 GetoptLong::REQUIRED_ARGUMENT],
+                   ['--log',                      GetoptLong::NO_ARGUMENT],
+                   ['--clog',                     GetoptLong::NO_ARGUMENT],
+                   ['--slog',                     GetoptLong::NO_ARGUMENT],
+                   ['--log_cycle',                GetoptLong::REQUIRED_ARGUMENT],
+                   ['--clog_cycle',               GetoptLong::REQUIRED_ARGUMENT],
+                   ['--slog_cycle',               GetoptLong::REQUIRED_ARGUMENT],
                    ['--tone',                     GetoptLong::REQUIRED_ARGUMENT],
                    ['--udsfmt',                   GetoptLong::REQUIRED_ARGUMENT],
+                   ['--colorbar_pos',             GetoptLong::REQUIRED_ARGUMENT],
                    ['--nocolorbar',               GetoptLong::NO_ARGUMENT],
                    ['--nozero',                   GetoptLong::NO_ARGUMENT],
                    ['--help',                     GetoptLong::NO_ARGUMENT]
@@ -801,7 +883,7 @@ end
 
 ## Print out help message
 if ($OPT_help)
-  help 
+  help
   exit(1)
 end
 
@@ -811,7 +893,7 @@ DCL::swlset('lwait', false) if ($OPT_nowait    || $OPT_Gw || $OPT_smooth || $OPT
 DCL::swlset('lalt',  true)  if ($OPT_alternate || $OPT_Ga || $OPT_smooth || $OPT_Gaw)
                                            # set backing store option
 if ($OPT_noannotate)
-  $annotate = false 
+  $annotate = false
 else
   $annotate = true
 end
@@ -820,12 +902,12 @@ $Overplot_max = ( $OPT_overplot.to_i || 1 )
 $Overplot = 1
 
 if ($OPT_nocolorbar)
-  $colorbar = false 
+  $colorbar = false
 else
   $colorbar = true
 end
 
-if ($OPT_tone) 
+if ($OPT_tone)
   case $OPT_tone
   when "a"
     $auto, $tonf, $tonb, $tonc = true,false,false,false
@@ -855,11 +937,9 @@ udsfmt = ($OPT_udsfmt||DCL.udqfmt())
 DCL.udsfmt(udsfmt)
 
 ## draw figure
-loopdim   = ( $OPT_animate || $OPT_anim ) 
+loopdim   = ( $OPT_animate || $OPT_anim )
 loopdim = loopdim.to_i if loopdim.to_i.to_s == loopdim
 
-### set colormap
-DCL.sgscmn($OPT_clrmap||1)
 ## open work station
 if NumRu::DCL::DCLVERSION.to_f >= 6
   DCL.gropn($OPT_wsn||1)
@@ -867,15 +947,19 @@ else
   DCL.gropn($OPT_wsn||4)
 end
 
+### set colormap
+DCL.sgscmn($OPT_clrmap||1)
+
 ## open netcdf variables
 
 while ARGV[0] do
   gturl = ARGV[0]
   gturl = gturl+'@'+$OPT_var if $OPT_var
+  gturl = gturl+','+$OPT_slice if $OPT_slice
 
   gp = GPhys::IO.open_gturl(gturl)
   print "  Reading #{gturl}\n"
-  ARGV.shift 
+  ARGV.shift
 
   ## for case of calculating difference of two gphys object
   if ($OPT_diff)
@@ -887,15 +971,21 @@ while ARGV[0] do
     gp = gp - GPhys::IO.open_gturl(gturl)
     print "  Taking difference: #{prev_gturl} - #{gturl}\n"
     gturl = prev_gturl+" - "+gturl
-    ARGV.shift 
+    ARGV.shift
   end
-  
+
   ## mean along any axis
   if ($OPT_mean)
     dims_mean = ($OPT_mean).split(/\s*,\s*/)
     dims_mean = dims_mean.map{|dim| dim.to_i.to_s == dim ? dim.to_i : dim}
   end
 
+  ## sum up along any axis
+  if ($OPT_sum)
+    dims_sum = ($OPT_sum).split(/\s*,\s*/)
+    dims_sum = dims_sum.map{|dim| dim.to_i.to_s == dim ? dim.to_i : dim}
+  end
+
   ## standard deviation along any axis
   if ($OPT_stddev)
     dims_stddev = ($OPT_stddev).split(/\s*,\s*/)
@@ -921,6 +1011,13 @@ while ARGV[0] do
         }
       end
 
+      ## sum up along any axis
+      if ($OPT_sum)
+        dims_sum.each{|dim|
+          g = g.sum(dim)
+        }
+      end
+
       ## standard deviation along any axis
       if ($OPT_stddev)
         dims_stddev.each{|dim|
@@ -947,14 +1044,14 @@ while ARGV[0] do
           kind_of_fig = "mark"
         elsif ($OPT_line || gp_rank == 1)
           kind_of_fig = "line"
-        elsif (!$OPT_line && gp_rank >= 2) && !$OPT_noshade && $OPT_nocont 
+        elsif (!$OPT_line && gp_rank >= 2) && !$OPT_noshade && $OPT_nocont
           kind_of_fig = "nocont"
-        elsif (!$OPT_line && gp_rank >= 2) && $OPT_noshade && !$OPT_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  
+      end
 
       draw(g, kind_of_fig)
   end