getEquilibrium Function

private pure function getEquilibrium(density, velocity, layout) result(equil)

Calculate the equilibrium distribution function in all directions

The equilibrium distribution function is:\n \n where is the weight in each direction,\n is the macroscopic value of density,\n is the speed of sound,\n is the lattice unit velocity in each direction,\n is the macroscopic value of velocity.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: density
real(kind=rk), intent(in) :: velocity(3)
type(mus_scheme_layout_type), intent(in) :: layout

scheme layout

Return Value real(kind=rk)(layout%fStencil%QQ)

return value


Called by

proc~~getequilibrium~~CalledByGraph proc~getequilibrium getEquilibrium proc~fillmyghostsfromfiner_quadraticcompact3d fillMyGhostsFromFiner_quadraticCompact3d proc~fillmyghostsfromfiner_quadraticcompact3d->proc~getequilibrium

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: ucx
real(kind=rk), private :: usq
integer, private :: iDir