diff --git a/mlair/run_modules/post_processing.py b/mlair/run_modules/post_processing.py index 304f004c7e4c904cc75f62f294920562b87258a7..a351dd92e70166dc5a4917a5bf7c3d5ef2b9686f 100644 --- a/mlair/run_modules/post_processing.py +++ b/mlair/run_modules/post_processing.py @@ -969,10 +969,16 @@ class PostProcessing(RunEnvironment): wind_sectors = self.data_store.get("wind_sectors", "general") h_sector_skill_scores = [] for sec in wind_sectors: - h_sector_skill_scores.append(statistics.skill_score_based_on_mse( - ds.where(self.upstream_wind_sector.squeeze() == sec), - obs_name=self.observation_indicator, pred_name=self.model_display_name, - ref_name=ref_name).assign_coords({"sector": sec}) + h_sector_skill_scores.append( + # statistics.SkillScores(None).general_skill_score(ds.where(self.upstream_wind_sector.squeeze() == sec), + # forecast_name=self.model_display_name, + # reference_name=ref_name, + # observation_name=self.observation_indicator) + statistics.skill_score_based_on_mse( + ds.where(self.upstream_wind_sector.squeeze() == sec).dropna(dim=self.index_dim), + obs_name=self.observation_indicator, pred_name=self.model_display_name, + ref_name=ref_name).assign_coords({"sector": sec} + ) ) sector_skill_scores = xr.concat(h_sector_skill_scores, dim="sector") sector_skill_scores = sector_skill_scores.assign_attrs({f"reference_model": ref_name})