compass.ocean.tests.tides.dem.dem_pixel.blend_front

compass.ocean.tests.tides.dem.dem_pixel.blend_front(e1st, i1st, e2nd, halo, sdev)[source]

Create a mask of linear weights to ‘blend’ two elev. fun at the ice sheet/shelf front:

ELEV = (1.-MASK) * E1ST + (0.+MASK) * E2ND

Elev. data assoc. with the 1st data-set is preserved in pixels where ice-thickness is non-zero. The two datasets are blended over a distance of approx. HALO pixels. The mask is additionally smoothed via a Gaussian filter with standard-deviation sigma=SDEV.