1. 導入：なぜデータが正しく読み込めないのか
数値計算の現場では、Fortranで計算した結果をC言語やPythonで読み込んだり、その逆を行ったりすることが頻繁にあります。しかし、バイナリデータをファイル経由でやり取りする際、「数値が全く別の値になってしまう」「読み込みの途中でエラーが発生する」というトラブルを経験したことはないでしょうか。その原因の多くは、データの書き込み形式における「アライメント（配置）」の不一致にあります。この記事では、この厄介な罠を回避し、異なる言語間で安全にデータをやり取りする方法を解説します。
2. 基礎知識：レコードマーカーとは何か
Fortranの従来の書式なし出力（`WRITE(UNIT) A`）では、データの前後を囲むように「レコードマーカー」と呼ばれる4バイトの制御情報が付与されます。これはデータブロックのサイズを示すもので、Fortran内部で読み書きする分には便利な仕組みですが、C言語の `fread` やPythonの `numpy.fromfile` などで直接ファイルを読み込むと、このマーカーまでデータとして扱われてしまい、計算結果が壊れる原因となります。現代の数値計算では、このマーカーを排除した「ストリームアクセス」を用いるのが鉄則です。
3. 実装/解決策：ACCESS=’STREAM’の活用
Fortranにおいて、レコードマーカーを出力させないためには、`OPEN` 文のオプションで `ACCESS=’STREAM’` を指定します。これにより、C言語の `fwrite` と全く同じように、メモリ上のバイナリデータをそのままファイルに書き出すことができます。これにより、データの先頭番地がズレることなく、他言語からの直接読み込みが可能になります。
4. サンプルプログラム：ストリームアクセスによるバイナリ出力
以下のコードは、`ACCESS=’STREAM’` を使用して安全にバイナリデータを書き出す例です。
PROGRAM write_binary
IMPLICIT NONE
INTEGER, PARAMETER :: dp = KIND(0.0D0)
REAL(KIND=dp) :: data(3) = [1.1D0, 2.2D0, 3.3D0]
INTEGER :: i
! ACCESS=’STREAM’を指定することで、レコードマーカーを付与せずに書き込む
OPEN(UNIT=10, FILE=’output.bin’, ACCESS=’STREAM’, STATUS=’REPLACE’)
! 配列をそのままバイナリとして出力
WRITE(10) data
CLOSE(10)
PRINT , “バイナリファイルの出力が完了しました。”
END PROGRAM write_binary
5. 応用・注意点：現場で役立つアドバイス
注意点1：エンディアンの確認
アライメントの問題を解決しても、計算機環境（アーキテクチャ）が異なると「エンディアン（バイト順）」の違いで数値が化けることがあります。異なる環境間でやり取りする場合は、ビッグエンディアンかリトルエンディアンかを確認し、必要に応じて変換処理を挟む必要があります。
注意点2：構造体のパディング
C言語側の構造体（struct）とFortranの配列を対応させる場合、コンパイラの「パディング（詰め物）」に注意してください。構造体のメンバ間に隙間が挿入される設定になっていると、単純なバイナリ読み込みでは位置がズレます。基本的には「パッキング」を有効にするか、単純な配列構造でデータをやり取りするのが最も安全です。
これらのポイントを押さえることで、言語間のデータ受け渡しにおけるデバッグ時間を大幅に削減できるはずです。
コメント