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. |