关于mask和区域面积加权的顺序的问题

旺仔女孩

请问大家有没有过数据先mask了陆地数据，然后再进行区域加权平均之后数据都是缺省值的情况，可以请大家帮我看看是程序的问题吗，还是区域加权平均只能在mask之前进行呢
;.......滤波...........
e_TS!0="member"
e_TS@missing_value=e_TS@_FillValue
wgts_decade = (/1,6,19,42,71,96,106,96,71,42,19,6,1/)*1.0
wgts_decade = wgts_decade/sum(wgts_decade)
kopt        = 1  
e_TS_lp= wgt_runave_n_Wrap(e_TS, wgts_decade, kopt,1)
;.......mask...........
c             = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r")
lsdata        = c->LSMASK
lsm           = landsea_mask(lsdata,TS&lat,TS&lon)
e_TS_mask= mask(e_TS_lp,lsm.eq.1,False)
copy_VarMeta(e_TS,e_TS_mask)

TSNP=e_TS_mask(:,:,{-60:70},{110:290})
;......................
;.......trend..........
ea_TS=dim_avg_n_Wrap(TSNP,0)
;.......区域加权平均........
lat = f ->lat({-60:70})
lon = f ->lon({110:290})
jlat1  = dimsizes(lat)
rad1   = 4.0*atan(1.0)/180.0
re1    = 6371220
rr1    = re1*rad1
dlon1  = abs(lon(2)-lon(1))*rr1
dx1    = dlon1*cos(lat*rad1)
dy1    = new(jlat1,typeof(dx1))
dy1(0) = abs(lat(2) - lat(1)) * rr1
dy1(1:jlat1-2) = abs(lat(2:jlat1-1) - lat(0:jlat1-3)) * rr1 * 0.5
dy1(jlat1-1)   = abs(lat(jlat1 - 1) - lat(jlat1 - 2)) * rr1
area1 = dx1 * dy1
e_TSNP1 = wgt_areaave(ea_TS, area1, 1.0, 1)
printMinMax(e_TSNP1,1)
然后e_TSNP1输出之后min=9.96921e+36   max=9.96921e+36

sususir

opt
If opt = 0, the area average is calculated using available non-missing data. If opt = 1, then if any point in q is missing, the area average is not computed. In this case, it will be set to the missing value, which is indicated by q@_FillValue, or the default missing value if q@_FillValue is not set.

旺仔女孩

sususir 发表于 2021-3-16 14:18
opt
If opt = 0, the area average is calculated using available non-missing data. If opt = 1, then i ...

谢谢前辈回答我的问题，问题已经解决了，我没有仔细看官网上的说明，下次会注意