maskout方法的使用以及maskout文件制作

xuchao51

谢谢分享！学习

isalcj

很有用的帖子，支持！！

blui

不错哦。。。

blui

画站点图画得漂亮还是很难的。。。。

blui

为什么我用站点图叠加风场画画不到的。。

wc619lf

感谢楼主无私的分享

学而无术

学习了，谢谢

北五味

学习了，感谢大神的指点啊，新人有活路了，呜呜~~~~

平流层的萝卜

多谢楼主分享！但我把多个省份的底图输出后却出现了错误，求指导：
这个是输出的maskout的ctl描述文件：
DSET E:\gdd\plot\huabei.dat
TITLE Mask data
UNDEF -9999
XDEF 20 LINEAR 110 0.6842105
YDEF 20 LINEAR 31 0.6315789
ZDEF 1 LINEAR 1 1
tdef 248 linear  02z01jan2013 3hr
VARS 1
mask 0 99  background mask data
ENDVARS
我的站点数据描述文件如下：
dset  E:\gdd\plot\vis.grd
dtype    station
stnmap   E:\gdd\plot\vis.map
undef    -999.0
title   visibility
tdef 248 linear  02z01jan2013 3hr
vars 1
vi  0 99  visiblity
endvars
最终用来出图的gs描述文件如下：
'open E:\gdd\plot\vis.ctl'
'open E:\gdd\plot\huabei.ctl'
'enable print E:\gdd\plot\micaps_vis.gmf'

i=1
while(i<=248)
'set grads off'
'set lon 100 130'
'set lat 20 50'
'set lev 850'

'set t 'i''
'set mpdset cnriver '
'set map 1 1 9'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'

*'define a=oacres(mask.2,vi.1)'
*'define a1=maskout(a,g.2-0.5)'
*'define aa=smth9(a1)'

'set gxout shaded'
'set clevs 0.1 0.3 1 10 15 20 25 30'
'd  maskout(oacres(mask.2,vi.1),mask.2)'
*'d a'
*'d vi.1'
*'cbarn.gs'

'q time'
x=subwrd(result,3)
'draw title 'x''

'print'
'c'
i=i+1
endwhile
'disable print'
'reinit'
;
出现的错误如图：

出图错误

一共要出248张图，但是只出了第一张图，后边全是空白，这是为什么呢？
第一张图应该是正确的。

第一张图

站点数据都没问题，之前用1*1的格点插值都是正常出图。
还请指教，多谢啦！

小小人

可是，有些缺测的地方会出的图不算美观。不过还是谢谢版主。