How do I set the proper boundary conditions in solve_bvp?

Question

I am trying to solve the following boundary value problem:

D*(dS^2/dz^2) - v*(dS/dz) - S/(S+k_s) = 0

Where D, v and k are constants (see code below). The bound