code: def gaussfill(dem, size=3, newmask=None):
      
      smooth = gauss_fltr_astropy(dem, size=size)
      smooth[~dem.mask] = dem[~dem.mask]
      if newmask is not None:
          smooth = np.ma.array(smooth, mask=newmask)
      return smooth
docs: Gaussian filter with filling