Added target boundaries to bed and surface, which fixed non-convergence of Stokes solver