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

vera

我想画夏天850的平均风场和位势高度场，用楼上的方法生成了xx.dat、xx.ctl文件，gs文件如下：
'reinit'
'sdfopen e:\uwnd.mon.mean.nc'
'sdfopen e:\vwnd.mon.mean.nc'
'sdfopen e:\hgt.mon.mean.nc'
'open D:\grads\tibetmaskout.ctl'

'enable print e:\1999summer850.gmf'

'set mpdset cnworld cnriver'
'set map 15 2 5'
'set lat 20 50'
'set lon 90 130'
'set xlint 5'
'set ylint 5'
'set xlopts 1 2 0.2'
'set ylopts 1 2 0.2'

'set gxout shaded'      
'set clevs 1500'      
'set ccols 0 5 0'
'set grads off'
'define hgtave=ave(hgt.3,t=618,t=620)'
'd maskout(hgtave,mask.4)'

'set z 3'
'set gxout vector'
'set arrowhead 0.15'
'set ccolor 1'
'set arrscl 0.3'
'define uwndave=ave(uwnd.1,t=618,t=620)'
'define vwndave=ave(vwnd.2,t=618,t=620)'
'd maskout(uwndave,mask.4);maskout(vwndave,mask.4)'

'set font 2 '
'set string 2 c 6 0'
'set strsiz 0.3 0.2'
'draw string 4.46 2.77 D'
'draw string 3.73 5.79 D'
*'draw string 4.25 2.35 D'

'set string 4 c 6 0'
'draw string 6.7 4.8 G'


'set font 0 '
'set string 1 c 2 0'
'set strsiz 0.2 0.3'
'draw string 5.5 8.0 1999 summer 850hPa Heights and Wind'

'print'
'disable print'

提示的错误信息是Operation error:  Incompatable grids
  Dimension ranges aren't equivalent
Operation Error:  Error from maskout function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = maskout(hgtave,mask.4)


可否请楼主指点一下？

mofangbao

vera 发表于 2011-8-17 20:32

                               
登录/注册后可看大图

我想画夏天850的平均风场和位势高度场，用楼上的方法生成了xx.dat、xx.ctl文件，gs文件如下：
'reinit'
' ...

说出你的问题吧，直接贴个gs都不知道你的问题在哪里啊，如果是自己没有思考过这个gs还是要先思考的啦，自己先努力哦

MeteoInfo

错误提示里已经写得很清楚了，maskout和数据网格不吻合

笨笨不笨

MeteoInfo 发表于 2011-8-17 23:19

                               
登录/注册后可看大图

错误提示里已经写得很清楚了，maskout和数据网格不吻合

借个地方，表达一下长期以来对您的膜拜

enso816

学会啦，谢谢！

luoziwuhui

学习一下啊！

以父之名

好高兴啊，来到这个网站，像掉进宝藏里一样，到处都是财富。非常支持版主，以后我也要多跟大家共享

张世轩

好厉害，楼主很强大

冬日情愫+

新手进来学习学习！！

Forcast

这个论坛很强大，多谢版主