def hotdays_calc(data, thr_p75):
    """
    Calculate number of hotdays.

    Return days with mean temperature above the 75th percentile
    of climatology.

    Parameters
    ----------
    data: array
        1D-array of temperature input timeseries
    thr_p75: float
        75th percentile daily mean value from climatology
    """
    hotdays = ((data > thr_p75)).sum()
    return hotdays