[dennou-ruby:004053] [gphys] ganalysis/eof.rb の修正

[Tsuyoshi Koshiro <tkoshiro@xxxxxxxxxxxxx>]

神代です。

gphys-1.5.4 での numru/ganalysis/eof.rb の変更で、対象のGPhys
オブジェクトが2次元の場合が想定されなくなっているようです。
このため、サンプルプログラム(eof.rbを実行)のような場合が動かな
くなっています。
次のようにするとよいのではないかと思います。
(最近のDCLに合わせて、装置番号を4から1に変えました)

--- eof.rb.old	2018-08-05 11:54:39.627147120 +0900
+++ eof.rb	2018-08-06 09:39:19.136014114 +0900
@@ -93,6 +93,7 @@
         w /= w.mean
         w.reshape!(n)
       else
+        if new_grid.rank > 1
         nc0 = new_grid.coord(0)
         nc1 = new_grid.coord(1)
         if !(opts[:disable_weight]||opts["disable_weight"]) &&
@@ -107,6 +108,9 @@
         else
           w = nil
         end
+        else
+          w = nil
+        end
       end
 
       ary = NArrayMiss.new(gphys.typecode, n_lost, n)
@@ -260,7 +264,7 @@
   eof,rate = gphys.eof("t")
 
   max = 5
-  DCL::gropn(4)
+  DCL::gropn(1)
   DCL::grfrm
   DCL::grsvpt(0.1,0.9,0.1,0.9)
   DCL::grswnd(-max,max,-max,max)


--
神代 剛 (こうしろ つよし)
気象研究所気候研究部第一研究室 リサーチアソシエイト