Merge pull request #887 from trhille/fix_landice_cull_distance_bug

Fix bug setting cell spacing outside cull distance

Add a necessary absolute value when setting large cell spacing outside of cull distance. Otherwise, this causes problems when using a negative cull distance to indicate culling to a geojson file.

Author: matthewhoffman

compass/landice/mesh.py

@@ -372,7 +372,7 @@ def set_cell_width(self, section_name, thk, bed=None, vx=None, vy=None,
     # that distance.
     if dist_to_edge is not None:
         mask = np.logical_and(
-            thk == 0.0, dist_to_edge > (3. * cull_distance))
+            thk == 0.0, dist_to_edge > np.abs(3. * cull_distance))
         logger.info('Setting cell_width in outer regions to max_spac '
                     f'for {mask.sum()} cells')
         cell_width[mask] = max_spac