Add time-varying forcing

Author: xylar

polaris/ocean/tasks/isomip_plus/__init__.py

@@ -73,7 +73,7 @@ def add_isomip_plus_tasks(component, mesh_type):
         )
 
         basic_expts = ['ocean0', 'ocean1', 'ocean2', 'ocean3', 'ocean4']
-        extra_expts = ['inception', 'wetting', 'drying']
+        extra_expts = ['wetting', 'drying']
         vert_coords = ['z-star', 'sigma']
 
         for experiment in basic_expts:
@@ -216,7 +216,7 @@ def _get_shared_steps(
     # ocean0 and ocean1 use the same topography
     shared_steps['ocean0'] = shared_steps['ocean1']
 
-    for experiment in ['inception', 'wetting', 'drying']:
+    for experiment in ['wetting', 'drying']:
         shared_steps[experiment] = dict(shared_steps['ocean1'])
         subdir = f'{resdir}/topo/scale/{experiment}'
         topo_scale = TopoScale(

polaris/ocean/tasks/isomip_plus/init/forcing.py

@@ -92,6 +92,7 @@ def run(self):
         Create restoring and other forcing files
         """
         self._compute_restoring()
+        self._write_time_varying_forcing()
 
     def _compute_restoring(self):
         config = self.config
@@ -182,3 +183,40 @@ def _compute_restoring(self):
             ds_forcing['tidalInputMask'] = xr.zeros_like(ds_mesh.xIsomipCell)
 
         write_netcdf(ds_forcing, 'init_mode_forcing_data.nc')
+
+    @staticmethod
+    def _write_time_varying_forcing():
+        """
+        Write time-varying land-ice forcing and update the initial condition
+        """
+
+        ds_topo = xr.open_dataset('topo.nc')
+
+        if 'Time' not in ds_topo.dims:
+            # no time-varying forcing needed
+            return
+
+        ds_out = xr.Dataset()
+        ds_out['xtime'] = ds_topo.xtime
+
+        ds_out['landIceDraftForcing'] = ds_topo.landIceDraft
+        ds_out.landIceDraftForcing.attrs['units'] = 'm'
+        ds_out.landIceDraftForcing.attrs['long_name'] = (
+            'The approximate elevation of the land ice-ocean interface'
+        )
+        ds_out['landIcePressureForcing'] = ds_topo.landIcePressure
+        ds_out.landIcePressureForcing.attrs['units'] = 'm'
+        ds_out.landIcePressureForcing.attrs['long_name'] = (
+            'Pressure from the weight of land ice at the ice-ocean interface'
+        )
+        ds_out['landIceFractionForcing'] = ds_topo.landIceFraction
+        ds_out.landIceFractionForcing.attrs['long_name'] = (
+            'The fraction of each cell covered by land ice'
+        )
+        ds_out['landIceFloatingFractionForcing'] = (
+            ds_topo.landIceFloatingFraction
+        )
+        ds_out.landIceFloatingFractionForcing.attrs['long_name'] = (
+            'The fraction of each cell covered by floating land ice'
+        )
+        write_netcdf(ds_out, 'land_ice_forcing.nc')

polaris/ocean/tasks/isomip_plus/isomip_plus_topo.cfg

@@ -27,17 +27,12 @@ expand_factor = 2.0
 # simple thickening and thinning experiments that involve scaling the Ocean1
 # landIcePressure and landIceDraft over time
 #
-# "inception" reference dates
-inception_dates = 0001-01-01_00:00:00, 0002-01-01_00:00:00, 0003-01-01_00:00:00
-# scaling at each date
-inception_scales = 0.0, 1.0, 1.0
-
 # "drying" reference dates
 drying_dates = 0001-01-01_00:00:00, 0002-01-01_00:00:00, 0003-01-01_00:00:00
 # scaling at each date
-drying_scales = 1.0, 2.0, 2.0
+drying_scales = 1.0, 1.05, 1.10
 
 # "wetting" reference dates
 wetting_dates = 0001-01-01_00:00:00, 0002-01-01_00:00:00, 0003-01-01_00:00:00
 # scaling at each date
-wetting_scales = 1.0, 0.0, 0.0
+wetting_scales = 1.0, 0.95, 0.9

polaris/ocean/tasks/isomip_plus/topo/scale.py

@@ -14,7 +14,7 @@ class TopoScale(Step):
 
     Attributes
     ----------
-    experiment : {'inception', 'drying', 'wetting'}
+    experiment : {'drying', 'wetting'}
         The scaling experiment
     """
 
@@ -36,10 +36,10 @@ def __init__(self, component, subdir, config, topo_remap, experiment):
         topo_remap : polaris.ocean.tasks.isomip_plus.topo_map.TopoRemap
             The step that produced the remapped topography to be scaled
 
-        experiment : {'inception', 'drying', 'wetting'}
+        experiment : {'drying', 'wetting'}
             The scaling experiment
         """
-        if experiment not in ['inception', 'drying', 'wetting']:
+        if experiment not in ['drying', 'wetting']:
             raise ValueError(f'Unexpected experiment {experiment}')
 
         super().__init__(component=component, name='topo_scale', subdir=subdir)