stnmap处理出错

kenever

1. 
   
2. real vec(160)
   
3. open(1,file='C:\Users\Administrator\Desktop\r1601.txt',status='old')
   
4. read(1,*)(vec(i),i=1,160)
   
5. close(1)
   
6. call stntogrd(vec)
   
7. end
   
8. subroutine stntogrd(x)
   
9. real lat(160),lon(160),x(160)
   
10. character*8 stid(160)
    
11. 
    
12. open(2,file='C:\Users\Administrator\Desktop\33.dat')
    
13. do 20 k=1,160
    
14. 20 read(2,*)lat(k),lon(k)
    
15. close(2)
    
16. do 2 i=1,160
    
17. 2 stid(i)=char(i)
    
18. open(3,file='C:\Users\Administrator\Desktop\000.grd',form='binary')
    
19. tim=0.0
    
20. nlev=1
    
21. nflg=1
    
22. do 40 i=1,160
    
23. write(3) stid(i),lat(i),lon(i)
    
24. #,tim,nlev,nflag,x(i)
    
25. 40 continue
    
26. nlev=0
    
27. write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
    
28. close(3)
    
29. return
    
30. end
    

复制代码

根据马红云老师主编的grads教程第八章
有全国160个站点和某月平均降水
将其处理成grd数据
程序为上述程序
数据为
51.72 126.65
48.77 121.92
49.22 119.75
50.45 121.7
49.17 125.23
……
160个站点的经纬度 保存为33.dat
   78   30    7   70   73   87   96   78  108  184   88  135   21   45  153   76  145  206   76   89
   99   34   20   49   33   62   76   48    8    5   23   12   53   92   47   24   62   77   28   60
   48   73   43   75   29   78   25  104  116   53   81   86  118  143   95  135  234  176  298  215
  129  137  206  162  100  190  224  386  340  303  173  330   63  358  239  160  187  151  225  137
  125  188  238  332  243  356  494  419  294  341  401  278  362  462  266  213  132  156  192  199
  208  101   44  101   55  136  158  150   90  115  117  169  304  232  228  178  376  168  185  124
  203  125   75   98  112   10   47   25   57   79   72   75   45   54    8   33   38   34   21    3
   47   44   56  101   13   24   29   10    3    2   25    9    3    2    8   22   14   21   28   45
全国160站某月平均降水 保存为r1601.txt

处理后可以正确出000.grd

之后将编写000.ctl
dset D:\111\000.grd
dtype station
stnmap D:\111\000.map
undef -999.0
title rain
tdef 1 linear Jun1951 1mo
vars 1
p 0 99 rainfall data
endvars

路径正确 运行！stnmap 后输入ctl时出现


我是新手。。请问大神们是哪里错了 处理了一天了 不知道错误出来哪里 非常感谢

zhuzhejl

同问   你找到原因了吗

zhuzhejl

同问   你找到原因了吗 能告诉我一下吗

rongli1015

同问！！！！！！！！！

porpoft

同问＋1

小渚

同问+1

阿坑坑坑阿坑坑

同问+1 求助！！！