grads加fortran画站点年均降水图

yqy1204

Enter stn ctl filename:  E:\taizhan\temp\zhandian\111p.ctl
  Name of binary data set: E:\taizhan\temp\zhandian\111p.grd
  Number of times in the data set: 50
  Number of surface variables: 1
  Number of level dependent variables: 0
Starting scan of station data binary file.
Processing time = 1
Binary data file open: E:\taizhan\temp\zhandian\111p.grd
  Invalid station hdr found in station binary file
  Possible causes:  Invalid level count in hdr
                    Descriptor file mismatch
                    File not station data
                    Invalid relative time
    levs = 1  flag = 1  time = 1961
  File name = E:\taizhan\temp\zhandian\111p.grd
这是在grads运行ctl的结果，说是二进制文件有问题。下面是原始txt被转为二进制的fortran程序
program sta2grd
implicit none
character*8 stid
real :: lon,lat,pre,year
integer :: nlev=1
integer :: flag=1
integer i
open(10,file='E:\taizhan\temp\zhandian\111p.txt',status='old',action='read')
open(11,file='E:\taizhan\temp\zhandian\111p.grd',status='replace',action='write',access='direct',form='unformatted',&
        recl=2500*5)

do i=1,2500
read(10,*)  stid,lat,lon,year,pre
!if(i==10 .or. i==11) then
write(11,rec=i) stid,lat,lon,year,nlev,flag,pre
!write(*,*) stid,lat,lon,tim,nlev,flag,rain
!endif
enddo
close(10)
nlev=0
write(11,rec=i+1) stid,lat,lon,year,nlev,flag,pre
close(11)
end program
求指点，到底是哪里错了？

男紫汗

不是有好几种可能，描述性文件如何？

yqy1204

男紫汗 发表于 2017-5-1 10:26
不是有好几种可能，描述性文件如何？

C:\Users\la\Desktop\11.png描述性文件
dset E:\taizhan\temp\zhandian\111p.grd
options template   
title 111p data
dtype station
stnmap E:\taizhan\temp\zhandian\111p.map
undef 1e+20
xdef 57 linear  70 2.5
ydef 17 linear  15 2.5
zdef 1 linear 1000 1
tdef 50 linear jan1961 1yr
vars 1
p   0  99  pre
endvars

我刚才打了打开dat看居然只有一行

男紫汗

照着这个在检查检查

蔷薇花烬

楼主最后问题解决了吗？同样遇到了这个问题~