複数のファイルに出力

ここまでに示した例では 1 つのファイルに全ての変数と属性の出力を行いま した. しかしながら複数のファイルへ出力することもできます.

ここでは 2 次元の熱対流を計算するプログラムを例に解説します. 出力する 変数は 2 つ(流線関数と温度)です. 1 つのファイルに出力を行うプログラム が benard_1.f90, 2 つのファイルに分けて出力を行うプログラムが benard_2.f90 です.

benard_2.f90 では 34 行目で

type(GT_HISTORY) :: hst_psi, hst_temp

とし, それぞれ流線関数と温度を出力するファイルを示す GT_HISTORY 型構造体変数 hst_psi と hst_temp を宣言します.

次に 77 行目と 95 行目の HistoryCreate で出力ファイル psi.nc と temp.nc を作成します. このとき最後の引数には 34 行目で宣言した hst_psi と hst_temp を代入します.

call HistoryCreate( &                             ! ヒストリー作成
  & file='psi.nc', title='convection in rotating annulus', &
                    : 
  & history=hst_psi ) 

                    :

call HistoryCreate( &                             ! ヒストリー作成
  & file='temp.nc', title='convection in rotating annulus', &
                    : 
  & history=hst_temp )

117 行目から 118 行目までが出力です. HistoryPut の第 3 引数には HistoryCreate で指定したそれぞれの出力先ファイルに対応する GT_HISTORY 型構造体変数を代入します.

call HistoryPut('psi',  trans_g(psi),  hst_psi)
call HistoryPut('temp', trans_g(temp), hst_temp)

では, 実際に nc ファイルを作成してみましょう. benard_1.f90 と benard_2.f90 は c2pack モジュールを利用するので, まず c2pack.f90 を コンパイルします.

c2pack.f90 と ispack_snip.f のコンパイル

c2pack.f90 と ispack_snip.f を ダウンロードし, 以下のようにコンパイルしてください.

$ gt5frt -c c2pack.f90 ispack_snip.f

ispack_snip.f は ISPACK の一部のパッケージを抜粋したもので, c2pack モジュールに必要な サブルーチンが含まれています.

コンパイルの結果, c2pack.o, c2pack.mod, ispack_snip.o が作成されれば OK です.

benard_1.f90 のコンパイル

では, benard_1.f90 をコンパイルしてみましょう.

$ gt5frt benard_1.f90 c2pack.o ispack_snip.o

すると, 実行ファイルを a.out が作成されます. そこで,

$ ./a.out

と実行させると, benard_1.nc が作成されます.

benard_2.f90 のコンパイル

では, 今度は benard_2.f90 をコンパイルしてみましょう. コンパイルの仕方は benard_1.f90 と同様です.

$ gt5frt benard_2.f90 c2pack.o ispack_snip.o

できた実行ファイル a.out を実行してみましょう.

$ ./a.out

すると今度は, 流線関数と温度のデータがそれぞれ, psi.nc と temp.nc に出力されます.