Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rp), | intent(in) | :: | g | |||
real(kind=rp), | intent(in) | :: | k | |||
real(kind=rp), | intent(in) | :: | l |
FUNCTION chic(g,k,l)
IMPLICIT NONE
REAL(rp) :: chic
REAL(rp), INTENT(IN) :: g
REAL(rp), INTENT(IN) :: k
REAL(rp), INTENT(IN) :: l
REAL(rp) :: D
REAL(rp) :: xi
xi = 2.0_rp*C_PI/(3.0_rp*l*k*g**3)
D = (0.5_rp*(SQRT(4.0_rp + (C_PI/xi)**2) - C_PI/xi))**(1.0_rp/3.0_rp)
chic = (1.0_rp/D - D)/g
END FUNCTION chic