NC文件转dat后用fortran读取出问题了

一颗梅子在南京

我要从nc文件中提取一下某一经纬度、某一时刻，17层的u，v风速，我的gs文件如下：
'reinit'
'sdfopen e:\nc\uwnd.1971.nc'
'sdfopen e:\nc\vwnd.1971.nc'
'set gxout fwrite'
'set fwrite e:\nc\uvwnd.1971.dat'
'set lon 0'
'set lat 52.5'
'set t 1'
z=1
while(z<=17)
'set z 'z
'd uwnd.1'
'd vwnd.2'
z=z+1
endwhile
'disable fwrite'

然后我想用fortran转化一下看看数据，运行的时候出错了：
forrtl: severe (257): formatted I/O to unit open for unformatted transfers, unit
1, file e:\nc\uvwnd.1971.dat

我的fortran程序是这样的：
program main
implicit none
real wind(17,2)
integer i,j
open(1,file='e:\nc\uvwnd.1971.dat',form='binary')
open(2,file='e:\nc\1971.txt')
do i=1,17
   read(1,*)(wind(i,j),j=1,2)
   write(2,*)(wind(i,j),j=1,2)
enddo
end

这是哪错了呢？

做个霸气的木头

binary数据的读取是只加文件号的，没有格式，read错了，再者你是u、v都写到文件里了，你只读一个吗

言深深

看看这个gs文件可for文件吧http://bbs.06climate.com/forum.php?mod=viewthread&tid=10848
无格式数据读取方式的时候需要注意隐循环

一颗梅子在南京

做个霸气的木头 发表于 2012-11-18 15:20

                               
登录/注册后可看大图

binary数据的读取是只加文件号的，没有格式，read错了，再者你是u、v都写到文件里了，你只读一个吗

我这个wind是二维的呀，就包括u,v了呀

一颗梅子在南京

一颗梅子在南京 发表于 2012-11-18 18:37

                               
登录/注册后可看大图

我这个wind是二维的呀，就包括u,v了呀

可以了，谢谢斑竹

一颗梅子在南京

言深深 发表于 2012-11-18 17:28

                               
登录/注册后可看大图

看看这个gs文件可for文件吧http://bbs.06climate.com/forum.php?mod=viewthread&tid=10848
无格式数据读取 ...

多谢言斑竹

言深深

一颗梅子在南京 发表于 2012-11-18 18:47

                               
登录/注册后可看大图

多谢言斑竹

出来就好哈

做个霸气的木头

一颗梅子在南京 发表于 2012-11-18 18:37

                               
登录/注册后可看大图

我这个wind是二维的呀，就包括u,v了呀

我表示看错了

blui

为什么我提取的每一层数据都不符合的。。。

一颗梅子在南京

blui 发表于 2013-4-8 11:47

                               
登录/注册后可看大图

为什么我提取的每一层数据都不符合的。。。

不知道你怎么弄的，多去论坛看看相关的帖子吧！