fortran数据读入问题 求指点

心语

这是楼主的.gs文件
'reinit'
'sdfopen e:/lw/netcdf2.nc'
'set fwrite e:/lw/quv.grd'
'set gxout fwrite'
'set t 1'
'd u'
‘d v'

'disable fwrite'
;

下图是netcdf2.nc的资料

mofangbao

把fortran代码贴一下

lqouc

你想输出什么？这个gs写的很有问题。
fortran的代码在哪里？报错是什么？
建议楼主再看看‘提问的智慧’。

心语

mofangbao 发表于 2014-4-20 10:13
把fortran代码贴一下

program main
implicit none
integer:: m,n,z
parameter (m=240,n=121,z=27)
integer i,j,k
real u(m,n,z),v(m,n,z),q(m,n,z)
open (1,file='e:/lw/quv.grd',form='binary')
do k=1,z
read(1,*)((u(i,j,k),i=1,m),j=1,n)
enddo
do k=1,z
read(1,*)((v(i,j,k),i=1,m),j=1,n)
enddo
end

心语

                               
登录/注册后可看大图

心语

lqouc 发表于 2014-4-20 10:13
你想输出什么？这个gs写的很有问题。
fortran的代码在哪里？报错是什么？
建议楼主再看看‘提问的智慧’ ...

楼主小白，fortran代码在上一楼哈，报错：formatted I/O to unit open for unformatted transfers,unit 1,file e:/...，

lqouc

心语 发表于 2014-4-20 10:23
楼主小白，fortran代码在上一楼哈，报错：formatted I/O to unit open for unformatted transfe ...

read里面的缺省格式设置（那个*）去掉（binary不需要format），gs里面维度最好set完整，zdef需要用循环输出。

心语

lqouc 发表于 2014-4-20 10:31
read里面的缺省格式设置（那个*）去掉（binary不需要format），gs里面维度最好set完整，zdef需要用循环输 ...

gs文件改了一下，请问这样可以吗？
'sdfopen e:/lw/netcdf2.nc'
'set fwrite e:/lw/quv.grd'
'set gxout fwrite'
'set t 1'
m=1
while(m<=27)
'set z 'm''
'd u'
'd v'
m=m+1
endwhile
'disable fwrite'
;

fortran代码
open (1,file='e:/lw/quv.grd')
read(1)(((u(i,j,k),i=1,m),j=1,n),k=1,z)
read(1)(((v(i,j,k),i=1,m),j=1,n),k=1,z)  运行了还是不能读入

lqouc

心语 发表于 2014-4-20 10:53
gs文件改了一下，请问这样可以吗？
'sdfopen e:/lw/netcdf2.nc'
'set fwrite e:/lw/quv.grd'

不能把binary的选项去掉，默认的格式是formatted，楼主要仔细看i/o部分的书啊。
且不说这个，你的输出存放顺序和读取的存放顺序不一致，读出来肯定也不对。fwrite之前想清楚你输出的数据是什么样子的。
open (1,file='e:/lw/quv.grd')
do k=1,z
read(1)((u(i,j,k),i=1,m),j=1,n)
read(1)((v(i,j,k),i=1,m),j=1,n)
end do

传说中的谁

搞不清楚的话就分别输出U,V吧，这样不容易混