はじめに

全球の気象データを提供している NOAA/GFS について、少し調べたので備忘録として記録しておく。専門外の人間が調べた結果を書いているので、もし理解に間違いがあれば是非コメント下さい！

予報データの入手先（2か所）

データをダウンロードできるページは二つあり、ファイル名が違うものの中身は全く一緒である。

https://nomads.ncep.noaa.gov/：過去10日間のデータが置いてあり更新は早い
https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs：過去半年前までのデータがあるが更新は遅め

1番のサイトについて

データのダウンロード方法

こちらのページ(https://nomads.ncep.noaa.gov/)の中段のテーブルの Global Models のGFSと名のついたものを選べばよい。分解能は 0.25 Degree, 0.50 Degree, 1.00 Degree の三種類ある。表の右側の「grib filter」、「https」、「OpenDAP」からダウンロードできる。

以下、https を選択した場合、モデル名＋年月日のディレクトリ(gfs.YYYYMMDD/)が表示される。それをクリックすると、予報の発表時刻(00, 06, 12, 18 の四種類)が表示される。それをクリックすると、atom(大気のデータ)、wave(海象のデータ)が表示される。それをクリックすると大量のファイルが表示される。

ファイル名の意味

ファイル名は

res：分解能（0p25, 0p50, 1p00の三種類）
h：予報を発表した時
f：予報を発表した時刻から何時間後

とすると命名のルールは、

gfs.t<hh>z.pgrb2.<res>.f<fff>

となる。つまり、2021年3月26日の12時に発表した、分解能1度の9時間後の予報データは、gfs.20210326/atom/フォルダ内の

gfs.t12z.pgrb2.1p00.f009

となる。予報は6時間ごとに発表され、予報時間は0時間後から384時間後まで3時間ごとにデータがある状況である。

拡張子が無いデータ：実際の予報データが格納されているファイル
拡張子が.idx：予報項目が含まれているヘッダファイル

である。

2番のサイトについて

データのダウンロード方法

こちらのページ(https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs)の中段の「GFS Forecast」からデータをダウンロードできる。モデルは0.5度の分解能のもの（004）と、1度の分解能のもの（003）の二種類がある。表の右側の「Data Access Links」からHTTPS, TDS, AIRS の好きな方法でダウンロードできる。

以下、HTTPS を選択した場合、年月のディレクトリ（YYYYMM/）が表示される。それをクリックすると次は年月日のディレクトリ（YYYYMMDD/）が表示される。それをクリックすると、大量のファイルが表示される。

※因みに「GFS Analysis」は再解析データなので、今回は割愛する。

ファイル名の意味

ファイル名は

mode：3=1度の予報データ、4=0.5度の用法データ
Y：予報を発表した年
M：予報を発表した月
D：予報を発表した日
h：予報を発表した時
m：予報を発表した分
f：予報を発表した時刻から何時間後

とすると、命名ルールは、

gfs_<mode><YYYYMMDD><hhmm>_<fff>.grb2

となる。つまり、2021年3月26日の12時に発表した、分解能1度の9時間後の予報データは、

gfs_3_20210329_1200_009.grb2

となる。予報は6時間ごとに発表され、予報時間は0時間後から384時間後まで3時間ごとにデータがある状況である。因みに、

拡張子が.grb2：実際の予報データが格納されているファイル
拡張子が.grb2.inv：予報項目が含まれているヘッダファイル

である。

データの項目について

含まれている予報値（データ項目）は、ファイルごとによって変わる。例えば、gfs.t12z.pgrb2.1p00.f009 (gfs_3_20210329_1200_009.grb2) の中身を wgrib2 で確認すると、

$ wgrib2 gfs.t12z.pgrb2.1p00.f009
1:0:d=2021032612:PRMSL:mean sea level:9 hour fcst:
2:75246:d=2021032612:CLMR:1 hybrid level:9 hour fcst:
3:86296:d=2021032612:ICMR:1 hybrid level:9 hour fcst:
---中略---
741:45317259:d=2021032612:HGT:PV=-2e-06 (Km^2/kg/s) surface:9 hour fcst:
742:45396668:d=2021032612:PRES:PV=-2e-06 (Km^2/kg/s) surface:9 hour fcst:
743:45473804:d=2021032612:VWSH:PV=-2e-06 (Km^2/kg/s) surface:9 hour fcst:

743 項目のデータがあることが確認できる。データ項目の意味は、https://www.nco.ncep.noaa.gov/pmb/products/gfs/の各項目のページに書いてある。例えば、gfs.t12z.pgrb2.1p00.f009 (gfs_3_20210329_1200_009.grb2) の場合は、https://www.nco.ncep.noaa.gov/pmb/products/gfs/gfs.t00z.pgrb2.1p00.f003.shtmlに書いてある。

おわりに

これでGFSの予報データから必要な情報を取り出して扱うことができる。今のところ他の記事では地表面の気温（"TMP:surface"）しか使用してないが、他の項目（例えば、雲のデータ、降水量、気圧etc...）も調べて使っていきたい。