deriveEqIncomp_FromState Subroutine

public subroutine deriveEqIncomp_FromState(state, iField, nElems, varSys, layout, res)

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.

Arguments

TypeIntentOptionalAttributesName
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


Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: QQ
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