mus_calcAuxField_MS Subroutine

private subroutine mus_calcAuxField_MS(auxField, state, neigh, nSize, nSolve, iLevel, stencil, varSys, derVarPos)

This routine compute auxField density and momentum for each species for multicomponent models. The momentum computed here is only momentum of transformed PDF. The momentum of original PDF is computed by solving linear equation system in compute kernel and the momentum in auxField is updated there.


real(kind=rk), intent(inout) :: auxField(:)

output auxField array

real(kind=rk), intent(in) :: state(:)

input state array

integer, intent(in) :: neigh(:)

connectivity array

integer, intent(in) :: nSize

number of elements in the state array

integer, intent(in) :: nSolve

number of elements excluding halos

integer, intent(in) :: iLevel

current level

type(tem_stencilHeader_type), intent(in) :: stencil

stencil header

type(tem_varSys_type), intent(in) :: varSys

variable system definition

type(mus_derVarPos_type), intent(in) :: derVarPos(:)

position of derived quantities in varsys



integer, private :: iElem
integer, private :: iFld
integer, private :: iDir
integer, private :: elemOff
integer, private :: nScalars
integer, private :: QQ
integer, private :: nFields
integer, private :: dens_pos
integer, private :: mom_pos(3)
real(kind=rk), private :: pdf(stencil%QQ)
type(mus_varSys_data_type), private, pointer:: fPtr
type(mus_scheme_type), private, pointer:: scheme