From f38ec10fce4da3a5c97f8f78e461ff7519c6078b Mon Sep 17 00:00:00 2001 From: Carolyn Begeman Date: Wed, 25 Oct 2023 12:56:28 -0500 Subject: [PATCH] fixup thin film culling --- compass/ocean/tests/isomip_plus/geom.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/compass/ocean/tests/isomip_plus/geom.py b/compass/ocean/tests/isomip_plus/geom.py index 551f70b3d5..3c52630511 100755 --- a/compass/ocean/tests/isomip_plus/geom.py +++ b/compass/ocean/tests/isomip_plus/geom.py @@ -157,7 +157,7 @@ def _get_geom_fields(ds_geom, ds_mesh, thin_film_present): y_cell = ds_mesh.yIsomipCell.values if thin_film_present: - ocean_fraction = - ds_geom['landFraction'] + 1.0 + ocean_fraction = xarray.where(ds_geom['Z_bed'] > 0., 0., 1.) else: ocean_fraction = (ds_geom['landIceFloatingFraction'] + ds_geom['openOceanFraction'])