Calculate_kr_VG_Base Module Function

module function Calculate_kr_VG_Base(alpha1, n1, m1, l, h) result(kr)

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(in) :: alpha1
real(kind=real64), intent(in) :: n1
real(kind=real64), intent(in) :: m1
real(kind=real64), intent(in) :: l
real(kind=real64), intent(in) :: h

Return Value real(kind=real64)


Called by

proc~~calculate_kr_vg_base~~CalledByGraph proc~calculate_kr_vg_base Calculate_kr_VG_Base interface~calculate_kr_vg_base Calculate_kr_VG_Base interface~calculate_kr_vg_base->proc~calculate_kr_vg_base proc~calculate_kr_base_impedance_vg Calculate_kr_Base_Impedance_VG proc~calculate_kr_base_impedance_vg->interface~calculate_kr_vg_base proc~calculate_kr_base_impedance_viscosity_vg Calculate_kr_Base_Impedance_Viscosity_VG proc~calculate_kr_base_impedance_viscosity_vg->interface~calculate_kr_vg_base proc~calculate_kr_base_vg Calculate_kr_Base_VG proc~calculate_kr_base_vg->interface~calculate_kr_vg_base proc~calculate_kr_base_viscosity_vg Calculate_kr_Base_Viscosity_VG proc~calculate_kr_base_viscosity_vg->interface~calculate_kr_vg_base interface~calculate_kr_base_impedance_vg Type_HCF_Base_Impedance_VG%Calculate_kr_Base_Impedance_VG interface~calculate_kr_base_impedance_vg->proc~calculate_kr_base_impedance_vg interface~calculate_kr_base_impedance_viscosity_vg Type_HCF_Base_Impedance_Viscosity_VG%Calculate_kr_Base_Impedance_Viscosity_VG interface~calculate_kr_base_impedance_viscosity_vg->proc~calculate_kr_base_impedance_viscosity_vg interface~calculate_kr_base_vg Type_HCF_Base_VG%Calculate_kr_Base_VG interface~calculate_kr_base_vg->proc~calculate_kr_base_vg interface~calculate_kr_base_viscosity_vg Type_HCF_Base_Viscosity_VG%Calculate_kr_Base_Viscosity_VG interface~calculate_kr_base_viscosity_vg->proc~calculate_kr_base_viscosity_vg proc~calculate_kflh_base_impedance_vg Calculate_Kflh_Base_Impedance_VG proc~calculate_kflh_base_impedance_vg->interface~calculate_kr_base_impedance_vg proc~calculate_kflh_base_impedance_viscosity_vg Calculate_Kflh_Base_Impedance_Viscosity_VG proc~calculate_kflh_base_impedance_viscosity_vg->interface~calculate_kr_base_impedance_viscosity_vg proc~calculate_kflh_base_vg Calculate_Kflh_Base_VG proc~calculate_kflh_base_vg->interface~calculate_kr_base_vg proc~calculate_kflh_base_viscosity_vg Calculate_Kflh_Base_Viscosity_VG proc~calculate_kflh_base_viscosity_vg->interface~calculate_kr_base_viscosity_vg interface~calculate_kflh_base_impedance_vg Type_HCF_Base_Impedance_VG%Calculate_Kflh_Base_Impedance_VG interface~calculate_kflh_base_impedance_vg->proc~calculate_kflh_base_impedance_vg interface~calculate_kflh_base_impedance_viscosity_vg Type_HCF_Base_Impedance_Viscosity_VG%Calculate_Kflh_Base_Impedance_Viscosity_VG interface~calculate_kflh_base_impedance_viscosity_vg->proc~calculate_kflh_base_impedance_viscosity_vg interface~calculate_kflh_base_vg Type_HCF_Base_VG%Calculate_Kflh_Base_VG interface~calculate_kflh_base_vg->proc~calculate_kflh_base_vg interface~calculate_kflh_base_viscosity_vg Type_HCF_Base_Viscosity_VG%Calculate_Kflh_Base_Viscosity_VG interface~calculate_kflh_base_viscosity_vg->proc~calculate_kflh_base_viscosity_vg proc~update_kflh_base_impedance_vg Update_Kflh_Base_Impedance_VG proc~update_kflh_base_impedance_vg->interface~calculate_kflh_base_impedance_vg proc~update_kflh_base_impedance_viscosity_vg Update_Kflh_Base_Impedance_Viscosity_VG proc~update_kflh_base_impedance_viscosity_vg->interface~calculate_kflh_base_impedance_viscosity_vg proc~update_kflh_base_vg Update_Kflh_Base_VG proc~update_kflh_base_vg->interface~calculate_kflh_base_vg proc~update_kflh_base_viscosity_vg Update_Kflh_Base_Viscosity_VG proc~update_kflh_base_viscosity_vg->interface~calculate_kflh_base_viscosity_vg interface~update_kflh_base_impedance_vg Type_HCF_Base_Impedance_VG%Update_Kflh_Base_Impedance_VG interface~update_kflh_base_impedance_vg->proc~update_kflh_base_impedance_vg interface~update_kflh_base_impedance_viscosity_vg Type_HCF_Base_Impedance_Viscosity_VG%Update_Kflh_Base_Impedance_Viscosity_VG interface~update_kflh_base_impedance_viscosity_vg->proc~update_kflh_base_impedance_viscosity_vg interface~update_kflh_base_vg Type_HCF_Base_VG%Update_Kflh_Base_VG interface~update_kflh_base_vg->proc~update_kflh_base_vg interface~update_kflh_base_viscosity_vg Type_HCF_Base_Viscosity_VG%Update_Kflh_Base_Viscosity_VG interface~update_kflh_base_viscosity_vg->proc~update_kflh_base_viscosity_vg

Source Code

    module function Calculate_kr_VG_Base(alpha1, n1, m1, l, h) result(kr)
        !$omp declare simd uniform(alpha1, n1, m1, l, h)
        implicit none
        real(real64), intent(in) :: alpha1
        real(real64), intent(in) :: n1
        real(real64), intent(in) :: m1
        real(real64), intent(in) :: l
        real(real64), intent(in) :: h
        real(real64) :: kr
        real(real64) :: Sw

        if (h < 0.0d0) then
            Sw = (1.0d0 + (-alpha1 * h)**n1)**(-m1)
        else
            Sw = 1.0d0
        end if

        kr = Sw**l * (1.0d0 - (1.0d0 - Sw**(1.0d0 / m1))**m1)**2.0d0

    end function Calculate_kr_VG_Base