Skip to content

Commit

Permalink
update forcing_method 1 in chap_cases
Browse files Browse the repository at this point in the history
  • Loading branch information
grantfirl committed Aug 14, 2024
1 parent 0a58e4f commit c874ad7
Showing 1 changed file with 23 additions and 1 deletion.
24 changes: 23 additions & 1 deletion scm/doc/TechGuide/chap_cases.rst
Original file line number Diff line number Diff line change
Expand Up @@ -351,7 +351,8 @@ Optional arguments:

#. ``--save_comp (-sc)``: flag to create UFS reference file for comparison

#. ``--use_nearest (-near)``: flag to indicate using the nearest UFS history file gridpoint
#. ``--use_nearest (-near)``: flag to indicate using the nearest UFS history file gridpoint for calculation
of forcing; only valid for use with -fm=1 or -fm=3

#. ``--forcing_method (-fm)``: method used to calculate forcing (1=total tendencies from UFS dycore,
2=advective terms calculated from UFS history files, 3=total time tendency terms calculated), default=2
Expand All @@ -363,6 +364,27 @@ Optional arguments:

#. ``--geostrophic (-geos)``: flag to turn on geostrophic wind forcing

Notes Regarding Implemented Forcing Methods

The ``--forcing_method`` option hides some complexity that should be understood when running this script since
each method has a particular use case and produces potentially very different forcing terms. Forcing method 1 is
designed to be used in concert with the three-dimensional UFS. I.e., the UFS must be run with diagnostic tendencies
activated so that the `nophysics` term is calculated and output for all grid points. This diagnostic term
represents the tendency produced for each state variable by the UFS between calls to the "slow" physics. This
includes the tendency due to advection, but also any tendencies due to other non-physics processes, e.g. "fast"
physics, coupling to external components, data assimilation, etc. Within the SCM, this diagnostic is used as the
forcing term for each state variable. Although one can achieve results as close as possible between a UFS column
and the single column model using this method, it will NOT be bit-for-bit for many reasons. Some of these reasons
include: diagnostic output is not typically instantaneous for every timestep, the UFS' vertical coordinate is
semi-Lagrangian and includes a remapping step as the surface pressure changes for each column, whereas the SCM
uses a Eulerian vertical coordinate without the vertical remapping step, and some interpolation happens in the
UFS_case_gen.py script due to the UFS initial conditions and history files using different grids. This method
can only be used when the UFS has been configured and run with the anticipation of running the SCM using this
forcing method afterward because it requires considerable extra disk space for the additional output.




.. _`ufsforcingensemblegenerator`:

UFS_forcing_ensemble_generator.py
Expand Down

0 comments on commit c874ad7

Please sign in to comment.