deriveEquilMSGas Subroutine

private recursive subroutine deriveEquilMSGas(fun, varSys, elempos, time, tree, nElems, nDofs, res)

Calculate the equlibrium of a given element number with the given input state vector.

The equilibrium distribution function is:\n $$f^{eq}_rest = \rho w_i( (3-2 mRatio_iField ) + 3(cx . u) + 9/2 (cx . u )^2 - 3/2 u^2 )$$ $$f^{eq}_i = \rho w_i( mRatio_iField + 3(cx . u) + 9/2 (cx . u )^2 - 3/2 u^2 )$$ 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.

The interface has to comply to the abstract interface tem_varSys_proc_element.

Arguments

Type	Intent	Optional	Attributes		Name
class(tem_varSys_op_type),	intent(in)			::	fun	Description of the method to obtain the variables, here some preset values might be stored, like the space time function to use or the required variables.
type(tem_varSys_type),	intent(in)			::	varSys	The variable system to obtain the variable from.
integer,	intent(in)			::	elempos(:)	Position of the TreeID of the element to get the variable for in the global treeID list.
type(tem_time_type),	intent(in)			::	time	Point in time at which to evaluate the variable.
type(treelmesh_type),	intent(in)			::	tree	global treelm mesh info
integer,	intent(in)			::	nElems	Number of values to obtain for this variable (vectorized access).
integer,	intent(in)			::	nDofs	Number of degrees of freedom within an element.
real(kind=rk),	intent(out)			::	res(:)	Resulting values for the requested variable. Linearized array dimension: (n requested entries) x (nComponents of this variable) x (nDegrees of freedom) Access: (iElem-1)fun%nComponentsnDofs + (iDof-1)*fun%nComponents + iComp