Remap conserve with smoothing to base mesh

We were not doing this before because the base mesh might be periodic, but periodicity is now being removed.
E3SM-Project · Jul 15, 2024 · 3407555 · 3407555
1 parent ea5f99d
commit 3407555
Showing 1 changed file with 5 additions and 5 deletions.
diff --git a/polaris/ocean/tasks/isomip_plus/__init__.py b/polaris/ocean/tasks/isomip_plus/__init__.py
@@ -75,18 +75,18 @@ def _get_shared_steps(mesh_type, resolution, mesh_name, resdir, component,
         base_mesh.set_shared_config(config, link='isomip_plus_topo.cfg')
 
     subdir = f'{resdir}/topo/map_base'
-    # we remap the topography onto the base mesh without smoothing and with
-    # the bilinear method because the conserve method doesn't work with
-    # periodic meshes
+    # we remap the topography onto the base mesh with smoothing, which
+    # requires the conserve method.  This only works because we have already
+    # culled planar meshes to remove periodicity
     topo_map_base = TopoMap(component=component,
                             name='topo_map_base',
                             subdir=subdir,
                             config=config,
                             mesh_name=mesh_name,
                             mesh_step=base_mesh,
                             mesh_filename='base_mesh.nc',
-                            method='bilinear',
-                            smooth=False)
+                            method='conserve',
+                            smooth=True)
 
     subdir = f'{resdir}/topo/remap_base'
     topo_remap_base = TopoRemap(component=component,