Skip to content

Commit

Permalink
only check if no nan values within the model domain. This worked prev…
Browse files Browse the repository at this point in the history
…iously because we only have been simulating inland subbasins
  • Loading branch information
jensdebruijn committed Sep 19, 2024
1 parent 132e60c commit f42fb53
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion geb/HRUs.py
Original file line number Diff line number Diff line change
Expand Up @@ -443,7 +443,7 @@ def spei_uncompressed(self):
if not hasattr(self, "spei_ds"):
self.spei_ds = self.load_forcing_ds("spei")
spei = self.load_forcing(self.spei_ds, self.model.current_time, compress=False)
assert not np.isnan(spei).any()
assert not np.isnan(spei[~self.mask]).any() # no nan values in non-masked cells
return spei

@property
Expand Down

0 comments on commit f42fb53

Please sign in to comment.