Source code for compass.ocean.tests.isomip_plus.ssh_adjustment
-from compass.ocean.iceshelf import adjust_ssh
-from compass.ocean.tests.isomip_plus.forward import get_time_steps
+import time
+
+from compass.ocean.iceshelf import adjust_ssh
from compass.step import Step
@@ -140,6 +141,8 @@ Source code for compass.ocean.tests.isomip_plus.ssh_adjustment
ntasks=None, min_tasks=None,
openmp_threads=None)
+ self.resolution = resolution
+
# generate the namelist, replacing a few default options
# start with the same namelist settings as the forward run
self.add_namelist_file('compass.ocean.tests.isomip_plus',
@@ -154,8 +157,6 @@ Source code for compass.ocean.tests.isomip_plus.ssh_adjustment
# we don't want the global stats AM for this run
options = dict()
- if not thin_film_present:
- options = get_time_steps(resolution)
options['config_AM_globalStats_enable'] = '.false.'
self.add_namelist_options(options)
@@ -195,6 +196,21 @@ Source code for compass.ocean.tests.isomip_plus.ssh_adjustment
Run this step of the test case
"""
config = self.config
+ resolution = self.resolution
+
+ dt_per_km = config.getfloat('isomip_plus', 'dt_per_km')
+ dt_btr_per_km = config.getfloat('isomip_plus', 'dt_btr_per_km')
+
+ # https://stackoverflow.com/a/1384565/7728169
+ # Note: this will drop any fractional seconds, which is usually okay
+ dt = time.strftime('%H:%M:%S', time.gmtime(dt_per_km * resolution))
+ btr_dt = time.strftime(
+ '%H:%M:%S', time.gmtime(dt_btr_per_km * resolution))
+
+ options = dict(config_dt="'{}'".format(dt),
+ config_btr_dt="'{}'".format(btr_dt))
+ self.update_namelist_at_runtime(options)
+
iteration_count = config.getint('ssh_adjustment', 'iterations')
adjust_ssh(variable='landIcePressure', iteration_count=iteration_count,
step=self)
diff --git a/latest/_modules/compass/ocean/tests/isomip_plus/viz.html b/latest/_modules/compass/ocean/tests/isomip_plus/viz.html
index 27f8512caa..83bae6fbad 100644
--- a/latest/_modules/compass/ocean/tests/isomip_plus/viz.html
+++ b/latest/_modules/compass/ocean/tests/isomip_plus/viz.html
@@ -115,6 +115,8 @@
Source code for compass.ocean.tests.isomip_plus.viz
import os
+import matplotlib.pyplot as plt
+import numpy as np
import xarray
from mpas_tools.io import write_netcdf
@@ -218,6 +220,24 @@ Source code for compass.ocean.tests.isomip_plus.viz
plotter.plot_horiz_series(delssh, 'delssh', 'delssh', True,
cmap='cmo.curl', vmin=-1, vmax=1)
+ wct = dsOut.ssh + dsMesh.bottomDepth
+ idx_thin = np.logical_and(wct[0, :] > 1e-1,
+ wct[0, :] < 1)
+ wct_thin = wct[:, idx_thin]
+ wct_mean = wct_thin.mean(dim='nCells').values
+ time = dsOut.daysSinceStartOfSim.values
+ fig = plt.figure()
+ plt.plot(time, wct_mean, '.')
+ fig.set_xlabel('Time (days)')
+ fig.set_ylabel('Mean thickness of thin film (m)')
+ plt.savefig('wct_thin_t.png')
+ plt.close()
+
+ plotter.plot_horiz_series(wct, 'H', 'H',
+ True, vmin=min_column_thickness + 1e-10,
+ vmax=700, cmap_set_under='r',
+ cmap_scale='log')
+
if os.path.exists(f'{sim_dir}/timeSeriesStatsMonthly.0001-01-01.nc'):
ds = xarray.open_mfdataset(
'{}/timeSeriesStatsMonthly*.nc'.format(sim_dir),
diff --git a/latest/developers_guide/ocean/generated/compass.ocean.tests.isomip_plus.isomip_plus_test.IsomipPlusTest.html b/latest/developers_guide/ocean/generated/compass.ocean.tests.isomip_plus.isomip_plus_test.IsomipPlusTest.html
index 3c21311aaf..6878f99c43 100644
--- a/latest/developers_guide/ocean/generated/compass.ocean.tests.isomip_plus.isomip_plus_test.IsomipPlusTest.html
+++ b/latest/developers_guide/ocean/generated/compass.ocean.tests.isomip_plus.isomip_plus_test.IsomipPlusTest.html
@@ -140,6 +140,7 @@ compass.ocean.tests.isomip_plus.isomip_plus_test.IsomipPlusTestfloat) – The horizontal resolution (km) of the test case
experiment ({'Ocean0', 'Ocean1', 'Ocean2'}) – The ISOMIP+ experiment
vertical_coordinate (str) – The type of vertical coordinate (z-star, z-level, etc.)
tidal_forcing (bool) – Whether the case has tidal forcing
time_varying_forcing (bool) – Whether the run includes time-varying land-ice forcing
thin_film_present (bool) – Whether a thin film is present under land ice
planar (bool, optional) – Whether the test case runs on a planar or a spherical mesh