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
let SSMI_mskNP=IF (abs(sic_GSFC[l=180]-NPhole) LE fac) THEN 1 ELSE 0
set var/bad=0 SSMI_mskNP
let conc_GSFC = IF (SSMI_mskNP[g=grd_GSFC@asn] EQ 1) THEN (`NPhole`) ELSE sic_GSFC
However, the north pole hole area in conc_GSFC is not expanded in the later part of the record as it should.
What did I do wrong?
Thanks,
--- Peng
--
