はじめに

この記事ではGTOOL3形式のファイルを(主にPythonで)扱うのに有用なオープンになっている情報をまとめる。

GTOOL3形式は沼口淳氏により策定された格子点データフォーマットである。¹現在では主に気候モデルMIROCの入出力に用いられており、関連する情報はMIROCグループ関係者限定のMIROC wikiにて共有されている。しかし最近ではMIROCと結合する陸面モデルMATSIRO/ILSや、MIROCの海洋コンポーネントCOCOと結合した非静力学モデルNICAM-COCO(NICOCO)など、MIROCグループ外にも利用が広まったことから、GTOOL3のデータ解析・可視化のための公開情報の需要が高まっている。

また、近年では大気海洋のデータ解析や可視化をPythonで行うようになり、Fortranで解析しGrADSで描画することを前提に開発されたツール群では不十分であると感じられるようになった。そのため、GTOOL3を読み込むためのPythonパッケージが複数の開発者により公開されている。

そこで本記事では、2024年11月現在、GTOOL3に関して公開されているツールや記事について、私の知っている情報をまとめて整理する。GTOOL3を初めて扱う初学者の道標とすることを目的としており、個別のツールの使い方の詳細については深くは立ち入らない。

なお、gtool4はNetCDFの規約として策定されたもので、GTOOL3とは異なるデータ形式である。そのため、gtool4についてはこの記事では扱わない。

データ形式

概要

GTOOL3のオリジナルのドキュメント「格子点データ解析ツールGTOOL3利用の手引き」は以下のリンクからアクセスできる。

GTOOL3のデータ形式は、ヘッダ部とデータ部の2つから構成される。ヘッダ部は16文字×64レコード、データ部は単純バイナリで格納される。

現在GTOOL3形式として流通しているデータは、ビックエンディアンとリトルエンディアンのどちらの場合もある。MIROCなどのモデル出力ではビックエンディアンが使われていることが多い。

通常はFortranのsequential accessで出力されヘッダ部とデータ部の前後にデータ長を表す4byteのmarkerが付くが、sequential accessでmarkerがつかない場合がある²。これは扱うデータサイズが大きくなり、配列の長さがsequential access のmarkerの4byte整数で表現できる長さを超えることに対処するためである。

軸情報ファイル

GTOOL3ではヘッダ部に緯度・経度・高度軸を表す文字列が書かれているが、それがどのような数値の配列であるかは軸情報ファイルという別のGTOOL3形式のファイルに記述されている。そのため、読み出したいファイルとは別に軸情報ファイルを用意しておく必要がある。

共用の計算機などですでにある場合はそのまま利用すれば良いが、ない場合は以下のリポジトリからダウンロードすることができる。³

利用するツールに応じて、環境変数GTAXDIRやGTAX_PATH等に軸情報ファイルのあるディレクトリを設定する必要がある。

なお、軸情報ファイルのうちGLON*, GGLA*, GLAT*については後述するngtmkaxコマンドを用いて作成することができる。

コマンドラインツール

gtool3ライブラリ (`gt*`コマンド)

沼口敦氏が作成したオリジナルのGTOOL3。ソースは公開されていないと思われる。

竹広真一氏がdcl-5.0に対応させたGTOOL3-DCL5が地球流体電脳倶楽部で公開されている。ライブラリの実態はほとんどオリジナルと変わっていないとのこと。

ドキュメントは以下にある。

libgtool3 (`ngt*`コマンド)

ngtから始まるコマンドラインツールが、GTOOL3ファイルの中身を確認したり、簡単な加工を行うのに有用である。上記gt*コマンドで対応されていない数値フォーマットにも対応している。

インストールの手順は以下の通り。

# 適当なディレクトリにクローンして中に入る
git clone https://github.com/kjoti/libgtool3.git
cd libgtool3

# ビルド & インストール
./configure
make
sudo make install

環境変数PATHに/usr/local/bin、LD_LIBRARY_PATHに/usr/local/libが含まれているか確認し、なければ追記する。

`ngtcf` `ncgt` (netCDF変換)

netCDFとの変換コマンドとしてngtcfやncgtがあるが、これらは上記リポジトリには含まれておらず、作成者や入手方法は把握できていないので情報求む。

extgtool3 (Extra-gtool3)

山下陽介氏が作成しているGTOOL3形式のデータを扱うためのコマンドラインツール集。オリジナルの機能拡張を目的としている。

TOUZA (楊座)

齋藤冬樹氏が開発しているDCLに依存しないGTOOL3のデータ加工ツール。

Pythonでの読み込み

Pythonでの読み込みパッケージは、私が確認している限り3つ存在する。公開された時系列順に以下に記す。

gtool io submodule

Kim Hyungjun氏作成。

解説記事はこちら。

xgtool3

橋本恵一氏作成。Xarrayで読み出せる。

TOUZA/Zbt

齋藤冬樹氏が開発している、Xarray engineとして動作するGTOOL3バックエンド。

終わりに

記事の修正やリンクの追加等の要望を歓迎します。

https://www.metsoc.jp/tenki/pdf/2001/2001_08_0577.pdf ↩
Fortranのバイナリ出力についての一般的な話は以下の記事を参考にされたい。 https://qiita.com/convection/items/f3be09bb280df030f980 ↩
この軸情報ファイルは私が勝手に公開しているものです。著作権はないと思いますがグレーですし、MIROCの公式で管理して頂けるならそちらに移管したいと考えています。 ↩