def snowmelt_rain_on_snow_heavily_forested(precipitation,
                                           temperatures,
                                           temperature_cutoff=32.0,
                                           rain_melt_coeff=0.007):
    """Calculate the amount of snowmelt rain-on-snow situations in
    heavily forested areas using a generalized equation for rain-on-snow
    situations in heavily forested areas (the mean canopy cover is greater
    than 80%). This snowmelt calculation is from the family of energy budget
    solutions.

    :param precipitation: Precipitation rates, in inches/day
    :type precipitation: numpy.ndarray
    :param temperatures: Temperatures, in degrees Fahrenheit
    :type temperatures: numpy.ndarray
    :param temperature_cutoff: Temperature when melt begins,
                               in degrees Fahrenheit
    :type temperature_cutoff: float
    :param rain_melt_coeff: Snowmelt coefficient when raining,
                             1/degrees Fahrenheit
    :type rain_melt_coeff: float
    :return snowmelt: Snowmelt values, in inches per day
    :rtype snowmelt: numpy.ndarray

    .. note::

    Equation:

    M = (0.074 + 0.007 * P_r) * (T_a - 32) + 0.05

    where

    M           snowmelt, inches/day
    P_r         rate of precipitation, inches/day
    T_a         temperature of saturated air, at 3 meters (10 ft) level,
                degrees Fahrenheit

    Reference:

    Engineering and Design - Runoff from Snowmelt
    U.S. Army Corps of Engineers
    Engineering Manual 1110-2-1406
    Chapter 5-3. Generalized Equations, Rain-on-Snow Situations, Equation 5-20
    https://www.wcc.nrcs.usda.gov/ftpref/wntsc/H&H/snow/COEemSnowmeltRunoff.pdf
    """
    return (
        (0.074 + rain_melt_coeff * precipitation)
        * (temperatures - temperature_cutoff) + 0.05
    )