竹広です.
GPhys レシピを参考に, 鉛直座標を sigma から高度に変換する計算を
手習いがてらいろいろ試しております. その過程でちょっと気になった
ことをご報告/質問させてください.
気になったことは, 4 次元データ [lon,lat,level,time] の鉛直座標をいっぺんに
変換させてから最後尾の時間切出しを行うのと, 先に時間だけ切り出しておいて
その時間断面だけ鉛直座標変換した結果がどうも一致しないようなのです.
その原因を追求する実力がないので識者にお伺いする次第です.
レシピの interpolate4_theta_coord.rb を修正してそのような例を
作ってみました. あやふやな理解で作成しているので意図していることと
違ったことをやっているかも知れません. そのようなところがありましたら
ご教示くだされば幸いです.
Takepiro(竹広真一)＠数理解析研究所. 京都大学
E-mail:takepiro@xxxxxxxxxxxxxx
takepiro@xxxxxxxxxxxxxxxxxxxx
Attachment:
interpo4_theta_coord_test.rb
Description: Binary data