maskout总是出错

我陪你一世流年

用maskout这个功能总是出错，请问是什么原因呢？我用的是2.5*2.5的ncep资料，用meteorInfo制作的maskout文件也是从105°E-110°E，25°N-30°N。  谢谢啦~~~






这是gs:
'reinit'
'open D:\qdn.ctl'
'sdfopen d:\grads\data\shum.2013.nc'
'set mpdset qdnbj'
t=711
'set lev 925'
'set t 't
'set lat 25 30'
'set lon 105 115'
'set ccolor 1'
'define a1=maskout(shum.2,mask.1)'
'define aa=smth9(a1)'
'd aa'
'draw shp qdnbj.shp'
'printim d:\grads\shum\shum.png white'
'c'





这是修改后的ctl：

DSET d:\qdn.dat
TITLE Mask data
UNDEF -9999.0
XDEF 5 LINEAR 105.0 2.5
YDEF 3 LINEAR 25.0 2.5
ZDEF 1 LINEAR 1.0 1.0
TDEF 1 LINEAR 12:00Z17Jun2013 6hr
VARS 1
  mask 0 99 background mask data (99)
ENDVARS

edwardli

历史太悠久了。。。我看主要矛盾（错误）是没有打开文件
最应该修改的是\，为/
再者，mask的TDEF应该和主数据一致
先试试吧，历史太悠久了，三四年不用GrADS了……

我陪你一世流年

edwardli 发表于 2021-8-17 21:08
历史太悠久了。。。我看主要矛盾（错误）是没有打开文件
最应该修改的是\，为/
再者，mask的TDEF应该和主 ...

我nc资料设置的时间是2013年6月27日20时，所以地图ctl为 TDEF 1 LINEAR 12:00Z17Jun2013 6hr  不知道对不对？斜杠貌似没有问题，我画其他图的时候也是用的\

四叶草

不用maskout好多年了，目前直接用province-basemap.gs

edwardli

我陪你一世流年 发表于 2021-8-18 09:12
我nc资料设置的时间是2013年6月27日20时，所以地图ctl为 TDEF 1 LINEAR 12:00Z17Jun2013 6hr  不知道对不 ...

那t 为什么是711……mask数据里面可没有711个时次啊

我陪你一世流年

edwardli 发表于 2021-8-18 13:02
那t 为什么是711……mask数据里面可没有711个时次啊

t=711 就是2013年6月27日20时  所以TDEF 1 LINEAR 12:00Z17Jun2013     不知道对不对

edwardli

我陪你一世流年 发表于 2021-8-19 14:15
t=711 就是2013年6月27日20时  所以TDEF 1 LINEAR 12:00Z17Jun2013     不知道对不对

在脚本中，t为711，那么对于任何一个打开的文件都会去找t=711，而不是找t＝2017062720
还是那句话，把mask描述文件的TDEF 从数据文件的ct中直接复制过来l

lightmoon

edwardli 发表于 2021-8-19 16:20
在脚本中，t为711，那么对于任何一个打开的文件都会去找t=711，而不是找t＝2017062720
还是那句话，把mas ...

GrADS对于时间的判断是这样的，t=711表示在当前文件的时间格式下（此处默认为qdn），找到所对应的绝对时间，也就是12Z01jun2013+710*6hr，然后对于其他所有文件，当前时间为该绝对时间。而文件qdn只在t=1时存在数据，所以画不出来，我觉得大概是这样的。

river

'define a1=maskout(shum.2,mask.1(t=1))'
另外确认一下那个路径下 是不是确实存在qdnbj 的底图文件