Minor fixes

Author: DamienIrving

unseen/eva.py

@@ -1176,7 +1176,7 @@ def gev_return_curve(
         np.isfinite(boot_event_return_periods)
     ]
     event_return_period_lower_ci = np.quantile(boot_event_return_periods, q)
-    event_return_period_upper_ci = np.quantile(boot_event_return_periods, q - 1)
+    event_return_period_upper_ci = np.quantile(boot_event_return_periods, 1 - q)
     event_data = (
         event_return_period,
         event_return_period_lower_ci,

unseen/fileio.py

@@ -253,6 +253,7 @@ def open_dataset(
     if no_leap_days:
         ds = ds.sel(time=~((ds[time_dim].dt.month == 2) & (ds[time_dim].dt.day == 29)))
     if rolling_sum_window:
+        ds = ds.compute()
         ds = ds.rolling({time_dim: rolling_sum_window}).sum(dim=time_dim)
     if time_freq:
         assert time_agg, "Provide a time_agg"

unseen/process_utils.py

@@ -60,7 +60,7 @@ def count_months(event_months):
     return counts
 
 
-def plot_event_seasonality(ds, ax=None, core_dim=None, outfile=None):
+def plot_event_seasonality(ds, ax=None, core_dim="time", outfile=None):
     """Plot event seasonality
 
     Parameters