Merge pull request #201 from xylar/isomip-plus-cull-periodic-base-mesh

Cull ISOMIP+ planar base mesh
E3SM-Project · Jul 15, 2024 · b97100b · b97100b
2 parents a87a3a4 + 3407555
commit b97100b
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 6 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,

diff --git a/polaris/ocean/tasks/isomip_plus/mesh/planar.py b/polaris/ocean/tasks/isomip_plus/mesh/planar.py
@@ -1,6 +1,7 @@
 import numpy as np
 import pyproj
 from mpas_tools.io import write_netcdf
+from mpas_tools.mesh.conversion import cull
 from mpas_tools.planar_hex import make_planar_hex_mesh
 from mpas_tools.translate import translate
 
@@ -79,4 +80,6 @@ def run(self):
             ds_mesh[f'lon{suffix}'] = (f'n{suffix}', np.deg2rad(lon))
 
         add_isomip_plus_xy(ds_mesh)
-        write_netcdf(ds_mesh, 'base_mesh.nc')
+        # cull the mesh to remove periodicity
+        ds_nonperiodic = cull(ds_mesh)
+        write_netcdf(ds_nonperiodic, 'base_mesh.nc')