Skip to content

Commit

Permalink
Rework documentation of glc_renormalize_smb
Browse files Browse the repository at this point in the history
Rework documentation of glc_renormalize_smb to be consistent with the
changes in ESCOMP/CMEPS#495
  • Loading branch information
billsacks committed Aug 9, 2024
1 parent 489304b commit ed947c3
Show file tree
Hide file tree
Showing 2 changed files with 10 additions and 16 deletions.
12 changes: 4 additions & 8 deletions doc/source/clm-cism-coupling.rst
Original file line number Diff line number Diff line change
Expand Up @@ -43,14 +43,10 @@ coupling becomes important. There are still some unhandled edge cases in this re
significant case being when CLM dictates an amount of glacial melt that exceeds the
available ice in a CISM grid cell), but we do take pains to achieve conservation in most
cases. An important mechanism to achieve this conservation is via a global renormalization
step done when mapping SMB from CLM to CISM, as described in :numref:`remapping_smb`. **By
default, this renormalization is done when running with a two-way-coupled ice sheet that
sends fluxes to other components (which is typically true in this case of an evolving ice
sheet), but is turned off in other cases. This leads to small differences in the remapped
SMB field in runs with an evolving vs. non-evolving ice sheet. This default can be
overridden via the driver namelist variable,** ``glc_renormalize_smb``, **keeping in mind
that setting this to** ``off`` **will break conservation for a configuration with an
evolving, two-way-coupled ice sheet.**
step done when mapping SMB from CLM to CISM, as described in :numref:`remapping_smb`.
(This renormalization is done by default, but can be turned off by setting the driver
namelist variable, ``glc_renormalize_smb`` to ``off``, keeping in mind that this will
break conservation in a fully-coupled configuration.)

Evolving ice sheet with one-way (diagnostic) coupling
-----------------------------------------------------
Expand Down
14 changes: 6 additions & 8 deletions doc/source/t-compsets.rst
Original file line number Diff line number Diff line change
Expand Up @@ -100,14 +100,12 @@ averages will appear in the cpl/hist directory within your archive space, with n
files are about 7 MB per year.)

A T compset run that later uses these coupler history files as forcing should give
*nearly* identical CISM results as the original run, **as long as you ensure that SMB
renormalization is done (or not done) in both cases, as described in**
:numref:`clm-cism-coupling`. **(For example, if the coupler history files were generated
from a fully-coupled case with an evolving, two-way-coupled ice sheet, then in the T case
you should set** ``glc_renormalize_smb = 'on'`` **in** ``user_nl_cpl`` **in order to
reproduce the results of the fully-coupled case.)** Small differences arise because these
forcing files are written with single precision, leading to roundoff error on the order of
10\ :sup:`-7`. If you need more precision, you can add the following to ``user_nl_cpl``:
*nearly* identical CISM results as the original run, as long as you ensure that SMB
renormalization is done (or not done) in both cases, as described in
:numref:`clm-cism-coupling` (based on the setting of ``glc_renormalize_smb``). Small
differences arise because these forcing files are written with single precision, leading
to roundoff error on the order of 10\ :sup:`-7`. If you need more precision, you can add
the following to ``user_nl_cpl``:

.. code-block:: console
Expand Down

0 comments on commit ed947c3

Please sign in to comment.