formatting

src/xsdba/_adjustment.py

@@ -22,7 +22,6 @@
 from .options import set_options
 from .processing import escore, jitter_under_thresh, reordering, standardize
 from .units import convert_units_to, units
-
 from .xclim_submodules.stats import _fitfunc_1d
 
 
@@ -716,7 +715,7 @@ def npdf_transform(ds: xr.Dataset, **kwargs) -> xr.Dataset:
     Notes
     -----
     If `n_escore` is negative, `escores` will be filled with NaNs.
-    
+
     """
     ref = ds.ref.rename(time_hist="time")
     hist = ds.hist.rename(time_hist="time")

src/xsdba/adjustment.py

@@ -20,11 +20,11 @@
 from xsdba.utils import uses_dask
 
 from ._adjustment import (
-    extremes_adjust, 
-    extremes_train,
     dqm_adjust,
     dqm_train,
     eqm_train,
+    extremes_adjust,
+    extremes_train,
     loci_adjust,
     loci_train,
     mbcn_adjust,
@@ -45,10 +45,8 @@
     pc_matrix,
     rand_rot_matrix,
 )
-
 from .xclim_submodules import stats
 
-
 __all__ = [
     "LOCI",
     "BaseAdjustment",

src/xsdba/units.py

@@ -222,7 +222,6 @@ def pint_multiply(
 }
 
 
-
 def ensure_absolute_temperature(units: str):
     """Convert temperature units to their absolute counterpart, assuming they represented a difference (delta).
 

tests/test_adjustment.py

@@ -8,10 +8,10 @@
 
 from xsdba import adjustment
 from xsdba.adjustment import (
-    ExtremeValues,
     LOCI,
     DetrendedQuantileMapping,
     EmpiricalQuantileMapping,
+    ExtremeValues,
     PrincipalComponents,
     QuantileDeltaMapping,
     Scaling,
@@ -765,14 +765,14 @@ def gen_testdata(c, s):
 
     @pytest.mark.slow
     def test_real_data(self, open_dataset):
-        dsim = open_dataset("sdba/CanESM2_1950-2100.nc")#.chunk()
-        dref = open_dataset("sdba/ahccd_1950-2013.nc")#.chunk()
+        dsim = open_dataset("sdba/CanESM2_1950-2100.nc")  # .chunk()
+        dref = open_dataset("sdba/ahccd_1950-2013.nc")  # .chunk()
         ref = dref.sel(time=slice("1950", "2009")).pr
         hist = dsim.sel(time=slice("1950", "2009")).pr
         # TODO: Do we want to include standard conversions in xsdba tests?
         # this is just convenient for now to keep those tests
         hist = pint_multiply(hist, "1e-03 m^3/kg")
-        hist = convert_units_to(hist,ref)
+        hist = convert_units_to(hist, ref)
 
         quantiles = np.linspace(0.01, 0.99, num=50)
 

tests/test_utils.py

@@ -15,7 +15,7 @@ def test_ecdf(timeseries, random):
     r = dist.rvs(10000, random_state=random)
     q = [0.01, 0.5, 0.99]
     x = xr.DataArray(dist.ppf(q), dims=("q",))
-    np.testing.assert_allclose(u.ecdf(timeseries(r, units="K" ), x), q, 3)
+    np.testing.assert_allclose(u.ecdf(timeseries(r, units="K"), x), q, 3)
 
     # With NaNs
     r[:2000] = np.nan
@@ -31,8 +31,8 @@ def test_map_cdf(timeseries, random):
     x_value = u.map_cdf(
         xr.Dataset(
             dict(
-                x=timeseries(xd.rvs(n, random_state=random), units = "kg / m^2 / s"),
-                y=timeseries(yd.rvs(n, random_state=random), units = "kg / m^2 / s"),
+                x=timeseries(xd.rvs(n, random_state=random), units="kg / m^2 / s"),
+                y=timeseries(yd.rvs(n, random_state=random), units="kg / m^2 / s"),
             )
         ),
         y_value=yd.ppf(q),