求助！用grads处理micaps第三类数据出错！！

diandian

根据baggio的帖子：
micaps第三类数据画图，自动生成ctl和gs文件，造福咯！
http://bbs.06climate.com/forum.php?mod=viewthread&tid=12511&fromuid=41850
(出处: 气象家园)
处理并画出了24小时累计降水量，但与micaps打开的图2相比，中心降水明显偏少（图1），还有为什么高原西部没有降水，micaps打开的图可以看到有降水啊

图1：自动生成的ctl、gs画的图

图2：MICAPS打开的图

图3：各个站点降水量

diandian

所用gs如下：
'reinit'
'open E:\micaps\rain\6\prc5.ctl'
'open E:\micaps\rain\6\data1.ctl'
'set grads off'
'set grid off'
'set parea 1 8.5 1 8.5'
'set mpdset cnworld'
'set xlopts 1 1 0.2'
'set ylopts 1 1 0.2'
'set lon 78 105'
'set lat 25 42'
'set clab forced'

'define a=oacres(g,rr1.2)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set gxout shaded'

'set cmin 0.1'
'set clevs  1 5 10 25 50 '
'e:\micaps\rain\output1.gs'
'd aa'
*'set gxout contour'
*'d aa'
'cbarn'
*'draw title 2011/08/18 08:00rr1'
'printim E:\micaps\rain\6\data1.png white'

碎琼乱玉

{:eb315:}{:eb315:}

pcc5201314

louzhu楼主这个问题 我也遇到过

diandian

pcc5201314 发表于 2015-4-9 21:35
louzhu楼主这个问题 我也遇到过

你解决了么？去掉'define aa=smth9(a1)'9点平滑效果好一点

ymm9106

刚试了，可以出图，对比了micaps数据，没发现错误。