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