Fortran Coder

标题: Fortran用stream流读取二进制文件问题 [打印本页]

作者: loveboyzmj 时间: 2015-7-30 18:23
标题: Fortran用stream流读取二进制文件问题
第一次在fcode发帖，初学fortran,如有不足请指教。
现有两个文件内容都是一样的，一个是有格式文件(*.FGRID)一个是无格式二进制文件(*.GRID)，其中有格式文件我用Editplus打开，内容如下：

'DIMENS  '          3 'INTE'
      10       10       25

无格式二进制文件,我想用Fortran的stream流打开，但是未能得到跟上边同样的结果，运行后我得到结果是：

DIMENS    50331648 INTE
268435456 201326592 167772160

主要是我要的整数 3 10 10 25 这几个数怎么一下变得那么大了呢？不明白，怎么解决呢？！
请大牛指教啊》》》》》》》》》》》》》》》》
百思不得其解，我的fortran代码如下：

[Fortran] 纯文本查看 复制代码

Program www_fcode_cn

  Use, Intrinsic :: ISO_C_BINDING      !// 使用内部模块，保持与 C 语言变量类型一致

  Implicit None

  character(len=4) :: flag0             !// 4字节

  character(len=8) :: flag1             !// 8字节

  character(len=4) :: dtype             !// 数据类型

  integer*4 :: i,j,k

  integer*4 :: ii,jj,kk

  Open( 12 , File = "INJECTWATER_FORWUNAN.GRID" , access="stream" , form = "unformatted" )

  Read( 12 ) flag0,flag1,i,dtype

  Write( * , * ) flag1,i,dtype

  Read( 12 ) ii,jj,kk

  Write( * , * ) ii,jj,kk

 !  Inquire( 12 , Pos = i )

 !  Write( * , * ) i

  Close( 12 )

 End Program www_fcode_cn

以上为我的代码，是哪个地方出问题了呢？！望求解，谢谢。

作者: fcode 时间: 2015-7-30 19:19
把文件传上来

作者: loveboyzmj 时间: 2015-7-30 20:01

fcode 发表于 2015-7-30 19:19
把文件传上来

共两个文件，见压缩包。
[attach]798[/attach]

作者: loveboyzmj 时间: 2015-7-30 20:04

fcode 发表于 2015-7-30 19:19
把文件传上来

Grid.rar (164.33 KB, 下载次数: 21) ，见压缩包，谢谢。

作者: fcode 时间: 2015-7-30 22:19
您的问题有两个：
1.您的数据是一个无格式顺序文件。所以如果你要用 stream 读取，需要在每个 read 前后各加一个 n。
之前你的代码在前面加了 flag0，但没有在后面加
2.您的数据是一个大端字节序（big-endian），而你的编译器是默认小端字节序的（little-endian），您需要在 Open 的时候指定为大端字节序。
当然，指定的方法不标准，因为语法没有规定大小端字节序。如果以下的代码在你的编译器上不支持，请告诉我你使用什么编译器，我再为你查找方法。

[Fortran] 纯文本查看 复制代码

Program www_fcode_cn

  Implicit None

  character(len=8) :: flag1             !// 8字节

  character(len=4) :: dtype             !// 数据类型

  integer(kind=4) :: i,j,k,n

  integer(kind=4) :: ii,jj,kk

  Open( 12 , File = "INJECTWATER_FORWUNAN.GRID" ,convert='BIG_ENDIAN', access="stream" , form = "unformatted" )

  Read( 12 ) n,flag1,i,dtype,n

  Write( * , * ) flag1,i,dtype

  Read( 12 ) n,ii,jj,kk,n

  Write( * , * ) ii,jj,kk

  Close( 12 )

 End Program www_fcode_cn

作者: loveboyzmj 时间: 2015-7-30 23:31
非常感谢fcode的耐心解答，明天我试一下。

作者: loveboyzmj 时间: 2015-7-31 09:28

fcode 发表于 2015-7-30 22:19
您的问题有两个：
1.您的数据是一个无格式顺序文件。所以如果你要用 stream 读取，需要在每个 read 前后各 ...

fcode你好

,感谢你的解答，用了你给我的代码，问题已解决，以后有不懂还得向你请教，太感谢了！
我现在用simply fortran，非常感谢论坛。

作者: loveboyzmj 时间: 2018-11-13 11:21
回顾当年的问答，现在依然获益匪浅！！！

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)