def METHOD_NAME(self, open_dataset):
    """Replicates doctest to facilitate debugging."""
    tas = open_dataset("ERA5/daily_surface_cancities_1990-1993.nc").tas
    t90 = percentile_doy(
        tas.sel(time=slice("1990-01-01", "1991-12-31")), window=5, per=90
    )
    tg90p(tas=tas, tas_per=t90.isel(percentiles=0), freq="YS", bootstrap=True)
    tg90p(tas=tas, tas_per=t90, freq="YS", bootstrap=True)