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

diandian

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

所用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'

小二

好的不要不要的

MeteoInfo

插值方法不同造成的。

diandian

小二 发表于 2015-4-9 08:26
好的不要不要的

你知道这个如何解决么？

diandian

MeteoInfo 发表于 2015-4-9 08:46
插值方法不同造成的。

老师你好，可以看到画出的图和micaps图降水带、强度都对应，但唯独中心降水（一个站点达到53mm）强度不够，请问老师如何修改？如何改变插值方法？这方面不了解，请老师指点一二！谢谢！

MeteoInfo

diandian 发表于 2015-4-9 09:54
老师你好，可以看到画出的图和micaps图降水带、强度都对应，但唯独中心降水（一个站点达到53mm）强度不够 ...

GrADS里只有Cressman插值方法，MICAPS中的插值方法我不了解。

在GrADS脚本里用了9点平滑算法，可能对高值有影响，可以去掉试试。

diandian

MeteoInfo 发表于 2015-4-9 11:07
GrADS里只有Cressman插值方法，MICAPS中的插值方法我不了解。

在GrADS脚本里用了9点平滑算法，可能对 ...

感谢老师,去掉9点平滑之后，出图效果提高，micaps里有Cressman、Barnes、Trngle差值方法，一般micaps降水用的是Barnes
E:\micaps\rain\6\data2.png

小二

天哪，我刚开始以为是分享帖呢
- -又上了一课

tr179667208

楼主棒棒哒

烟雨孤舟客

去论坛多搜索，很多都有解答