- 积分
- 3087
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-16
- 最后登录
- 1970-1-1
|
发表于 2012-11-16 21:47:35
|
显示全部楼层
版主,你好,我要从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
这是哪错了呢? |
|