def extr_hotdays_calc(data, thr_p95):
    """
    Calculate number of extreme hotdays.

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

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