站点数据插值到格点上

冷水舟

我按照马红云李丽萍老师那本书第八章的 步骤来做的   只不过换成了我自己的数据和区域经纬度，插值成2x2的网格~然后最后一步RUN那个gs文件的时候就出现问题了提示如下

不知道是不是fortran程序出错了~我的fortran程序如下
real vec(19)
cccc 从文件"51-st.txt"中读入数据 cccc
     open（1,file='I:\51-st.txt',status='new'）
     read(1,*) (vec(i),i=1,19)
     close(1)
cccc 将文本记录转换为GrADS所支持的二进制记录 cccc
     call stntogrd(x)
  end
cccc stntogrd子程序 cccc
     subroutine stntogrd(x)
  real lat(19),lon(19),x(19)
  character*8 stid(19)
cccc 从文件"china.dat"中读入19站的经纬度 cccc
     open(2,file='I:\china.dat')
  do 20 k=1,19
20   read(2,*) lat(k),lon(k)
     close(2)
cccc 给出各站的站号 cccc
     do 2 i=1,19
2  stid(i)=char(i)
cccc 按顺序向二进制文件“51-st.grd”中写入记录 cccc
     open (3,file='I:\51-st.grd',form='binary')
  tim=0.0
  nlev=1
  nflag=1
  do 40 i=1,19
  write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
               
40       continue
cccc 在文件最后，给出一个特别记录，表示这个时间组的记录结束ccccc
      nlev=0
      write(3) stid(i-1),lat（i-1），lon(i-1),time,nlev,nflag
      close(3)
      return
      end
大家帮我看看为什么出错了~~~谢谢~谢谢~

liyf

http://bbs.06climate.com/forum.php?mod=viewthread&tid=4903

船头尺

好资料，多谢分享！