dynamic_viscosity Function

private pure function dynamic_viscosity(Z, T)

Arguments

TypeIntentOptionalAttributesName
real(kind=wp), intent(in) :: Z

altitude (m)

real(kind=wp), intent(in) :: T

temperature (K)

Return Value real(kind=wp)


Contents

Source Code


Source Code

pure real(wp) function dynamic_viscosity(Z,T)
    real(wp),intent(in) :: Z !! altitude (m)
    real(wp),intent(in) :: T !! temperature (K)
    if (Z > 86000.0_wp) error stop 'unable to compute dynamic viscosity above 86km'
    dynamic_viscosity = 1.458e-6_wp * T**(3.0_wp/2.0_wp) / (T + 110.4_wp)
end function dynamic_viscosity