trt_advRel_d3q27 Subroutine

  1. mus_compute_d3q27_module.f90
  2. mus_d3q27_module
  3. trt_advRel_d3q27

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

No comment yet!

TODO add comment

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

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

Contents

Variables


Variables

TypeVisibilityAttributesNameInitial
integer, private :: iElem
integer, private :: nScalars
real(kind=rk), private :: f(-1:1,-1:1,-1:1)
real(kind=rk), private :: u
real(kind=rk), private :: v
real(kind=rk), private :: w
real(kind=rk), private :: u2
real(kind=rk), private :: v2
real(kind=rk), private :: w2
real(kind=rk), private :: rho
real(kind=rk), private :: inv_rho
real(kind=rk), private :: wP
real(kind=rk), private :: wN
real(kind=rk), private :: p_part
real(kind=rk), private :: n_part
real(kind=rk), private :: X0
real(kind=rk), private :: X1
real(kind=rk), private :: XN
real(kind=rk), private :: Y0
real(kind=rk), private :: Y1
real(kind=rk), private :: YN
real(kind=rk), private :: Z0
real(kind=rk), private :: Z1
real(kind=rk), private :: ZN
integer, private :: dens_pos
integer, private :: vel_pos(3)
integer, private :: elemOff