NCL读取二进制文件时的问题

花楹

读取二进制文件时，只出现第一行的数据，其他行全部为999，不知道哪出错了，下面是我的ncl脚本：
请各位帮忙看看，二进制文件 1961 (211.5 KB, 下载次数: 16)

2012-4-17 19:06 上传

点击文件名下载附件

hzwjy

顺序存储和直接存储的数据分开讲，
1.
FORTRAN程序写直接存储无格式数据的代码是：
OPEN(21,file=filename, form='unformatted', access='direct', action='read', recl=....)
对应的NCL读取函数为fbindirread

现在假定有组28层1度*1度的数据T，大小为(360*181*28)
FORTRAN写的时候可以这样
OPEN(21,file="file1.dat", form='unformatted', access='direct', recl=360*181*28*4)
WRITE(21,rec=1) T
一次把数据都写入文件

也可以
OPEN(21,file="file2.dat", form='unformatted', access='direct', recl=360*181*4)
DO i = 1, 28
  WRITE(21,rec=i) T(:,:,i)
ENDDO
写个循环，一次写一层数据

这两种方法出来的文件结果会是完全一样的（顺序存储就不是了）！
相应的，无论你用什么方式写这个数据，NCL读取时也可以有多种读取方法：
  T = fbindirread(filename, -1, (/28,181,360/), "float")
中间的-1表示一次将整个文件读入

也可以写个循环，一次只读取一层

2.
FORTRAN顺序存储的无格式数据的代码是：
OPEN(21,file=filename, form='unformatted', access='sequential')
还是那个数据T

OPEN(21,file="file3.dat", form='unformatted', access='sequential')
WRITE(21) T
一次写完所有数据

OPEN(21,file="file4.dat", form='unformatted', access='sequential')
DO i = 1, 28
  WRITE(21) T(:,:,i)
ENDDO
一次写入一层数据
这两种写法得到的文件是不同的，文件大小就是不一样的！
因为在写顺序存储的二进制文件时，编译器会在每个数据记录添加额外的记录，无格式顺序文件中的第一个和最后一个字节是保留字节，第一个字节的值为75，最后一个字节的值为130。Fortran使用这些字节作为错误检测和文件结束的判断。
（详细参考FORTRAN教材）

对应的读取数据的NCL函数为fbinrecread，由于记录块大小有特殊要求，编译FORTRAN代码时要加参数：
gfortran foo.f90 -frecord-marker=4 -o foo.exe
运行foo.exe生成的数据就是NCL可以读取的。
一次写入28层的顺序二进制文件的读取：
  T = fbinrecread("file3.dat", 0, (/28,181,360/), "float")

一次写入1层的顺序二进制文件的读取：
  T = new((/28, 181, 360/), "float")
  do i = 0, 27
    T(i,:,:) = fbinrecread("file4.dat", i, (/181,360/), "float")
  end do

不熟悉的话，这个函数很容易用错。而且一般情况下FORTRAN程序的编译都不会用 -frecord-marker=4（以后NCL可能会添加额外的功能，就不用设置这个编译参数了）。
所以，FORTRAN程序最好用直接存储的二进制格式来写。或者用netCDF格式，NCL操作netCDF格式的数据是很方便的。
===============

这样说明够详细了吧。

花楹

hzwjy 发表于 2012-4-19 19:28

                               
登录/注册后可看大图

顺序存储和直接存储的数据分开讲，
1.
FORTRAN程序写直接存储无格式数据的代码是：

太强大了，非常感谢啊~我试试哈，就是还有一个疑问，你说 这个编译我在fortran里么有见过哎，在哪个地方嘞？万分感谢啊~

mofangbao

你直接把脚本贴出来吧 截图别人还得一点点的敲  
虽然我不会ncl，但是提问说 不知道哪儿错了是最容易不想帮你看的！  自己最起码先一点点的调试出错误在哪一句话产生的，自己对这个问题的已有的了解，这都是为了让别人更好的帮你，请楼主设身处地的想一想

mofangbao

你直接把脚本贴出来吧 截图别人还得一点点的敲  
虽然我不会ncl，但是提问说 不知道哪儿错了是最容易不想帮你看的！  自己最起码先一点点的调试出错误在哪一句话产生的，自己对这个问题的已有的了解，这都是为了让别人更好的帮你，请楼主设身处地的想一想

花楹

mofangbao 发表于 2012-4-17 19:29

                               
登录/注册后可看大图

你直接把脚本贴出来吧 截图别人还得一点点的敲  
虽然我不会ncl，但是提问说 不知道哪儿错了是最容易不想帮 ...

我调试了一下午了，弄不出来才发的帖子，我也是刚学啊~不会一弄不出来就发的。。。。也没找到论坛上相近的东西

花楹

grd.ncl (107 Bytes, 下载次数: 108)

2012-4-17 20:36 上传

点击文件名下载附件

mofangbao

花楹 发表于 2012-4-17 20:34

                               
登录/注册后可看大图

我调试了一下午了，弄不出来才发的帖子，我也是刚学啊~不会一弄不出来就发的。。。。也没找到论坛上相近的 ...

我的意思是你要把你已经做的描述一下啊，这样让别人有个切入点

花楹

mofangbao 发表于 2012-4-17 20:37

                               
登录/注册后可看大图

我的意思是你要把你已经做的描述一下啊，这样让别人有个切入点

就五行命令啊，第三行就是读取的，应该是（94,192）这有问题，但是怎么改都改不对

hzwjy

这个数据是顺序存储还是直接存储？
不同的存储方式所用的NCL读取函数是不同的

fbinrecread是用来读取顺序存储的无格式FORTRAN数据的
我个人不太喜欢用这个函数读取，因为这个函数读取数据的记录标记（record markers）为
4-bytes，而linux下常用的gfortran等编译器默认为8-byte，如果真有需要用这个NCL函数读取二进制数据，生成数据用的FORTRAN代码需要用"-frecord-marker=4"参数进行编译。
这个函数的说明里提到了这些。
=======

我建议用直接存储格式的二进制数据，然后用fbindirread读取

花楹

hzwjy 发表于 2012-4-19 10:29

                               
登录/注册后可看大图

这个数据是顺序存储还是直接存储？
不同的存储方式所用的NCL读取函数是不同的

你好 我就是顺序读取的无格式二进制数据呢，不知道有什么问题呢？不知道您能帮忙看看嘛 我在附件中都给出来了哦 多谢啊

花楹

hzwjy 发表于 2012-4-19 10:29

                               
登录/注册后可看大图

这个数据是顺序存储还是直接存储？
不同的存储方式所用的NCL读取函数是不同的

您好，我用您说的fbindirread函数做了，是对的，谢谢哈~但是我按照您说的“生成数据用的FORTRAN代码需要用"-frecord-marker=4"参数进行编译。”然后用fbinrecread读取就不对，还是出错呢？你的意思是说fortran程序要这样写吗？open(2,file='d:\essay2\grid_ave\1961'.grd',form='unformatted',recl=4)