求助，关于站点数据处理无法生成.map文件

liusun2010

参照坛子里《GrADS站点文件作图详细解决方案》，处理站点数据，原数据为123个站一段时间的累积降水，信息如下：经度，纬度，降水量，站点编号
fortran程序如下：
program sta2grd
implicit none
  character(len=8) stid
  real lon,lat,rain,tim
  integer nlev,flag

  tim=0.0
  nlev=1
  flag=1
  open(1,file='rr_2019,08,10,08_2019,08,12,13.txt',status='old')
  open(2,file='sta.grd',status='replace',form='binary')
   read(1,*)
   read(1,*)
   10 read(1,*)lon,lat,rain,stid
   print*,stid,lat,lon,rain
   write(2)stid,lat,lon,tim,nlev,flag,rain
   goto 10
   100 continue
   close(1)
   nlev=0
   write(2)stid,lat,lon,tim,nlev,flag
   close(2)
   end

生成了相应的sta.grd文件，对应的ctl为：
DSET G:/sta/sta.grd
DTYPE station
STNMAP G:/sta/sta.map
UNDEF 9999
TITLE Rain Data
TDEF 1 linear 29mar2019 1mo
VARS 1
rain 0 99 Rain Data
ENDVARS


grads生成.map文件出错，具体为：
ga-> !stnmap -i g:/sta/sta.ctl
  Name of binary data set: G:/sta/sta.grd
  Number of times in the data set: 1
  Number of surface variables: 1
  Number of level dependent variables: 0

Starting scan of station data binary file.
Binary data file open: G:/2019Lekima/sta/sta.grd

Processing time step 1
  Low Level I/O Error:  Read error on data file
  Error reading 28 bytes at location 3936
  Possible cause: Premature EOF
ga->


想问下是否有人知道是为啥，感谢

liusun2010

解决了，10 read(1,*,end=100)lon,lat,rain,stid改成这样就好了