[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [dennou-ruby:001827] Re: GPhys::EP_Flux
堀之内様
塚原です.
> 「流れ」ベクトル用のメソッドを作りました。ggraph.rb の DCLExtモ
> ジュールに入れてあります(CVSのみ)。ugvect の代りに使うものですが、
> 前に話した通り、factor で全体の長さを微調整できます。
早速作ってくださったのですね. ありがとうございます.
# というか僕が作るはずだったのにすみません. 仕事が遅くて...
> さらに言えば、itr=1 でないと、ちゃんと流れになりません。
> 将来的には ugvect を使わないで地図投影に対応するものも欲しいです
> ね。
そうですね. 上記の問題以外にも ugvect は不便なことが多いので,
あるといいなと思います.
> 以下、テストプログラムです。U座標の縦横比にかかわらず、流れが
> 追えるのがわかると思います。
早速走らせて見ました. ugvect だと x 成分が強く出てしまうところ,
ug_flow_vect では望ましい方向を向いています. 素晴らしいです.
--------------------------------------
北海道大学院理学研究科 地球惑星科学専攻
地球流体力学研究室 M2 塚原大輔
email::daktu32@xxxxxx
--------------------------------------
> 塚原様：
>
> 「流れ」ベクトル用のメソッドを作りました。ggraph.rb の DCLExtモ
> ジュールに入れてあります(CVSのみ)。ugvect の代りに使うものですが、
> 前に話した通り、factor で全体の長さを微調整できます。
>
> また、xintv, yintv で指定する間引き機能をつけました。
> ugvect を呼んでますので、矩形座標専用です(ugvect)。
> さらに言えば、itr=1 でないと、ちゃんと流れになりません。
> 将来的には ugvect を使わないで地図投影に対応するものも欲しいです
> ね。
>
> factor=1.0は x, y 成分の長さの最大が平均グリッド間隔になるよう
> にという ugvect 的仕様です。流れを表すため、ベクトルの縦横比を、
> （U 座標の縦横比／V座標の縦横比）で調整します。速度場やフラック
> スの表示に好適。
>
> def ug_flow_vect( fx, fy, factor=1.0, xintv=1, yintv=1)
> raise ArgumentError,"Expect 2D arrays" if fx.rank != 2 || fy.rank != 2
> raise ArgumentError,"fx.shape != fy.shape" if fx.shape != fy.shape
> raise ArgumentError,"xintv must be a positive integer" if xintv < 0
> raise ArgumentError,"yintv must be a positive integer" if yintv < 0
> nx, ny = fx.shape
> if xintv >= 2
> idx = NArray.int(nx/xintv).indgen!*xintv # [0,xintv,2*xintv,..]
> fx = fx[idx, true]
> fy = fy[idx, true]
> end
> if yintv >= 2
> idx = NArray.int(ny/yintv).indgen!*yintv # [0,yintv,2*yintv,..]
> fx = fx[true, idx]
> fy = fy[true, idx]
> end
> nx, ny = fx.shape # again, because of xintv & yintv
> vx0,vx1,vy0,vy1 = DCL.sgqvpt
> ux0,ux1,uy0,uy1 = DCL.sgqwnd
> dvx = (vx1-vx0)/nx
> dvy = (vy1-vy0)/ny
> ax = (vx1-vx0)/(ux1-ux0) # factor to convert from U to V coordinate
> ay = (vy1-vy0)/(uy1-uy0) # factor to convert from U to V coordinate
> fxmx = fx.abs.max
> fymx = fy.abs.max
> raise "fx has no data or all zero" if fxmx == 0
> raise "fy has no data or all zero" if fymx == 0
> cn = [ dvx/(ax*fxmx), dvy/(ay*fymx) ].min # normarization constant
> lnrmal = DCL.ugpget('LNRMAL') # save
> DCL.ugpset('LNRMAL', false)
> DCL.ugvect( factor*cn*ax*fx, factor*cn*ay*fy )
> DCL.ugpset('LNRMAL', lnrmal) # recover
> end
>
> 以下、テストプログラムです。U座標の縦横比にかかわらず、流れが
> 追えるのがわかると思います。
>
> ---
> require "numru/ggraph"
> include NumRu
>
> nx = 20
> ny = 10
> ux = 10.0
> uy = 1.0
> fx = ux*NArray.sfloat(nx,ny).indgen!
> fy = uy*NArray.sfloat(nx,ny).indgen!
>
> DCL.gropn(1)
>
> DCL.grfrm
> DCL.grswnd(0,ux,0,uy)
> DCL.grsvpt(0.2,0.8,0.2,0.8)
> DCL.grstrn(1)
> DCL.grstrf
> DCL.usdaxs
> DCLExt.ug_flow_vect(fx, fy, 2.0)
>
> DCL.grfrm
> DCL.grswnd(0,ux,0,uy)
> DCL.grsvpt(0.2,0.8,0.2,0.8)
> DCL.grstrn(1)
> DCL.grstrf
> DCL.usdaxs
> DCLExt.ug_flow_vect(fx, fy, 1.0, 2, 1)
>
> DCL.grcls
> ---
>
> 堀之内
>