getEquilibrium_forPdfSubset Function

private pure function getEquilibrium_forPdfSubset(subset, layout, varPos) result(equil)

Calculate the equilibrium distribution function in all directions

The equilibrim distribution function is:\n $$f^{eq}_i = w_i \rho ( 1 + \frac{\vec c_i \cdot \vec u}{c^2_s} + \frac{ {(\vec c_i \cdot \vec u)}^2}{2c^4_s} - \frac{\vec u \cdot \vec u}{2c^2_s})$$ \n where $w_i$ is the weight in each direction,\n $\rho$ is the macroscopic value of density,\n $c_s$ is the speed of sound,\n $\vec c_i$ is the lattice unit velocity in each direction,\n $\vec u$ is the macroscopic value of velocity.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in)			::	subset(:)
type(mus_scheme_layout_type),	intent(in)			::	layout
integer,	intent(in)			::	varPos(:)

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