This routine computes equilibirium from state array Here it is assumed that rho0 = 1.0
This subroutine's interface must match the abstract interface definition derive_FromState in derived/mus_derVarPos_module.f90 in order to be callable via velFromState, equilFromState, momFromState, velocitiesFromState, and momentaFromState function pointers.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=rk), | intent(in) | :: | state(:) | Array of state n * layout%fStencil%QQ * nFields |
||
| integer, | intent(in) | :: | iField | Current field |
||
| integer, | intent(in) | :: | nElems | number of elements |
||
| type(tem_varSys_type), | intent(in) | :: | varSys | variable system which is required to access fieldProp information via variable method data c_ptr |
||
| type(mus_scheme_layout_type), | intent(in) | :: | layout | scheme layout contains stencil definition and lattice weights |
||
| real(kind=rk), | intent(out) | :: | res(:) | Output of this routine Dimension: n * nComponents of res |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | private | :: | |||||
| integer, | private | :: | iElem | ||||
| integer, | private | :: | iDir | ||||
| real(kind=rk), | private | :: | f(layout%fStencil%QQ) | ||||
| real(kind=rk), | private | :: | rho | ||||
| real(kind=rk), | private | :: | vel(3) | ||||
| real(kind=rk), | private | :: | usq | ||||
| real(kind=rk), | private | :: | ucx |