This subroutine calculates the flux of the Nernst-Planck equation on the reference cubic element.
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. |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
real(kind=rk), | private | :: | diffusivitySqrt |