gfortranで環境変数(FORT10等)を使って、ファイルとUNITを結合するためのパッチを作成しました(xTAPP.gfortran.diff)。
Debian (wheezy)でビルドする際の手順は以下の通りです。
なお、このI/O結合ライブラリの元ソースやテストコードは https://github.com/wistaria/giounit で公開しています。
早速ですが、先ほどのパッチは元素の種類が一種類の時しか上手く動きませんでした。giounit_open.Cに数行追加した修正版のパッチを添付します。Siの例とquartzの例について動作確認しました。元素が6種類までであれば動くと思います。
コードを作成してくださり、ありがとうございます。
取り込む方向で検討します。