Merge pull request #730 from trhille/landice/fix_antarctica_mesh_gen

Fix `landice/antarctica/mesh_gen`. This reverts commit 8b327bd, which inadvertently broke the landice/antarctica/mesh_gen case and could lead to unexpected behavior elsewhere.
MPAS-Dev · Nov 13, 2023 · e5c0e82 · e5c0e82
2 parents 5327d22 + 6ff2656
commit e5c0e82
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/compass/landice/mesh.py b/compass/landice/mesh.py
@@ -529,10 +529,10 @@ def build_cell_width(self, section_name, gridded_dataset,
     geom_points, geom_edges = set_rectangular_geom_points_and_edges(*bnds)
 
     # Remove ice not connected to the ice sheet.
-    flood_mask = gridded_flood_fill(thk) == 0
-    thk[flood_mask] = 0.0
-    vx[flood_mask] = 0.0
-    vy[flood_mask] = 0.0
+    flood_mask = gridded_flood_fill(thk)
+    thk[flood_mask == 0] = 0.0
+    vx[flood_mask == 0] = 0.0
+    vy[flood_mask == 0] = 0.0
 
     # Calculate distance from each grid point to ice edge
     # and grounding line, for use in cell spacing functions.