subroutine read_file( file_name, nx, ny, rec_num, var ) ! 出力結果読み取りルーチン ! 本ルーチンでは, ダイレクトアクセスを読み出す際, 1 変数のバイト数を 4 バイト ! と仮定して読み出すので, 4 バイト以外のファイルを読み出すときは, ! やはりソースファイルの書き換えが必要となる(要修正) implicit none integer, intent(in) :: nx ! データの x 方向の個数 integer, intent(in) :: ny ! データの y 方向の個数 integer, intent(in) :: rec_num ! 読み出すデータのレコード番号 character(*), intent(in) :: file_name ! 読み出すデータファイル名 real, intent(inout) :: var(nx,ny) ! 読み出すデータ integer :: i, j, er, k ! 作業用配列 open(unit=11, file=file_name, access='direct', recl=4*nx*ny, status='old') read(11,rec=rec_num) ((var(i,j),i=1,nx),j=1,ny) close(unit=11, status='keep') end subroutine read_file