Hello everyone,
I have a monthly sea ice time series of 30 years that has three different north pole mask areas (see the purple shaded areas in the attached figure).
I would like to use the mask in the upper right for the whole time period, i.e., the mask area in the lower left will be expanded while that in the upper left will be shrank. The previous mask areas in or outside of this mask will be considered/filled as missing data.
For redefining the mask part, the relevant code is as follows:
let fac=1.e-5
let NPhole=-0.05
use "`root_dir`/`dta_dir`/`var_id`_`dta_typ`_`var_typ`_`hemi_id`_y`yrs`.des"
let sic_GSFC=`var`[d=1]
define grid/like=sic_GSFC grd_GSFC
