Fortran读写二进制文件问题

这个也有-

我用grads导出了两个计算出来的量，是二进制的，然后我打算用Fortran处理，我现在就直接读这个文件然后在把它输出成grads能够识别的二进制格式，但是发现输出后的文件和读入的文件的大小不一样了。
program main

integer i,j,t
integer,parameter::nx=101,ny=76,tt=125
real ri(nx,ny,tt), rhcha(nx,ny,tt)

open(10, file='E:\data\fnl201505\rirhcha.grd',form='binary',recl=nx*ny,access='direct')

irec=1
do t=1,tt
  read(10,rec=irec) ((ri(i,j,t), i=1,nx),j=1,ny)
  read(10,rec=irec) ((rhcha(i,j,t), i=1,nx), j=1,ny)
  irec=irec+1
end do



open(20, file='E:\data\fnl201505\rirhcha2.grd', form='binary', recl=nx*ny, access='direct')

irec=1
do t=1,tt
  write(20,rec=irec) ((ri(i,j,t), i=1,nx),j=1,ny)
  write(20,rec=irec) ((rhcha(i,j,t), i=1,nx), j=1,ny)
  irec=irec+1
end do

close(10)
close(20)

end
这是我的Fortran程序，中间没有进行计算，就是读和写，但是文件大小近差了好多。。。
使用grads画图也出现了错误。。。

王磊

http://v.fcode.cn/video-file_io_ascii.html 可以参考这期视屏，了解2进制文件的存储。

这个也有-

王磊 发表于 2015-11-30 23:34
http://v.fcode.cn/video-file_io_ascii.html 可以参考这期视屏，了解2进制文件的存储。

很好的视频，正在学习中。

这个也有-

问题已解决，在读和写的循环中再加入一句：irec=irec+1

do t=1,tt
  read(10,rec=irec) ((ri(i,j,t), i=1,nx),j=1,ny)
  irec=irec+1
  read(10,rec=irec) ((rhcha(i,j,t), i=1,nx), j=1,ny)
  irec=irec+1
end do
改成如上，就对了。