Subroutine to calculate the flux for Nernst-Planck equations.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | state(4) | State vector (as conservative variables). The order of this vector has to be where and denoted concentration (scalar) and diffusive fluxes (vector). |
||
real(kind=rk), | intent(in) | :: | diffusivity | The diffusivity in the cell. |
||
real(kind=rk), | intent(out) | :: | flux | The flux inside the cell. The order of this vector is the same as the input arguments. |