fortran读取二进制文件出错，求助！！！

weinihou

繁星一点 发表于 2017-7-27 20:28
我直接把recl那一句删掉了，出来的数据还是一样。
我就是要把二进制的数据提取出来，写成十进制的整数数 ...

要么试试 integer(kind = 4)  或者 integer*8  ？ 2这种短的我猜不用试

上面那个   ‘d hgt’  是提取整数的意思吧？我完全不懂那些

要么先计算下文件大小，和数据量匹配一下，看看到底提取出来的数据是哪种大小的

繁星一点

weinihou 发表于 2017-7-27 21:38
要么试试 integer(kind = 4)  或者 integer*8  ？ 2这种短的我猜不用试

上面那个   ‘d hgt’  是提取 ...

d hgt 就是提取hgt这个变量
我计算了提取出来的变量的个数是对的，只是值不对

繁星一点

river 发表于 2017-7-27 21:37
顺序没错啊，既然是用grads自带的函数写出来的文件没救一定是符合grads要求的格式的。你下载兰溪的grads2 ...

下载试了，果然成功了
瞬间有一种海阔天空的感觉
太谢谢了

weinihou

所以也不知道到底是哪里不合GrADS规还是有啥问题了
曲线救国绕开了Fortran二进制读取

river

繁星一点 发表于 2017-7-27 22:07
下载试了，果然成功了
瞬间有一种海阔天空的感觉
太谢谢了

兰溪那个脚本其实也是用了grads自带的命令，所以楼主还是要加强对grads的了解，grads还是很强大的

繁星一点

river 发表于 2017-7-28 11:58
兰溪那个脚本其实也是用了grads自带的命令，所以楼主还是要加强对grads的了解，grads还是很强大的

前两天看到过那个帖子，当时一根筋的想用自己的这张方法，所以钻牛角尖了

繁星一点

weinihou 发表于 2017-7-27 22:56
所以也不知道到底是哪里不合GrADS规还是有啥问题了
曲线救国绕开了Fortran二进制读取

是的，我怀疑是fortran顺序读取的问题。
画图做了检验，grd文件提取肯定没有出错。那唯一出错的就是读数据了，但是又查不出哪里错了

yfw07

很厉害，谢谢楼主

霹雳电光鼠

读取的地方这么改一下试试  
open(10,file='e:\nc\199806.grd',form='binary',access="direct",recl=nx*ny)  
         do t=1,nt
                read(10) ((hgt(x,y,nt),x=1.nx),y=1,ny)
         enddo

lgz

recl=nx*ny*nt