如何消除作差值后画图出现的由地形导致的颜色异常区域

13765676798

各位老师好！我在计算变量的时候设置了缺省值，可是打开两个文件做差值时右上角出现了异常的颜色块，请问如何在gs文件中消除右上角异常部分呢？请赐教，感激不尽！
部分gs文件如下：
'reinit'

'open G:\bs2\hua\hisvv_1.ctl'
'open G:\bs2\hua\vv_1.ctl'
'set parea 1 5.3 3.5 6.5'
'set gxout shaded'

'set rgb 21 0 0 133'
'set rgb 22 0 0 153'
'set rgb 23 0 0 174'
'set rgb 24 0 0 194'
'set rgb 25 0 0 214'
'set rgb 26 0 0 235'
'set rgb 27 0 0 255'
'set rgb 28 20 20 255'
'set rgb 29 41 41 255'
'set rgb 30 82 82 255'
'set rgb 31 102 102 255'
'set rgb 32 122 122 255'
'set rgb 33 163 163 255'
'set rgb 34 204 204 255'
'set rgb 35 224 224 255'
'set rgb 36 245 245 255'

'set rgb 37 255 255 255'

'set rgb 38 255 245 245'
'set rgb 39 255 224 224'
'set rgb 40 255 204 204'
'set rgb 41 255 163 163'
'set rgb 42 255 122 122'
'set rgb 43 255 102 102'
'set rgb 44 255 82 82'
'set rgb 45 255 41 41'
'set rgb 46 255 20 20'
'set rgb 47 255 0 0'
'set rgb 48 235 0 0'
'set rgb 49 214 0 0'
'set rgb 50 194 0 0'
'set rgb 51 174 0 0'
'set rgb 52 153 0 0'
'set ccols 35 31 33 35 40 41 43 45'
'set clevs -30 -20 -10 0 10 20 30'

'd vv.2-vv'
'draw title a.Win'

river

我认为正常情况下，GrADS会自动略过缺测值，不会参加计算的吧

13765676798

river 发表于 2020-5-3 19:04
我认为正常情况下，GrADS会自动略过缺测值，不会参加计算的吧

可是出现了这种情况老师有没有觉得可行的解决方法呢？谢谢

river

13765676798 发表于 2020-5-3 20:46
可是出现了这种情况老师有没有觉得可行的解决方法呢？谢谢

你的两个资料应该缺测值是一样的吧？两个资料是能见度资料？两个资料是时间上不同？建议把资料的描述文件也贴上

13765676798

对的，两个资料是850hPa滤波后的北太平洋的经向风方差，缺测值应该是一样的，都是一个模式输出的历史和未来的资料

13765676798

ctl文件是这样的：
dset G:\bs2\MRI-CGCM3-0000\rcp85\vv_850\vv_4.dat
title Daily NOAA-CIRES 20th Century Reanalysis
undef -9.96921e+36
xdef 49 linear 120.000000 2.500000
ydef 17 linear 20.000000 2.500000
tdef 1 linear 00Z01Jan2075 1dy
zdef 1 levels 850
vars 1
vv 1 99
endvars

278803532

有可能是你的缺测的描述不对，undef -9.96921e+36，就是这么一句，虽然你描述了缺测值，但是你图中的缺测的值不是这个。自己找出来从新定义就好了