[dennou-ruby:003212] gpview へのパッチ

[Takayuki Tanigawa <tanigawa@xxxxxxxxxxxxxxxxxxxx>]

北大低温研の谷川と申します。

gpview に添付のパッチの様なものを取りこんでもらうことは
可能でしょうか？

具体的には、オプションに数学関数を指定すると、その演算を
施した結果を作図する、というものです。

使用例：
    % gpview --operation log10 data.nc@temp,lon=0

個人的には log10 しか使ってないのですが、数学関数のうち
引数が一つのもの（例：sqrt, sin など）は全て問題なく行えると
思います。
オプションの名前 (operation) は他の方が良いかもしれません。

ご検討頂けたら幸いです。

-- 
Takayuki Tanigawa
Center for Planetary Science / ILTS, Hokkaido University

--- gpview.org	2009-01-16 05:57:04.000000000 +0900
+++ gpview	2010-11-15 18:51:31.000000000 +0900
@@ -90,7 +90,10 @@
                            coast_world, border_world, plate_world, state_usa,
                            coast_japan, pref_japan  
 
-
+  --operation <math_func>: operation of the specified math function on
+                           the data.
+                           <math_func> should be a math function with one
+                           argument such as log10, sqrt, sin, etc.
 
 ===LINE/MARK OPTIONS
 
@@ -162,6 +165,7 @@
     % gpview --levels 0,10,20,30 --pattern 30999,40999,50999,60999,70999 data.nc@temp
     % gpview --overplot=3 --anim t data.nc@temp,lon=0
     % gpview --mark data.nc@temp,lon=0,lat=0
+    % gpview --operation log10 data.nc@temp,lon=0
 
   For multiple Gphys variables, 
     % gpview data1.nc@temp data2.nc@temp
@@ -553,6 +557,7 @@
                    ['--mean',                     GetoptLong::REQUIRED_ARGUMENT],
                    ['--map',                      GetoptLong::REQUIRED_ARGUMENT],
                    ['--m',                        GetoptLong::NO_ARGUMENT],
+                   ['--operation',                GetoptLong::REQUIRED_ARGUMENT],
                    ###     line/mark option    ###
                    ['--line',                     GetoptLong::NO_ARGUMENT],
                    ['--mark',                     GetoptLong::NO_ARGUMENT],
@@ -644,6 +649,11 @@
     }
   end
 
+## operation of a mathematical function
+  if ($OPT_operation)
+    eval "gp = gp.#{$OPT_operation}"
+  end
+  
   GGraph.margin_info($0, gturl) if $annotate  # draw margin infomation
   kind_of_fig = draw_setup(gp)                # determine figure kind