From 1a964f717f12068e4658d910f9f887748857493a Mon Sep 17 00:00:00 2001 From: sloosvel Date: Fri, 17 Jan 2020 10:04:54 +0100 Subject: [PATCH] Do not save basin if all values are 0 --- earthdiagnostics/ocean/siasiesiv.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/earthdiagnostics/ocean/siasiesiv.py b/earthdiagnostics/ocean/siasiesiv.py index 24cf3c0e..dd13da4e 100644 --- a/earthdiagnostics/ocean/siasiesiv.py +++ b/earthdiagnostics/ocean/siasiesiv.py @@ -210,7 +210,8 @@ class Siasiesiv(Diagnostic): else: var_res.units = 'm^2' for i, basin in enumerate(self.masks): - var_region[i, ...] = netCDF4.stringtoarr(str(basin), 50) - var_res[..., i] = res[i, ...] + if not np.all(res[i, ...]==0): + var_region[i, ...] = netCDF4.stringtoarr(str(basin), 50) + var_res[..., i] = res[i, ...] handler_temp.close() self.generated[var].set_local_file(temp, diagnostic=self) -- GitLab