mus_advRel_kFluid_rBGK_vStd_lD3Q19 Subroutine

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

Advection relaxation routine for the D3Q19 model with BGK.

$$f_\alpha(x_i+e_{\alpha,i},t+1) = f_\alpha(x_i,t) - \omega(f_\alpha(x_i,t)-f^{eq}_{\alpha}(x_i,t))$$

The number of floating point operation in this routine is 160 roughly.

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.

Arguments

Type	Intent	Optional	Attributes		Name
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