Introduce outlet glacier regions

Create region masks for each tidewater outlet glacier in domain

Author: alexolinhager

compass/landice/tests/uummannaq_disko/mesh.py

@@ -9,6 +9,7 @@
     build_mali_mesh,
     clean_up_after_interp,
     interp_gridded2mali,
+    make_region_masks,
 )
 from compass.model import make_graph_file
 from compass.step import Step
@@ -41,6 +42,8 @@ def __init__(self, test_case):
         self.mesh_filename = 'UummannaqDisko.nc'
         self.add_output_file(filename=self.mesh_filename)
         self.add_output_file(filename='graph.info')
+        self.add_output_file(
+            filename=f'{self.mesh_filename[:-3]}_regionMasks.nc')
         self.add_input_file(
             filename='greenland_1km_2024_01_29.epsg3413.icesheetonly.nc',
             target='greenland_1km_2024_01_29.epsg3413.icesheetonly.nc',
@@ -121,6 +124,13 @@ def run(self):
         make_graph_file(mesh_filename=self.mesh_filename,
                         graph_filename='graph.info')
 
+        mask_filename = f'{self.mesh_filename[:-3]}_regionMasks.nc'
+        make_region_masks(self, self.mesh_filename, mask_filename,
+                          self.cpus_per_task,
+                          tags=['UummannaqDisko'],
+                          component='landice',
+                          all_tags=False)
+
         # Do some final validation of the mesh
         ds = xr.open_dataset(self.mesh_filename)
         # Ensure basalHeatFlux is positive