江湖救急

20131301013

我的还是不行，为什么呢？求大神帮忙
program ex03
implicit none
integer i,j
real lat(160),lon(160),r1607(160,63)
character*8 stid(160)
real tim,nlev,nflag
open (1,file='d:\grads\grads\r1607.txt',status='old')
read(1,*)((r1607(i,j),i=1,160),j=1,63)
close(1)
open(2,file='d:\grads\grads\id.txt',status='old')
do i=1,160
read(2,*)stid(i),lat(i),lon(i)
enddo
close(2)
open (3,file='d:\grads\grads\rain.grd',form='binary',status='replace')
do j=1,63
tim=0.0
nlev=1
nflag=1
do i=1,160
write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,r1607(i,j)
enddo
nlev=0
nflag=0
write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
enddo
close(3)
end

dset d:/grads/grads/rain.grd
dtype station
stnmap d:/grads/grads/rain.map
undef -999.0
tdef 63 linear jul1951 1yr
vars 1
rain 0 99 rainfall data
endvars

qixiang_learner

去掉nflag=0试试

20131301013

qixiang_learner 发表于 2015-7-6 23:54
去掉nflag=0试试

还是不行