There are no bottom boundary conditions for CA and CC in equations 35 from L'Heureux, but py-pde demands left and right boundary conditions. This means that '"curvature" : 0' is a dummy boundary condition. We will make sure that it does not leak into our integration, by using only backward differencing for the spatial derivatives in the right-hand sides of the time derivative equations for CA and CC.
_I will submit a similar PR to merge into the Use_solve_ivp_without_py-pde_wrapper as well, since we want to maintain that branch._
There are no bottom boundary conditions for CA and CC in equations 35 from L'Heureux, but py-pde demands left and right boundary conditions. This means that '"curvature" : 0' is a dummy boundary condition. We will make sure that it does not leak into our integration, by using only backward differencing for the spatial derivatives in the right-hand sides of the time derivative equations for CA and CC.
_I will submit a similar PR to merge into the
Use_solve_ivp_without_py-pde_wrapperas well, since we want to maintain that branch._