気象データは GRIB 形式で得られることが多い。これを R で直接読み込む方法があればいいのだが、きっと無い。
google の検索結果 そこで、wgrib でデコード(符号化されたものを解くこと)することが必要になる。ちなみに、R は netCDF は読めるようになっている。
wgrib の使い方については、他のサイトにも書いてあるが、とりあえず、wgrib をコマンドラインのオプション無しで起動したときのメッセージを掲載しておく。
Portable Grib decoder for NCEP/NCAR Reanalysis etc.
it slices, dices v1.8.0.3c beta (1-17-03) Wesley Ebisuzaki
DWD-tables 2,201-203 (9-13-2001) Helmut P. Frank
spectral: Luis Kornblueh (MPI)
usage: wgrib [grib file] [options]
Inventory/diagnostic-output selections
-s/-v short/verbose inventory
-V diagnostic output (not inventory)
(none) regular inventory
Options
-PDS/-PDS10 print PDS in hex/decimal
-GDS/-GDS10 print GDS in hex/decimal
-verf print forecast verification time
-ncep_opn/-ncep_rean default T62 NCEP grib table
-4yr print year using 4 digits
Decoding GRIB selection
-d [record number|all] decode record number
-p [byte position] decode record at byte position
-i decode controlled by stdin (inventory list)
(none) no decoding
Options
-text/-ieee/-grib/-bin convert to text/ieee/grib/bin (default)
-nh/-h output will have no headers/headers (default)
-H output will include PDS and GDS (-bin/-ieee only)
-append append to output file
-o [file] output file name, 'dump' is default
このうち、Inventory の部分は、NCEP/NCARのデータを grib 形式のまま使うための補助的なファイルの作成に使うと考えられるので、(おそらく)関係ない。Decodeing GRIB の部分が関係ある。それほど量は多くないので、丁寧に見て欲しい。-d は「record number」を指定しなさい、とあるが、all でよかろう。(none) の部分は、オプションを何も書かないと decode しないことを説明している。