FORTRANによるモデル出力結果の読み込み

• 以下のようなコードでモデルの出力をFORTRANプログラムに読み込むことができる。

「daily meanで一年間のデータの場合」

integer idim, jdim, kdim　　　　　　　　　　　!次元の数( 経度、緯度、高度 )
parameter( idim=128, jdim=64, kdim=32 )
real*4 gdat( idim,jdim,kdim )　　　　　　　　!読み込む変数
character head( 64 )*16　　　　　　　　　　!データのヘッダ


open( 10, file="ファイル名", form='unformatted' )　!読み込みたいファイルを開く
do n = 1, 365
　read( 10 ) head　　　　　　　　　　　　　　!ヘッダの読み込み
　read( 10 ) gdat　　　　　　　　　　　　　　!データの読み込み
enddo

　* モデルデータは1日（1時間ステップ）ごとにヘッダ→データという構造なので、ヘッダとデータを交互に読む。

読み込んだその後は…

任意の一点のデータを抽出

全球総量の計算