Add remap topography step to tides mesh step

Author: sbrus89

compass/ocean/tests/tides/mesh/__init__.py

@@ -1,5 +1,6 @@
 from compass.mesh.spherical import IcosahedralMeshStep
 from compass.ocean.mesh.cull import CullMeshStep
+from compass.ocean.mesh.remap_topography import RemapTopography
 from compass.ocean.tests.tides.configure import configure_tides
 from compass.ocean.tests.tides.dem import CreatePixelFile
 from compass.ocean.tests.tides.mesh.vr45to5 import VRTidesMesh
@@ -52,9 +53,14 @@ def __init__(self, test_group, mesh_name):
 
         self.add_step(base_mesh_step)
 
+        remap_step = RemapTopography(test_case=self,
+                                     base_mesh_step=base_mesh_step,
+                                     mesh_name=mesh_name)
+        self.add_step(remap_step)
+
         self.add_step(CullMeshStep(
             test_case=self, base_mesh_step=base_mesh_step,
-            with_ice_shelf_cavities=True))
+            with_ice_shelf_cavities=True, remap_topography=remap_step))
 
     def configure(self):
         """

compass/ocean/tests/tides/mesh/icos7/icos7.cfg

@@ -5,7 +5,7 @@
 # number of cores to use
 cull_mesh_cpus_per_task = 18
 # minimum of cores, below which the step fails
-cull_mesh_min_cpus_per_task = 1 
+cull_mesh_min_cpus_per_task = 1
 # maximum memory usage allowed (in MB)
 cull_mesh_max_memory = 1000
 
@@ -36,11 +36,11 @@ min_depth = 1
 # number of cores to use
 init_ntasks = 36
 # minimum of cores, below which the step fails
-init_min_tasks = 8 
+init_min_tasks = 8
 # maximum memory usage allowed (in MB)
 init_max_memory = 1000
 # number of threads
-init_threads = 1 
+init_threads = 1
 
 ## config options related to the forward steps
 # number of cores to use
@@ -52,5 +52,40 @@ forward_max_memory = 1000
 # number of threads
 forward_threads = 1
 
-# TPXO version for validation 
+# TPXO version for validation
 tpxo_version = TPXO9
+
+
+# config options related to remapping topography to an MPAS-Ocean mesh
+[remap_topography]
+
+# the name of the topography file in the bathymetry database
+topo_filename = BedMachineAntarctica_v3_and_GEBCO_2023_0.0125_degree_20240828.nc
+
+# variable names in topo_filename
+lon_var = lon
+lat_var = lat
+bathymetry_var = bathymetry
+ice_thickness_var = thickness
+ice_frac_var = ice_mask
+grounded_ice_frac_var = grounded_mask
+ocean_frac_var = ocean_mask
+bathy_frac_var = bathymetry_mask
+
+# the description to include in metadata
+description = Bathymetry is from GEBCO 2023, combined with BedMachine
+              Antarctica v3 around Antarctica.
+
+# the target and minimum number of MPI tasks to use in remapping
+ntasks = 4096
+min_tasks = 360
+
+# remapping method {'bilinear', 'neareststod', 'conserve'}
+method = conserve
+
+# threshold of what fraction of an MPAS cell must contain ocean in order to
+# perform renormalization of elevation variables
+renorm_threshold = 0.01
+
+# the density of land ice from MALI (kg/m^3)
+ice_density = 910.0

compass/ocean/tests/tides/mesh/vr45to5/vr45to5.cfg

@@ -60,3 +60,37 @@ forward_threads = 1
 
 # TPXO version for validation
 tpxo_version = TPXO9
+
+# config options related to remapping topography to an MPAS-Ocean mesh
+[remap_topography]
+
+# the name of the topography file in the bathymetry database
+topo_filename = BedMachineAntarctica_v3_and_GEBCO_2023_0.0125_degree_20240828.nc
+
+# variable names in topo_filename
+lon_var = lon
+lat_var = lat
+bathymetry_var = bathymetry
+ice_thickness_var = thickness
+ice_frac_var = ice_mask
+grounded_ice_frac_var = grounded_mask
+ocean_frac_var = ocean_mask
+bathy_frac_var = bathymetry_mask
+
+# the description to include in metadata
+description = Bathymetry is from GEBCO 2023, combined with BedMachine
+              Antarctica v3 around Antarctica.
+
+# the target and minimum number of MPI tasks to use in remapping
+ntasks = 4096
+min_tasks = 360
+
+# remapping method {'bilinear', 'neareststod', 'conserve'}
+method = conserve
+
+# threshold of what fraction of an MPAS cell must contain ocean in order to
+# perform renormalization of elevation variables
+renorm_threshold = 0.01
+
+# the density of land ice from MALI (kg/m^3)
+ice_density = 910.0