mus_calcAuxField_fluidIncomp_d2q9 Subroutine

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

This routine compute auxFields density and velocity for incompressible d2q9 model for fluid and nGhostFromCoarser elements


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 :: dens_pos
integer, private :: vel_pos(3)
integer, private :: elemOff
real(kind=rk), private :: rho
real(kind=rk), private :: pdf(9)
real(kind=rk), private :: u_x
real(kind=rk), private :: u_y
integer, private :: iElem
integer, private :: QQ
integer, private :: nScalars