站点转成格点画图遇到麻烦了 求大神帮帮忙！

立志天下

我的数据是85个站点的平均气温资料 站点转换成格点已经成功，但是在最后一步画图出现问题。
这是画图的gs 帮我看看

1. 'reinit'
   
2. 'open e:\cc\fyw\grid.ctl'
   
3. 'open e:\cc\fyw\t.ctl'
   
4. 'set grads off'
   
5. 
   
6. 'set lon 73 140'
   
7. 'set lat 17 55'
   
8. 'set grads off'
   
9. 
   
10. 'set mpdset cn'
    
11. 'define a=oacres(p.1,t.2)'
    
12. 
    
13. 'define aa=smth9(a)'
    
14. 'd aa'
    
15. 
    
16. 'printim e:\cc\fyw\t.png white'
    
17. ;
    

复制代码

画图出现的问题

立志天下

我的数据是85个站点的平均气温资料 站点转换成格点已经成功，但是在最后一步画图出现问题。
这是画图的gs 帮我看看

1. 'reinit'
   
2. 'open e:\cc\fyw\grid.ctl'
   
3. 'open e:\cc\fyw\t.ctl'
   
4. 'set grads off'
   
5. 
   
6. 'set lon 73 140'
   
7. 'set lat 17 55'
   
8. 'set grads off'
   
9. 
   
10. 'set mpdset cn'
    
11. 'define a=oacres(p.1,t.2)'
    
12. 
    
13. 'define aa=smth9(a)'
    
14. 'd aa'
    
15. 
    
16. 'printim e:\cc\fyw\t.png white'
    
17. ;
    

复制代码

画图出现的问题

立志天下

在贴上我的ctl

1. DSET e:\cc\fyw\ave-t_7.grd
   
2. DTYPE STATION
   
3. STNMAP e:\cc\fyw\ave-t_7.map
   
4. UNDEF 9999.9
   
5. TITLE MONTHLY tmp
   
6. TDEF 1 LINEAR JUL1950 1yr
   
7. VARS 1
   
8. t 0 99 MONTHLY tmp DATA
   
9. ENDVARS
   
10. 
    

复制代码

river

另一个ctl 呢？

lqouc

这样的问题，多半是Fortran就不对，map文件生成时是不是就说了只有一个站点啊？

立志天下

river 发表于 2014-3-7 11:51

                               
登录/注册后可看大图

另一个ctl 呢？

这个 ctl 没问题的吧

1. dset  e:\cc\fyw\grid.grd
   
2. title  grid
   
3. undef 99999.9
   
4. xdef 71 linear 70.0 1.0
   
5. ydef 41 linear 15.0  1.0
   
6. zdef 1  levels 1000   
   
7. tdef    1 linear JUN1961 1mo
   
8. vars 1
   
9. p  0 99 grid data1
   
10. endvars
    

复制代码

立志天下

lqouc 发表于 2014-3-7 12:41

                               
登录/注册后可看大图

这样的问题，多半是Fortran就不对，map文件生成时是不是就说了只有一个站点啊？

生成map 是这样的，请问哪里出问题了吗？

lqouc

立志天下 发表于 2014-3-7 13:36

                               
登录/注册后可看大图

生成map 是这样的，请问哪里出问题了吗？

看上去问题不大，话说你的两个ctl里面tdef的参数最好一致。

立志天下

lqouc 发表于 2014-3-7 13:40

                               
登录/注册后可看大图

看上去问题不大，话说你的两个ctl里面tdef的参数最好一致。

版主  果然是时间这个没设置好，万分感谢！现在 弄出来了。谢谢！！

river

立志天下 发表于 2014-3-7 13:30

                               
登录/注册后可看大图

这个 ctl 没问题的吧

楼下的已经说过了，问题就是两个ctl 时间设置得一致