central_moment_split Function

  1. mus_compute_cumulant_module.f90
  2. mus_compute_cumulant_module
  3. central_moment_split

public pure function central_moment_split(f, a, b, g, ux, uy, uz) result(kappa)

Calculating central moment by spliting among directions. This follows equations 43, 44, 45 in cumulent paper (Geier .et al 2015) We first do x direction for better performance.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: f(-1:1,-1:1,-1:1)

PDF
integer, intent(in) :: a

order of central moments
integer, intent(in) :: b

order of central moments
integer, intent(in) :: g

order of central moments
real(kind=rk), intent(in) :: ux
real(kind=rk), intent(in) :: uy
real(kind=rk), intent(in) :: uz

Return Value real(kind=rk)


Called by

proc~~central_moment_split~~CalledByGraph proc~central_moment_split central_moment_split proc~cumulant_d3q27 cumulant_d3q27 proc~cumulant_d3q27->proc~central_moment_split proc~cascaded_d3q27 cascaded_d3q27 proc~cascaded_d3q27->proc~central_moment_split

Contents

Variables


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: ka(-1:1,-1:1)
real(kind=rk), private :: kb(-1:1)
integer, private :: ii
integer, private :: jj
integer, private :: kk
real(kind=rk), private, parameter:: ii_rk(-1:1) =[-1._rk, 0._rk, 1._rk]