deriveEquilIncomp_FromMacro_d3q19 Subroutine

public subroutine deriveEquilIncomp_FromMacro_d3q19(density, velocity, iField, nElems, varSys, layout, res)

This routine computes equilbrium from density and velocity

This subroutine's interface must match the abstract interface definition derive_FromMacro in derived/mus_derVarPos_module.f90 in order to be callable via equilFromMacro function pointer.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: density(:)

Array of density. Single species: dens_1, dens_2 .. dens_n multi-species: dens_1_sp1, dens_1_sp2, dens_2_sp1, dens_2_sp2 ... dens_n_sp1, dens_n_sp2

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

Array of velocity. Size: dimension 1: n*nFields. dimension 2: 3 (nComp) 1st dimension arrangement for multi-species is same as density

integer, intent(in) :: iField

Current field

integer, intent(in) :: nElems

number of elements

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

variable system which is required to access fieldProp information via variable method data c_ptr

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

scheme layout contains stencil definition and lattice weights

real(kind=rk), intent(out) :: res(:)

Output of this routine Dimension: n*nComponents of res


Contents


Variables

TypeVisibilityAttributesNameInitial
integer, private :: iElem
real(kind=rk), private :: fEq(19)
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 :: usqn
real(kind=rk), private :: usqn_o2
real(kind=rk), private :: coeff_1
real(kind=rk), private :: coeff_2
real(kind=rk), private :: ui1
real(kind=rk), private :: ui3
real(kind=rk), private :: ui10
real(kind=rk), private :: ui11
real(kind=rk), private :: ui12
real(kind=rk), private :: ui13
real(kind=rk), private :: fac_1
real(kind=rk), private :: fac_2
real(kind=rk), private :: fac_3
real(kind=rk), private :: fac_4
real(kind=rk), private :: fac_9
real(kind=rk), private :: fac_10
real(kind=rk), private :: fac_11
real(kind=rk), private :: fac_12
real(kind=rk), private :: fac_13
real(kind=rk), private :: sum1_1
real(kind=rk), private :: sum1_2
real(kind=rk), private :: sum2_1
real(kind=rk), private :: sum2_2
real(kind=rk), private :: sum3_1
real(kind=rk), private :: sum3_2
real(kind=rk), private :: sum4_1
real(kind=rk), private :: sum4_2
real(kind=rk), private :: sum9_1
real(kind=rk), private :: sum9_2
real(kind=rk), private :: sum10_1
real(kind=rk), private :: sum10_2
real(kind=rk), private :: sum11_1
real(kind=rk), private :: sum11_2
real(kind=rk), private :: sum12_1
real(kind=rk), private :: sum12_2
real(kind=rk), private :: sum13_1
real(kind=rk), private :: sum13_2