mrt_advRel_d3q27_incomp Subroutine

public subroutine mrt_advRel_d3q27_incomp(fieldProp, inState, outState, auxField, neigh, nElems, nSolve, level, layout, params, varSys, derVarPos)

Semi-optimized explicit implementation for incompressible model

This subroutine interface must match the abstract interface definition kernel in scheme/mus_scheme_type_module.f90 in order to be callable via compute function pointer.

First load all local values into temp array

Arguments

TypeIntentOptionalAttributesName
type(mus_field_prop_type), intent(in) :: fieldProp(:)

Array of field properties (fluid or species)

real(kind=rk), intent(in) :: inState(nElems*varSys%nScalars)

input pdf vector

real(kind=rk), intent(out) :: outState(nElems*varSys%nScalars)

output pdf vector

real(kind=rk), intent(inout) :: auxField(nElems*varSys%nAuxScalars)

Auxiliary field computed from pre-collision state Is updated with correct velocity field for multicomponent models

integer, intent(in) :: neigh(nElems*layout%fStencil%QQ)

connectivity vector

integer, intent(in) :: nElems

number of elements in state Array

integer, intent(in) :: nSolve

number of elements solved in kernel

integer, intent(in) :: level

current level

type(mus_scheme_layout_type), intent(in) :: layout

current layout

type(mus_param_type), intent(in) :: params

global parameters

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

variable system definition

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

position of derived quantities in varsys for all fields


Calls

proc~~mrt_advrel_d3q27_incomp~~CallsGraph proc~mrt_advrel_d3q27_incomp mrt_advRel_d3q27_incomp proc~mrt_d3q27 mrt_d3q27 proc~mrt_advrel_d3q27_incomp->proc~mrt_d3q27

Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: iElem
integer, private :: iDir
integer, private :: nScalars
real(kind=rk), private :: f000
real(kind=rk), private :: fN00
real(kind=rk), private :: f0N0
real(kind=rk), private :: f00N
real(kind=rk), private :: f100
real(kind=rk), private :: f010
real(kind=rk), private :: f001
real(kind=rk), private :: f0NN
real(kind=rk), private :: f0N1
real(kind=rk), private :: f01N
real(kind=rk), private :: f011
real(kind=rk), private :: fN0N
real(kind=rk), private :: f10N
real(kind=rk), private :: fN01
real(kind=rk), private :: f101
real(kind=rk), private :: fNN0
real(kind=rk), private :: fN10
real(kind=rk), private :: f1N0
real(kind=rk), private :: f110
real(kind=rk), private :: fNNN
real(kind=rk), private :: fNN1
real(kind=rk), private :: fN1N
real(kind=rk), private :: fN11
real(kind=rk), private :: f1NN
real(kind=rk), private :: f1N1
real(kind=rk), private :: f11N
real(kind=rk), private :: f111
real(kind=rk), private :: rho
real(kind=rk), private :: u_x
real(kind=rk), private :: u_y
real(kind=rk), private :: u_z
real(kind=rk), private :: usq
real(kind=rk), private :: sum_e1
real(kind=rk), private :: sum_e2
real(kind=rk), private :: sum_e3
real(kind=rk), private :: sum_ux1
real(kind=rk), private :: sum_ux2
real(kind=rk), private :: sum_ux3
real(kind=rk), private :: sum_uy1
real(kind=rk), private :: sum_uy2
real(kind=rk), private :: sum_uy3
real(kind=rk), private :: sum_uz1
real(kind=rk), private :: sum_uz2
real(kind=rk), private :: sum_uz3
real(kind=rk), private :: s_mrt(QQ)
real(kind=rk), private :: meq(QQ)
real(kind=rk), private :: mneq(QQ)
real(kind=rk), private :: mom(QQ)
real(kind=rk), private :: fneq(QQ)
integer, private :: dens_pos
integer, private :: vel_pos(3)
integer, private :: elemOff
real(kind=rk), private :: omegaBulk