bc_gauss function

Apply and subtract a Gaussian smoother in the spectral domain.