ptfkit.rawls1982

def calc_full_wrc_rawls1982(
    *,
    sand: float | ArrayLike,
    organic_matter: float | ArrayLike,
    bulk_density: float | ArrayLike,
    theta_33: float | ArrayLike,
    theta_1500: float | ArrayLike,
    out: Rawls1982PTFResult[NDArray[floating]] | None = None,
) -> Rawls1982PTFResult[floating] | Rawls1982PTFResult[NDArray[floating]]:
    r"""Estimate a twelve-point water-retention curve using theta_33 and theta_1500.

    Arguments:
        sand: Sand content. (%)
        organic_matter: Organic matter content. (%)
        bulk_density: Bulk density. (g/cm^3)
        theta_33: Measured or estimated volumetric water content at -33 kPa. (cm^3/cm^3)
        theta_1500: Measured or estimated volumetric water content at -1500 kPa. (cm^3/cm^3)
        out: Optional output arrays for in-place calculation.

    Returns:
        Rawls1982PTFResult: Results grouped by result attributes.

    Models:
        $h(\theta)$: Twelve point water-retention curve

    Notes:
        Prediction target: Volumetric water contents from -4 to -1500 kPa.
        The source offers three regression levels; this function uses measured theta_33 and
            theta_1500 for the intermediate points.

    Warning:
        The source value 0.8888 is retained literally for the theta_7 intercept.

    """
    values = _call(
        _calc_ptf_rawls1982_full_wrc,
        sand,
        organic_matter,
        bulk_density,
        theta_33,
        theta_1500,
        out=out,
    )

    return Rawls1982PTFResult(*values)

def calc_theta_1500_rawls1982(
    *,
    clay: float | ArrayLike,
    organic_matter: float | ArrayLike,
    out: NDArray[floating] | None = None,
) -> floating | NDArray[floating]:
    r"""Estimate volumetric water content at -1500 kPa.

    Arguments:
        clay: Clay content. (%)
        organic_matter: Organic matter content. (%)
        out: Optional output arrays for in-place calculation.

    Returns:
        theta_1500: Volumetric soil water content at -1500 kPa. (cm^3/cm^3)

    Models:
        $h(\theta)$: Point estimate at -1500 kPa

    Notes:
        Prediction target: Volumetric soil water content at -1500 kPa.

    """
    return _call(
        _calc_ptf_rawls1982_theta_1500,
        clay,
        organic_matter,
        out=out,
    )

def calc_theta_33_rawls1982(
    *,
    sand: float | ArrayLike,
    organic_matter: float | ArrayLike,
    theta_1500: float | ArrayLike,
    out: NDArray[floating] | None = None,
) -> floating | NDArray[floating]:
    r"""Estimate volumetric water content at -33 kPa using measured or estimated theta_1500.

    Arguments:
        sand: Sand content. (%)
        organic_matter: Organic matter content. (%)
        theta_1500: Measured or estimated volumetric water content at -1500 kPa. (cm^3/cm^3)
        out: Optional output arrays for in-place calculation.

    Returns:
        theta_33: Volumetric soil water content at -33 kPa. (cm^3/cm^3)

    Models:
        $h(\theta)$: Point estimate at -33 kPa

    Notes:
        Prediction target: Volumetric soil water content at -33 kPa.

    """
    return _call(
        _calc_ptf_rawls1982_theta_33,
        sand,
        organic_matter,
        theta_1500,
        out=out,
    )