fortran 77 代码write问题

fortran中的write用法大致如下：write(设备号，格式)。Fortran对这个设备号的要求是要求在设备号上的数据是整型。11是整型，自己定义的整型变量也是整型。

举个栗子：

integer::n

n=11

write(n,*)"hello world!"

的意思就是说再文件号为11的文件中写入“hello world！”这句话。

以一个 m x n 矩阵 A 为例，介绍Fortran二进制文件的读写

1、打开与关闭文件：open/close

open(unit=99,file='data.dat',access='direct',form='unformatted',recl=8)

close(99)

其中：

unit=99 ! 设备号，Fortran 77规定取值范围为1～99

file='data.dat' ! 文件名

access='direct'  ! 以直接访问的方式 *** 作文件

form='unformatted' ! 以无格式方式进行 *** 作，当access='direct'时，form缺省为'unformated'

recl=8  ! 每个记录的大小

上面比较重要的是 recl，我们读写数据时，每次读写的是一个记录，长度就由recl指定，

一般表示字节数，如recl=8表示记录长度为8个字节。

但在某些旧的编译器下以4个字节为单位，如recl=2表示记录长度为8字节，

这个问题需要参考编译器手册，我们这里是以字节为单位。

2、写数据：write

write(99,rec=k) A(i,j)  ! 将 A(i,j) 作为第 k 条记录写入到文件中

这里的 rec 表示记录的序号，每次只能写入一个记录，但长度可以自己定。

将一个双精度矩阵写入文件，可以有下面三种方式实现：

(1) 将每个元素看成一个记录，recl=8

      open(unit=99,file='data.dat',access='direct',form='unformatted',

     &     recl=8)

      k = 0

      do j=1,n

         do i= 1,m

            k=k+1

            write(99,rec=k) A(i,j) 

         enddo

      enddo

      close(unit=99)

(2) 把一列看成一个记录，长度为 8*m

 

 open(unit=99,file='data.dat',access='direct',form='unformatted',

     &     recl=8*m)

      do j=1,n

         write(99,rec=j) A(1,j) 

      enddo

      close(unit=99)

(3) 把整个看成一个记录，长度为 8*m*n

      open(unit=99,file='data.dat',access='direct',form='unformatted',

     &     recl=8*m*n)

      write(99,rec=1) A 

      close(unit=99)

3、读数据：read

read(99,rec=k) B(i,j) ! 将第k个记录读入到 B(i,j) 中

我们可以通过下面三种方式读一个矩阵

(1) 将每个元素看成一个记录，recl=8

   open(unit=99,file='data1.dat',access='direct',form='unformatted',

     &     recl=8)

      k = 0

      do j=1,n

         do i= 1,m

            k=k+1

            read(99,rec=k,iostat=ferr) A(i,j) 

         enddo

      enddo

      close(unit=99)

(2) 把一列看成一个记录，长度为 8*m

 

 open(unit=99,file='data1.dat',access='direct',form='unformatted',

     &     recl=8*m)

      do j=1,n

         read(99,rec=j,iostat=ferr) (A(i,j), i=1,m)

      enddo

      close(unit=99)

(3) 把整个看成一个记录，长度为 8*m*n

     open(unit=99,file='data1.dat',access='direct',form='unformatted',

     &     recl=8*m*n)

      read(99,rec=1,iostat=ferr) ((A(i,j), i=1,m), j=1,n)

      close(unit=99)

4、读取整个文件

     open(unit=99,file='data.dat',access='direct',form='unformatted',

     &     recl=8)

      k = 0

      do while (.true.)

         k = k + 1

         read(99,rec=k,iostat=ferr) x(k)

         if (ferr/=0) exit

      enddo

1x表示输出一个空格，以此类推。

估计你这个是固定格式的源程序，那个“+”号是续行符号。

如果是这样输出结果就应该是：

------------------------------------------------------------------

t(s) Pa(kPa)Pw(kPa)Set(cm)

------------------------------------------------------------------

---