gfortran与Fortran 6.x 中二进制问题(form='binary'/'unformatted')的解决

/xin村儿/

gfortran与Fortran 6.x 中二进制问题(form='binary'/'unformatted')的解决

纠结了多年的’心病‘，今天偶然蹦出’经验‘了——

因为gfortran不支持binary格式，二进制只能以unformatted表达，有帖子说使用gfortran编译时直接把binary换成unformatted即可，但实践中还是出现错误，为什么呢？请往下看：

在Fortran6.6中，binary和unformatted均可使用，但用的时候有区别：      
OPEN(10,FILE='路径/AAAAA.grd',FORM="unformatted",status='old',ACCESS='DIRECT', RECL=NX*NY)   !!!注意不*4
等价于
OPEN(10,FILE='路径/AAAAA.grd',FORM="binary",        status='old',ACCESS='DIRECT', RECL=NX*NY*4)   !!!注意要*4


而在gfortran中，二进制只能使用unformatted，但用的时候和Fortran6.6又有区别：      
OPEN(10,FILE='路径/AAAAA.grd',FORM="unformatted",status='old',ACCESS='DIRECT', RECL=NX*NY*4)   !!!注意要*4
等价于Fortran6.6中的
OPEN(10,FILE='路径/AAAAA.grd',FORM="unformatted",status='old',ACCESS='DIRECT', RECL=NX*NY)   !!!注意不*4

读写时均为： READ/WRITE(10, rec=IREC) ARR      !!!假设数组为ARR(I,J), 注意有格式时不能直接使用数组名ARR
个人经验，有待更多验证！
补充：http://bbs.06climate.com/forum.php?mod=viewthread&tid=27204











来自群组: 成于大气 信达天下

zhouling

非常感谢楼主的经验分享，非常有帮助~

兰溪之水

这个乘4或不乘4不能随便这么下结论的。。因为这些默认字节长度编译器可以设置的

/xin村儿/

兰溪之水 发表于 2014-11-3 15:39
这个乘4或不乘4不能随便这么下结论的。。因为这些默认字节长度编译器可以设置的

谢谢兰溪哥指点，这些问题没有深究过，或许正因如此，所以不懂老出错走弯路，只是自己总结经验在这里

自导自演

兰溪之水 发表于 2014-11-3 15:39
这个乘4或不乘4不能随便这么下结论的。。因为这些默认字节长度编译器可以设置的

兰溪老师，可不可以详细点讲一讲*4或不*4的意义

看云淡风清

兰溪之水 发表于 2014-11-3 15:39
这个乘4或不乘4不能随便这么下结论的。。因为这些默认字节长度编译器可以设置的

@兰溪之水 大大，我是在gfortran下编译的，我的问题是读取出错但是不知道怎么改可以帮忙看看嘛，我的程序如下：
program main
integer,parameter::nxy=144*73,nt=54,ct=0
real a(nxy,nt),ind(nt),r(nxy),bb(nxy),cc,t,n(nxy),aa(nxy),dd(nxy),ee(nxy),ff,gg(nxy)
real x1(nt),y1(nxy,nt)
open(3,file='hgt_winter_1961_2014.grd',form='unformatted',status='old',access='direct',recl=nxy*nt*4)
open(4,file='tp.grd',form='unformatted',status='old',access='direct',recl=nxy*nt*4)
open(7,file='rt.grd',form='unformatted')
open(8,file='tt.grd',form='unformatted')
read(3)((a(i,it),i=1,nxy),it=1,nt)
read(4)(ind(it),it=1,nt)

老泥巴

自己摸索了一下，写了一个小例子：
program ex20220421
    implicit none
    real::a=10.23,b
    !gfortran不支持binary，这边写入.grd需要指定access='stream'
    open(unit=1,file='svp.grd',access='stream')
    !切记写入的时候不可指定格式，也不可以写write(1,*)
    write(1) a
    close(1)
   
    !读取文件输出看看
    open(unit=2,file='svp.grd',access='stream')
    !注意二进制文件无格式，所以读取的时候直接read(2)，不可以写read(2,*)
    read(2) b
    print*,b
    close(2)
end program ex20220421