A two-stage numerical approach for the sparse initial source identification of a diffusion-advection equation