grads读取多个数据

wuyang0074

各位前辈好：
我在做grads读取多个数据，下面是脚本，请帮忙给看看，哪里有错误？谢谢了！
报错：
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
……同上，好多！…………………………………………………………
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Data Request Warning:  Request is completely outside file limits
Cannot contour grid - all undefined values
Cannot contour grid - all undefined values
                                                  脚本如下：
'clear'
'reinit'

'open 2000_09.ctl'
!'open 2005_09.ctl'

'set lon 86 127'
'set lat 37 54'
'set grads off'
'set timelab off'
'set grid off'
'set lev 1'
'set mpdraw off'
'set xlint 10'
'set ylint 5'

'set gxout shaded'
'set csmooth on'
'set clevs 0.01 0.05 0.09 0.13 0.17 0.21 0.25 0.29 0.33 0.37 0.41 0.45 0.49 0.53 0.57 0.61 0.65 0.69 0.74 0.8 0.9 1'
'set csmooth on'
'set t s'
'a=ave(smois.1,t=1,t=117)'
'd a'


!'open 2000_09.ctl'
'open 2005_09.ctl'

'set lon 86 127'
'set lat 37 54'
'set grads off'
'set timelab off'
'set grid off'
'set lev 1'
'set mpdraw off'
'set xlint 10'
'set ylint 5'

'set gxout shaded'
'set csmooth on'
'set clevs 0.01 0.05 0.09 0.13 0.17 0.21 0.25 0.29 0.33 0.37 0.41 0.45 0.49 0.53 0.57 0.61 0.65 0.69 0.74 0.8 0.9 1'
'set csmooth on'
'set t s'
'b=ave(smois.2,t=1,t=117)'
'd b'

'define h=(a+b)/2'
'd h'
'set strsiz 0.2 0.2'

'draw shp MGGY.shp'
'set csmooth on'
'draw string 0.7 7 (2015-10)'
'cbarn 0.7 0 5.5 1.2'

小其其格

'set lev 1'   ***用lev的话要设置为气压层1000hPa或850hPa或700hPa等……用set z 1才设置为1、2、3……

还有你'set gxout shaded'绘制a、b、h在一张图上，那最后a、b都会被'd h'覆盖

wuyang0074

小其其格 发表于 2017-10-16 06:42
'set lev 1'   ***用lev的话要设置为气压层1000hPa或850hPa或700hPa等……用set z 1才设置为1、2、3……

...

我只要h就行了，您看我这个整体上写的对吗？

小其其格

你试试看下面这样能达到你要求么：

1. 'clear'
   
2. 'reinit'
   
3. 
   
4. 'open 2000_09.ctl'
   
5. 'open 2005_09.ctl'
   
6. 
   
7. 'set lon 86 127'
   
8. 'set lat 37 54'
   
9. 'set grads off'
   
10. 'set timelab off'
    
11. 'set grid off'
    
12. 'set lev 1000'
    
13. 'set mpdraw off'
    
14. 'set xlint 10'
    
15. 'set ylint 5'
    
16. 
    
17. 
    
18. 'a=ave(smois.1,t=1,t=117)'
    
19. 'b=ave(smois.2,t=1,t=117)'
    
20. 
    
21. 'define h=(a+b)/2'
    
22. 'd h'
    
23. 'set strsiz 0.2 0.2'
    
24. 
    
25. 'draw shp MGGY.shp'
    
26. 'set csmooth on'
    
27. 'draw string 0.7 7 (2015-10)'
    
28. 'cbarn 0.7 0 5.5 1.2'

复制代码

wuyang0074

小其其格 发表于 2017-10-16 17:24
你试试看下面这样能达到你要求么：

还是不行，您看，能不能加个QQ，这样能快点1506887954