Sunday, 26 January 2025

氷河合宿その4i:netcdfファイルが壊れてる?

 ERA5データのダウンロードに関して、トラブルがあるとの相談があり、対処法がわかりましたので共有します。

昨年末まで順調に稼働していたダウンロードスクリプトですが、今年に入って「ダウンロードしたnetcdfファイルが開けない」との相談が届きました。こちらで試したところ、同様の症状が発生することを確認しました。なお、pressure levelのファイルは問題なく開けます。

スクリプトからではなく、ECMWFのウェブサイトから直接ダウンロードしてみたところ、圧縮ファイルには二つのファイルが含まれており、それぞれ "u10, v10, t2m, d2m, ap" と "tp, ssrd, strd" が含まれていました。後者の要素は再解析データの中では「フラックス」として扱われているため、この混在が問題なのかもしれない、と考えました。

そこで、"u10, v10, t2m, d2m, ap" と "tp, ssrd, strd" を別々にダウンロードするスクリプトを用意し、ダウンロードしてみたところ、それぞれのファイルを問題なく(Pythonで)読み込むことができました。

ECMWFに何とかして欲しいところですが、とりあえずの対処療法です。欲しいポイントのデータを抽出するスクリプトを書き換える必要がありますが、どうしてもわからない人は相談に乗ります。