chic Function

private function chic(g, k, l)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in) :: g
real(kind=rp), intent(in) :: k
real(kind=rp), intent(in) :: l

Return Value real(kind=rp)


Contents

Source Code


Source Code

  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