求助，fortran读数据出问题

寄居蟹

program main


     PARAMETER (K_MAX=1, NT=37,IN=21,nV=2) !nv means the number of variables
! define 2D array:
       real,dimension(in,nt,nv) :: va
             character variable(2)*10
          data variable/'h','t'/
irec=1
do m=1,nv
! now read data
open(m,file=''//trim(variable(m))//'_0813.dat',form='binary',&
    recl=nt*in*4)
    read(m,rec=irec)((va(i,j,m),i=1,in),j=1,nt)
    irec=irec+1
        close(m)
end do
stop
end
出现问题如下：
我把两个文件分开一起读的时候， program main
     PARAMETER (K_MAX=1, NT=37,IN=21)
! define 2D array:
  real,dimension(in,nt) :: t,h
irec=1
! now read data
       open(21,file='T_0813.dat',form='unformatted',&
              recl=nt*in*4)
         read(21,rec=irec)((t(i,j),i=1,in),j=1,nt)
         irec=irec+1
        close(21)

       open(20,file='h_0813.dat',form='unformatted',&
              recl=nt*in*4)
         read(20,rec=irec)((h(i,j),i=1,in),j=1,nt)
         irec=irec+1
        close(20)
stop
end

也是读第二个文件就出现相同的错误。看到一个文章上说fortran错误代码36是一下解释
severe (36): Attempt to access non-existent record
FOR$IOS_ATTACCNON. A direct-access READ or FIND statement
attempted to access beyond the end of a relative file (or a sequential file
on disk with fixed-length records) or access a record that was previously
deleted from a relative file.
求各位指导下错误原因！谢谢！

维欻?

提示是没有找到第二个文件，你看看文件名是不是不对应

冬日情愫+

建议楼主不要使用10以内的数字当文件号吧，你可以用其他一个定值，反正你的循环结束前你都有关闭该文件。

寄居蟹

mofangbao 发表于 2013-9-3 14:25

                               
登录/注册后可看大图

close(m)放到循环外面去  赶紧多看基础知识吧  这种问题下次就自己解决

放到循环外面也还是相同的错误

寄居蟹

维欻? 发表于 2013-9-3 10:37

                               
登录/注册后可看大图

提示是没有找到第二个文件，你看看文件名是不是不对应

文件绝对存在，吧这两个文件顺序调换也还是读不了第二个文件

寄居蟹

冬日情愫+ 发表于 2013-9-3 11:12

                               
登录/注册后可看大图

建议楼主不要使用10以内的数字当文件号吧，你可以用其他一个定值，反正你的循环结束前你都有关闭该文件。

也不是这个原因呐

lqouc

你的每个文件是多大啊？

寄居蟹

lqouc 发表于 2013-9-3 15:14

                               
登录/注册后可看大图

你的每个文件是多大啊？

每个文件基本是4KB

lqouc

寄居蟹 发表于 2013-9-3 15:20

                               
登录/注册后可看大图

每个文件基本是4KB

晕，我是说你的文件是不是每个都只有in*nt那么大，如果不是的话你第二个文件开始record定位在第二个记录，但是你并没有第二个记录，所以你最开始的报错就会说你在尝试读取一个不存在的记录。
你读文件的目的是什么？把每个文件内的第一个记录都读出来么？不是很明白你为什么把文件和设备号还有记录数在一个循环里面转，这样很容易出错啊。

mofangbao

不好意思，之前没仔细看，你这个是因为你的irec在读取下一个文件的时候没有初始化为1，比如第一个文件从开头的地方读取了37*21*4个字节，你第二个文件如果不初始化irec=1的话，由于上面irec=irec+1，所以第二个文件是从37*21*4之后的字节开始读取的，当然超出了文件范围。
造成这个错误的原因可能是你只是看到了别人代码里面用到了irec=irec+1，却不知道这是干啥用的。。。（楼上的解释也是这个意思）