求解:程序出错

开始起飞1

我先把我的想法说下，txt文件里有一个时间序列的值，因为每年只有一个值，想把这个值铺满区域的每个格点中，得出一时间序列的区域文件，所以想把这个文件转成铺满整个区域格点的二进制文件，然后给其配ctl文件，然后用grads与其他资料作图。
用下面的Fortran程序把一txt文件转成我想要的二进制grd格式文件，但是给grd文件配上ctl后，在grads里输出只有第一年有数据，其余年份全是0，不知道为什么，求解！谢谢

program main
integer,parameter::nx=72, ny=36, na=40
real:: xinjifeng(na,ny,nx)
integer ::year
open(10,file='c:\grads19\data\xinjifeng.txt',form='formatted')
open(20,file='c:\grads19\data\xinjifeng.grd',form='unformatted',status='replace',access='direct',recl=72*36*4)
do i=1,na
read(10,*)year,xinjifeng(i,1,1)
do j=1,ny
   do k=1,nx
     xinjifeng(i,j,k)=xinjifeng(i,1,1)
   end do
  end do
print*,xinjifeng(i,1,1)
write(20,rec=i)((xinjifeng(i,j,k),k=1,nx),j=1,ny)
end do
1000 format('F')
end

txt文件是这样的

river

”每年只有一个值，想把这个值铺满区域的每个格点中，得出一时间序列的区域文件“ 求解释

言深深

你的read语句，为何要设置三维数组，每个都read到xinjifeng(i,1,1)
并且在jk循环的时候，只赋值了一个值，这是何解？还有你配套的ctl文件在哪里？

lqouc

言深深 发表于 2014-4-5 17:46
你的read语句，为何要设置三维数组，每个都read到xinjifeng(i,1,1)
并且在jk循环的时候，只赋值了一个值 ...

他是想让空间里每一个点都是i,1,1的值，得到一个常值空间，我感觉这段程序没啥错。应该自己print一下每一步输出的东西，而且把ctl，gs什么都贴上来。
话说做的什么饭啊？

言深深

lqouc 发表于 2014-4-5 18:34
他是想让空间里每一个点都是i,1,1的值，得到一个常值空间，我感觉这段程序没啥错。应该自己print一下每一 ...

逐步print出来也是一个不错的办法，不过整个程序看下来貌似问题不大···就不知道read的原文件已经描述的ctl文件有没有问题了···

开始起飞1

言深深 发表于 2014-4-5 20:34
逐步print出来也是一个不错的办法，不过整个程序看下来貌似问题不大···就不知道read的原文件已经描述 ...

dset c:\GrADS19\data\xinjifeng.grd
undef  -99.99
title sample grid data
xdef 72 linear -177.5 5.0
ydef 36 linear -87.5 5.0
zdef 1 linear 500 1
tdef 40 linear jan1971 1yr
vars 1
a1 0 99 grid data
endvars

开始起飞1

lqouc 发表于 2014-4-5 18:34
他是想让空间里每一个点都是i,1,1的值，得到一个常值空间，我感觉这段程序没啥错。应该自己print一下每一 ...

dset c:\GrADS19\data\xinjifeng.grd
undef  -99.99
title sample grid data
xdef 72 linear -177.5 5.0
ydef 36 linear -87.5 5.0
zdef 1 linear 500 1
tdef 40 linear jan1971 1yr
vars 1
a1 0 99 grid data
endvars
这是ctl文件

开始起飞1

言深深 发表于 2014-4-5 17:46
你的read语句，为何要设置三维数组，每个都read到xinjifeng(i,1,1)
并且在jk循环的时候，只赋值了一个值 ...

dset c:\GrADS19\data\xinjifeng.grd
undef  -99.99
title sample grid data
xdef 72 linear -177.5 5.0
ydef 36 linear -87.5 5.0
zdef 1 linear 500 1
tdef 40 linear jan1971 1yr
vars 1
a1 0 99 grid data
endvars

言深深

开始起飞1 发表于 2014-4-5 21:50
dset c:\GrADS19\data\xinjifeng.grd
undef  -99.99
title sample grid data

用下面这个程序看看，给你改了一下···
program main
integer,parameter::nx=72, ny=36, na=40
real:: xinjifeng(na,ny,nx)
integer ::year
open(10,file='c:\grads19\data\xinjifeng.txt')
do i=1,na
read(10,*)year,xin
do j=1,ny
   do k=1,nx
     xinjifeng(i,j,k)=xin
   end do
  end do
enddo
close(10)
open(20,file='c:\grads19\data\xinjifeng.grd',form='binary')
write(20)(((xinjifeng(i,j,k),k=1,nx),j=1,ny),i=1,na)
close(20)
end