[dennou-ruby:001869] Re: GPhys::EP_Flux ver.0.0.2.2 released

[Takeshi Horinouchi <horinout@xxxxxxxxxxxxxxxxxx>]

塚原さま
(cc: 西澤さま)

> GPhys::EP_Flux の 0.0.2.2 をリリースしました. 間隔おかずにすみません.
> 
>   http://ruby.gfd-dennou.org/products/ep_flux/ep_flux.0.0.2.2.tar.gz

デモで下記のエラーが出ます。原因は NArrayMiss は NArray を知って
るけど、その逆は知らないので、 NArray * NArrayMiss が出来ないと
いうことですね。従って現時点での対症療法は掛け算の順序を逆にする
です。現在の左側の項は軸起源なので、必ず NArray ですので（欠損な
し）。

さて、西澤さま。narray_miss.rb 内で NArray の coerce を再定義し
て、相手が NArrayMiss ならよろしくやってくれて、それ以外なら 
オリジナルのに頼むようにするというのはどうでしょう。

% ruby demo_NCEP_1.rb 
/usr/local/lib/ruby/site_ruby/1.8/numru/derivative.rb:143:in `*': NArrayMiss#to_f should return Float (TypeError)
        from /usr/local/lib/ruby/site_ruby/1.8/numru/derivative.rb:143:in `threepoint_O2nd_deriv'
        from /usr/local/lib/ruby/site_ruby/1.8/numru/gphys/derivative.rb:88:in `threepoint_O2nd_deriv'
        from /usr/local/lib/ruby/site_ruby/1.8/numru/gphys/ep_flux.rb:331
        from /usr/local/lib/ruby/site_ruby/1.8/numru/gphys/ep_flux.rb:330:in `call'
        from /usr/local/lib/ruby/site_ruby/1.8/numru/gphys/ep_flux.rb:386:in `deriv'
        from /usr/local/lib/ruby/site_ruby/1.8/numru/gphys/ep_flux.rb:424:in `ep_full_sphere'
        from demo_NCEP_1.rb:14

参考まで、依存ライブラリーのバージョンは numru-derivative.0.1.1,
narray_miss-1.1.0, gphys は CVS 版です。（例の8/26の
each_along_dims_write のデバグ以降のコミットはありません）


>   * ((<GGraph::divergence_on_merdional_section>)) を追加
>                                 (in ggraph_on_merdional_section.rb)
>     * EP Flux divergence が目で見てわかるようなベクトルを描画

名前について。divergence_on_merdional_section では、ベクトルでな
く divergence 自体の描画に見えます。実体は sinφ を横軸にしたベ
クトル図ですよね？ なお、横軸を sinφ にするなら、軸もそれで書か
ないと駄目です。今の場合、目盛と実際の位置がずれてるでしょ。

demo3 について。lcntl を false にするかわりに、下付き添字指定の
制御文字を '_' 以外の何かにするという手があります。縦軸はその方
が綺麗になりますね。

堀之内