求助，fortran读数据出问题

mofangbao

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

                               
登录/注册后可看大图

晕，我是说你的文件是不是每个都只有in*nt那么大，如果不是的话你第二个文件开始record定位在第二个记录， ...

也许他不明白irec是干啥用的~不过现在应该明白了

lqouc

mofangbao 发表于 2013-9-3 15:52

                               
登录/注册后可看大图

也许他不明白irec是干啥用的~不过现在应该明白了

其实吧，我也不是特别明白，嘿嘿，我都是瞎猜的。作为一个菜鸟，压力山大啊。

寄居蟹

mofangbao 发表于 2013-9-3 15:49

                               
登录/注册后可看大图

不好意思，之前没仔细看，你这个是因为你的irec在读取下一个文件的时候没有初始化为1，比如第一个文件从开头 ...

谢谢，果然FORTRAN自学的一知半解

寄居蟹

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

                               
登录/注册后可看大图

晕，我是说你的文件是不是每个都只有in*nt那么大，如果不是的话你第二个文件开始record定位在第二个记录， ...

谢谢详细解释

冬日情愫+

mofangbao 发表于 2013-9-3 15:49

                               
登录/注册后可看大图

不好意思，之前没仔细看，你这个是因为你的irec在读取下一个文件的时候没有初始化为1，比如第一个文件从开头 ...

是的，没发现呀，进入下个文件前应该把记录号重设为1.   至于37*21后面要不要乘以4貌似和编译器有关吧，我以前用地好像都没有乘以4.

mofangbao

冬日情愫+ 发表于 2013-9-3 17:39

                               
登录/注册后可看大图

是的，没发现呀，进入下个文件前应该把记录号重设为1.   至于37*21后面要不要乘以4貌似和编译器有关吧，我 ...

曾经写过一个子函数来测试是否需要乘以4,这个其实每个编译器都可以自己设置的,原理就是二进制写入一个4个字节的文件,然后用*4的方式去读取,捕获错误,如果默认不需要乘以4,但是你乘以了,是会出错的,这样就知道不应该乘以4了.函数如下:

1. subroutine testbyte(byt)
   
2. integer val,byt
   
3. byt=4
   
4. open(1,file='test.dat',status='replace',form='binary')
   
5. write(1)123
   
6. close(1)
   
7. open(1,file='test.dat',status='old',access='direct',form='unformatted',recl=byt)
   
8. read(1,rec=1,err=100)val
   
9. goto 200
   
10. 100 continue
    
11. byt=1
    
12. 200 continue
    
13. close(1)
    
14. call system("del test.dat")
    
15. endsubroutine
    

subroutine testbyte(byt)<br /> integer val,byt<br /> byt=4<br /> open(1,file='test.dat',status='replace',form='binary')<br /> write(1)123<br /> close(1)<br /> open(1,file='test.dat',status='old',access='direct',form='unformatted',recl=byt)<br /> read(1,rec=1,err=100)val<br /> goto 200<br /> 100 continue<br /> byt=1<br /> 200 continue<br /> close(1)<br /> call system("del test.dat")<br /> endsubroutine<br />