GrADS用meteoinfo做mask时出错

zch2010xb

dset E:\MODIS_out\AOD_China_tenyears\AOD_China_tenyears.grd
undef -9999
xdef 31 linear 109 1
ydef 31 linear 31  1
zdef 1  linear 1   1
tdef 1  linear jun1951 1yr
vars 1
AOD 0 1 10yearANG
endvars



DSET E:\MODIS_out\AOD_China_tenyears\northeast.dat
TITLE Mask data
UNDEF -9999.0
XDEF 1000 LINEAR 118.0 0.018018018
YDEF 1000 LINEAR 38.0 0.016016016
ZDEF 1 LINEAR 1.0 1.0
TDEF 0 LINEAR jun1951 1yr
VARS 1
  mask 0 99 background mask data (99)
ENDVARS



'reinit'
'open E:\MODIS_out\AOD_China_tenyears\AOD_China_tenyears.ctl'
'open E:\MODIS_out\AOD_China_tenyears\northeast.ctl'
'enable print E:\MODIS_out\AOD_China_tenyears\AOD_China_tenyears.gmf '
'set grid off'
'set grads off'
'set timelab off'
'set lon 116 135'
'set lat 38 54'
'set gxout grfill'
*'set gxout shaded'
'set mpdset cnworld'
*'set clevs 22 23 24 25 26 27 28'
'set xlopts 1 7 0.2'
'set ylopts 1 7 0.2'
'set xlint 4'
'set ylint 4'
*'E:\beijing1\rb.gs'
*'define XXX=smth9(AOD_year)'
*'d XXX'
'run define_colors.gs'
*'run C:\define_colors.gs'
*'color 0.4 1.8 0.2  -kind blue->mediumblue->white->indianred->red'
'set clevs 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4'
'set ccols 59 55 43 31 23 65 27 69 79'
'set gxout grfill'
'd maskout(AOD,mask)'
'cbarn 1 1 9.5 4.2'

lhaikun@163.com

你的mask文件的分辨率要与你的数据的分辨率相同才行吧，你改改试试

chongzika

lhaikun@163.com 发表于 2015-7-24 11:30
你的mask文件的分辨率要与你的数据的分辨率相同才行吧，你改改试试

用meteoinfo生成数据的时候就需要分辨率一致么？那样子的话边界岂不是太模糊了，不应该吧

lhaikun@163.com

chongzika 发表于 2015-7-24 16:11
用meteoinfo生成数据的时候就需要分辨率一致么？那样子的话边界岂不是太模糊了，不应该吧

找了半天找到这个 你瞅瞅

lhaikun@163.com

回复别人不能加附件 只能重传 抱歉

chongzika

http://bbs.06climate.com/forum.php?mod=viewthread&tid=8221这里有人有同样的问题http://bbs.06climate.com/forum.p ... p;extra=&page=2

balfulosa

一点修改意见，之前用mask的时候参考了论坛里的帖子http://bbs.06climate.com/forum.p ... =893&extra=page%3D1，直接用的类似下面的代码，个人理解也是先需要进行差值才能d变量。然后就是，你的mask变量作为第二个文件打开，mask应该是mask.2，试一试这样好用不。还是建议参考下上面的帖子。还有就是楼上说的，msak和你的AOD变量经纬度范围貌似不一致。
'define aa=oacres(mask.2,AOD,10,5)'
'define bb=smth9(maskout(aa,mask.2))'
'd bb'