central_moment Function

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

Calculating central moment. This follows equations 21 in cumulent paper ( Geier .et al 2015 )

Arguments

Type	Intent	Optional	Attributes		Name
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)

Contents

Variables

ii jj kk ii_rk

Variables

Type	Visibility	Attributes		Name		Initial
integer,	private		::	ii
integer,	private		::	jj
integer,	private		::	kk
real(kind=rk),	private,	parameter	::	ii_rk(-1:1)	=	[-1._rk, 0._rk, 1._rk]