Fortran读ncep资料转换成的.grd数据，读取不到时次信息

akun

编程小菜求助，将ncep资料用grads提取出变量保存的后缀名是.grd，我用Fortran将其转成txt的时候，程序运行总是提示遇到文件结束。我调试了一下程序，好像是程序读不到时次，我提取的是248个时次的数据，但是Fortran里面设置t=0是才能读取，而且读出来也写不出来。只要t>0，就会遇到文件结束，肿么回事啊？

提数据的gs和Fortran程序如下：
gs:

parameter(nx=12,ny=9,nz=1,nt=0)
dimension sst(nx,ny,nz,nt)                  
real sst

open(10,file='fnl201378_CWATclm.grd',form='binary')
open(20,file='fnl201378_CWATclm.txt')

do it=1,nt
        do iz=1,nz
                read(10)((sst(i,j,iz,it),i=1,nx),j=1,ny)                 
        enddo                                                                                       
enddo

do it=1,nt
        do iz=1,nz
                write(20,*)((sst(i,j,iz,it),i=1,nx),j=1,ny)   
        enddo
enddo

close(10)
close(20)

end

akun

晕，还没打完就点发出去了。
gs文件：
'reinit'
'open E:/zxl/fnl/fnl_201378.ctl'
'set fwrite E:/zxl/fnl/fnl_201378CWATclm.grd'
'set gxout fwrite'

'set lat 26 34'
'set lon 98 109'

tt=1
while(tt<=247)
'set t 'tt

zz=1
while(zz<=1)
'set z 'zz
'd CWATclm'
zz=zz+1
endwhile

tt=tt+1
endwhile
'disable fwrite'
;




程序：

parameter(nx=12,ny=9,nz=1,nt=0)
dimension sst(nx,ny,nz,nt)                  
real sst

open(10,file='fnl201378_CWATclm.grd',form='binary')
open(20,file='fnl201378_CWATclm.txt')

do it=1,nt
        do iz=1,nz
                read(10)((sst(i,j,iz,it),i=1,nx),j=1,ny)                 
        enddo                                                                                       
enddo

do it=1,nt
        do iz=1,nz
                write(20,*)((sst(i,j,iz,it),i=1,nx),j=1,ny)   
        enddo
enddo

close(10)
close(20)

end

言深深

将程序段：
do it=1,nt
        do iz=1,nz
                read(10)((sst(i,j,iz,it),i=1,nx),j=1,ny)                 
        enddo                                                                                       
enddo

修改为如下形式试试：


               read(10)((((sst(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)

akun

言深深 发表于 2014-3-9 20:50

                               
登录/注册后可看大图

将程序段：
do it=1,nt
        do iz=1,nz

还是不行。

river

akun 发表于 2014-3-9 20:45

                               
登录/注册后可看大图

晕，还没打完就点发出去了。
gs文件：
'reinit'

提取的资料你确定是对的吗？有没有画一下图试试呢。还有就是把set lat 和set lon 改成用set x 和set y，并且放在循环里面。