set_pdfAcoustic Function

public function set_pdfAcoustic(layout, omega, rho0, mom, incompressible) result(pdf)

Calculate the distribution function in all directions by using the fEq + fNeq

The input vector must look like 2D:

 mom = [ rho, ux, uy, Sxx, Syy, Sxy ]
 mom = [ rho, ux, uy, uz, Sxx, Syy, Szz, Sxy, Syz, Sxz ]

Arguments

TypeIntentOptionalAttributesName
type(mus_scheme_layout_type), intent(in) :: layout
real(kind=rk), intent(in) :: omega
real(kind=rk), intent(in) :: rho0
real(kind=rk), intent(in) :: mom(layout%fStencil%QQ)
logical, intent(in), optional :: incompressible

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


Calls

proc~~set_pdfacoustic~~CallsGraph proc~set_pdfacoustic set_pdfAcoustic proc~getneq_acoustic getNEq_acoustic proc~set_pdfacoustic->proc~getneq_acoustic proc~getequilibriumincomp getEquilibriumIncomp proc~set_pdfacoustic->proc~getequilibriumincomp interface~getequilibrium getEquilibrium proc~set_pdfacoustic->interface~getequilibrium proc~convprepost convPrePost proc~getneq_acoustic->proc~convprepost proc~getequilibrium_forpdfsubset getEquilibrium_forPdfSubset interface~getequilibrium->proc~getequilibrium_forpdfsubset proc~getequilibrium_forelemfromstate getEquilibrium_forElemfromState interface~getequilibrium->proc~getequilibrium_forelemfromstate proc~geteqbydensvel getEqByDensVel interface~getequilibrium->proc~geteqbydensvel proc~getvelocity_forpdfsubset getVelocity_forPdfSubset proc~getequilibrium_forpdfsubset->proc~getvelocity_forpdfsubset proc~getdensity_forpdfsubset getDensity_forPdfSubset proc~getequilibrium_forpdfsubset->proc~getdensity_forpdfsubset interface~getvelocity getVelocity proc~getequilibrium_forelemfromstate->interface~getvelocity interface~getdensity getDensity proc~getequilibrium_forelemfromstate->interface~getdensity interface~getvelocity->proc~getvelocity_forpdfsubset proc~getvelocity_forelemfromstate_noforce getVelocity_forElemFromState_noForce interface~getvelocity->proc~getvelocity_forelemfromstate_noforce interface~getdensity->proc~getdensity_forpdfsubset proc~getdensity_forelemfromstate getDensity_forElemFromState interface~getdensity->proc~getdensity_forelemfromstate

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: fEq(layout%fStencil%QQ)
real(kind=rk), private :: fnEq(layout%fStencil%QQ)
real(kind=rk), private :: moms(layout%fStencil%QQ)
logical, private :: incomp =.false.