導入：データ不整合によるプログラム停止を防ぐ
数値計算エンジニアにとって、外部ファイルからのデータ読み込みは避けて通れない作業です。しかし、入力データに予期せぬ文字が含まれていたり、フォーマットが崩れていたりすると、プログラムは即座に強制終了してしまいます。大規模なシミュレーションにおいて、数時間かけた計算の末にデータ不整合で停止するのは避けたいものです。今回紹介する「READ文のERR=分岐」は、異常なデータ行を検知し、プログラムを止めずに「スキップ」や「ログ記録」を行うための非常に重要なテクニックです。
基礎知識：ERR=分岐とは何か
FortranにおけるREAD文のERR引数は、入出力操作中にエラーが発生した場合に、指定されたラベル番号へ処理をジャンプさせる機能です。
通常、数値データ（F10.0など）を読み込む場所に文字列が混入すると、Fortranのランタイムライブラリは致命的なエラーを出し停止します。ERR=888 を指定することで、エラー発生時に888番の行へ制御が移り、そこで適切なエラーハンドリング（変数のリセットやエラーメッセージの出力）を行うことが可能になります。
実装・解決策：エラーを「無視」ではなく「管理」する
単にエラーを無視するのではなく、エラーが発生した行番号や内容をログに残すことで、データクレンジングの効率が劇的に向上します。ループ内でREAD文を使用し、ERRラベルでエラー処理ルーチンへ飛ばすことで、読み込み不能な行をスキップして次の行の読み込みへ進む「耐障害性の高いインポーター」を構築できます。
サンプルプログラム
以下のコードは、ファイル内の数値データを読み込み、エラー発生時にその旨を通知して処理を継続する例です。
program read_robust
implicit none
integer :: i, ios, err_count
real :: x
err_count = 0
open(10, file='data.txt', status='old')
do i = 1, 10
! ERR=888 を指定し、読み込み失敗時に888行目へジャンプさせる
read(10, '(F10.0)', err=888) x
print , "読み込み成功:", x
cycle ! 成功時は次のループへ
888 continue ! エラー発生時のジャンプ先
err_count = err_count + 1
print , "エラー発生: ", i, "行目のデータが不正です。"
end do
close(10)
print , "総エラー数:", err_count
end program read_robust
応用・注意点：現場で役立つアドバイス
1. IOSTATとの併用：
ERR=分岐はエラーの発生有無を判定するのには便利ですが、何が起きたかを知るためには IOSTAT= 変数 を併用するのがベストプラクティスです。IOSTATはエラーコードを返すため、原因の切り分けが容易になります。
2. 無限ループの回避：
ERR=分岐で飛ばした際、ファイルポインタが正しく進まない（異常な行に留まる）と、同じエラーで無限ループに陥る可能性があります。エラー発生時は、必要に応じてスキップ処理（行全体を読み飛ばすなど）を併用してください。
3. デバッグの重要性：
実運用ではエラーを無視するだけでなく、どのファイル、どの行でエラーが起きたかを必ずログに出力してください。データ作成者のミスを迅速に指摘できる環境を作ることが、エンジニアとしての信頼性に繋がります。
コメント