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.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | f(-1:1,-1:1,-1:1) |
|
||
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 |