From 012685ce27fdf6211bf1a1f3e7ceb91bb45e3a66 Mon Sep 17 00:00:00 2001 From: Javier Vegas-Regidor Date: Thu, 10 Dec 2020 11:52:36 +0100 Subject: [PATCH] Fix invalid values in density --- earthdiagnostics/ocean/density.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/earthdiagnostics/ocean/density.py b/earthdiagnostics/ocean/density.py index d4d4aea4..69b070c4 100644 --- a/earthdiagnostics/ocean/density.py +++ b/earthdiagnostics/ocean/density.py @@ -134,7 +134,7 @@ class Density(Diagnostic): np.full(bigthetao.shape[1:], ref_pressure, dtype=np.float32) )) - sigma_values = np.stack(sigma_values) + sigma_values = np.ma.masked_invalid(np.stack(sigma_values)) sigma_cube = bigthetao.copy(sigma_values) sigma_cube.var_name = f'sigma{sigma}' sigma_cube.standard_name = 'sea_water_sigma_theta' @@ -144,7 +144,7 @@ class Density(Diagnostic): ) sigma_cube.units = 'kg m-3' temp = TempFile.get() - iris.save(sigma_cube, temp, zlib=True) + iris.save(sigma_cube, temp, zlib=True, fill_value=1e20) del sigma_cube del sigma_values self.sigma[sigma].set_local_file(temp) -- GitLab