Re: [jmadata-users33] software to read GPV (and to convert intoNetCDF)

[Shin-ichi Suzuki <ssuzuki@xxxxxx>]

堀之内様
(Bjmadata-users の皆様
(B
(B防災科研の鈴木です。
(B
(B> > まず、鉛直方向の関数が p になっているのを level に変えてみました。
(B> COARDS では圧力の座標名は level でなくてもいいのではないでしょうか。
(B
(Bそうなんですね。手元にあった netcdf のデータが NCEP のやつだったもので、
(B単純にそれを真似してみました。
(B
(B> 時間はちゃんと書き出されるけど、その内容が「1年1月1日00Z から3時
(B> 間おき」なのですね（1年1月1日からでなく各年の1月1日から？）。そ
(B> の場合の単位は、"days since 2002-01-01 00:00+00" などとして 
(B> since 以下でオフセットを表すことになります
(B> (http://www.unidata.ucar.edu/packages/udunits/)。単位はもしかし
(B> て hours でしょうか。1年の始めからとなると、days に変えた方が分
(B> かりやすいですね。
(B
(B実は、投稿した変更だと、GrADS では時間軸方向にはちゃんと
(B表示できていませんでした。未確認のまま投稿してしまいました。
(B堀之内さんのおっしゃる通りなのですが、その後、gpv2nc.rb の時間の部分は
(B--------------------------------------------------------------ここから
(B    output.def_var('time', "sfloat", [dimt])
(B    output.var('time').put_att('units', 'hours since 2002-08-30 18:00:0.0')
(B    output.var('time').put_att('long_name', 'Time')
(B    output.var('time').put_att('actual_range' , '0, 18')
(B--------------------------------------------------------------ここまで
(Bというようにしました。これで作成された netcdf ファイルを ncdump で
(B見ると、時間に関する定義の部分は
(B--------------------------------------------------------------ここから
(B        float time(time) ;
(B                time:units = "hours since 2002-08-31 18:00:0.0" ;
(B                time:long_name = "Time" ;
(B                time:actual_range = "0, 18" ;
(B--------------------------------------------------------------ここまで
(Bというようになり、データの部分も
(B--------------------------------------------------------------ここから
(B time = 0, 3, 6, 9, 12, 15, 18 ;
(B--------------------------------------------------------------ここまで
(Bでしたので（これはもちろん最初からちゃんと出ていたのですが）、
(Bうまく GrADS で見えるようになりました。
(B
(Bもともとのプログラムの場合、gpv2nc.rb は /usr/bin/gpv2nv として
(Bインストールされますが、今のところ、gpv2nc.rb を自分のディレクトリに
(Bコピーして、ファイル毎に時間の部分を書き換え、
(B% ruby ./gpv2nc.rb -o 2002AUG31.MSM18PLM018_2.nc MSM18PLM018_2
(Bというように使っています。
(B
(B> > ＃ もちろん、直接改良してくれればありがたいですが……
(B> そのつもりですが、少し時間がかかりますので、もしこれまでのご変更
(B> で不都合がなければしばらくそれをお使い下さい。
(B
(Bお手数おかけします。ありがとうございます。
(Bあまり急いではおりません。
(B
(B--------------------------------------------------
(B 鈴木 真一 ／ 防災科学技術研究所 総合防災研究部門
(B e-mail : ssuzuki@xxxxxx  
(B--------------------------------------------------