Dear Georg,
when you set Dirichlet boundary conditions the equation gets replaced by
the Dirichlet boundary condition. I think that's what's clear.
When you set Neumann boundary conditions in Dumux you set a value for
the entire flux term (advective + diffusive part). So when you set a
Hi Georg,
If you have grad(S) somewhere in your equation yes. Otherwise no.
If you want to have a boundary condition like grad(S)*n = 0, calculate the
residual (the whole flux term, i.e. everything with a divergence in front) in
solDependentNeumann and just leave out the part of the flux that
