台风路径与降水叠加

zshen001

    前段时间发了个台风碧利斯的路径和降水叠加的帖子，http://bbs.06climate.com/forum.p ... mp;page=2#pid704076，后来发现了里面的一些问题，首先是降水与实况不符，原来是gs里面设置新的时间变量以后，没有进行新的插值，还有地形边界以外没有mask掉，现在贴出修改之后的gs和数据文件，希望管理员大大能把我上一篇帖子删掉，以免造成误导，所用的数据来源为micaps第三类降水数据经过处理和台风最佳路径集数据。再贴一下修改前和修改后的图片对比。

修改后

zshen001

原来是'set t 196'
'define ro=oacres(g.1,p.2,50,30,20,10,5)'
'define yy=maskout(ro,ro-0)'
'define zz1=smth9(ro)'
'set t 197'
'define zz2=smth9(ro)'
'set t 198'
'define zz3=smth9(ro)'
'set t 199'
'define zz4=smth9(ro)'
修改后：'set t 196'
'define ro=oacres(g.1,p.2,50,30,20,10,5)'
'define yy=maskout(ro,ro-0)'
'define zz1=smth9(ro)'
'set t 197'
'define ro=oacres(g.1,p.2,50,30,20,10,5)'
'define yy=maskout(ro,ro-0)'
'define zz2=smth9(ro)'
'set t 198'
'define ro=oacres(g.1,p.2,50,30,20,10,5)'
'define yy=maskout(ro,ro-0)'
'define zz3=smth9(ro)'
'set t 199'
'define ro=oacres(g.1,p.2,50,30,20,10,5)'
'define yy=maskout(ro,ro-0)'
'define zz4=smth9(ro)'

yinger

谢谢分享代码！学习了

yinger

代码运行时出现如下错误：


这是什么原因？

yinger

画出来了，跟楼主一样的图形！非常感谢分享代码

yinger

请教一下，grid2006，pre2006，result06.map，这3个是什么文件？怎么做出来的？谢谢

zshen001

yinger 发表于 2016-10-7 18:33
请教一下，grid2006，pre2006，result06.map，这3个是什么文件？怎么做出来的？谢谢

grid2006地图网格化文件，pre2006是2006年降水数据，result06.map是打开pre2006.ctl的索引文件，建议你看下站点数据转格点数据相关资料

zshen001

yinger 发表于 2016-10-7 17:43
代码运行时出现如下错误：

你可能是改了里面的代码吧，正常情况下运行应该不会有什么问题，你再看看路径有没有问题，版本是2.1气象家园整合版

yinger

zshen001 发表于 2016-10-7 20:03
grid2006地图网格化文件，pre2006是2006年降水数据，result06.map是打开pre2006.ctl的索引文件，建议你看 ...

谢谢告知！

yinger

zshen001 发表于 2016-10-7 20:09
你可能是改了里面的代码吧，正常情况下运行应该不会有什么问题，你再看看路径有没有问题，版本是2.1气象 ...

已成功出图，谢谢！