-
Notifications
You must be signed in to change notification settings - Fork 19
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Stochastic GM+E #280
Stochastic GM+E #280
Conversation
This commit adds a Stochastic Kinetic Energy Backscatter Scheme (SKEBS) where the amplitude of the backscatter is based on either the GM work rate and/or the viscous work rate (FrictWork). Each of these can be multiplied by a coefficient so that, e.g. the backscatter rate could be 50% of the GM work rate plus 20% of the viscous work rate. The vertical structure of the backscatter rate associated with GM has vertical structure given by the EBT profile. This code was developed starting from Phil Pegion's branch, and it builds on the stochastic physics package (external) developed by Phil Pegion, Niraj Agarwal, and collaborators. This package allows the length and time scales of the backscatter to be set via namelist parameters. This commit may break the stochastic EPBL and SPPT schemes developed by P. Pegion.
@iangrooms, can you please resolve the conflicts listed below? |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
It's passing aux_mom tests.
@gustavo-marques Can we tag this version of |
@mnlevy1981; yes, I think it's fine to tag this version to work out the requirement for the external |
This PR adds a Stochastic Kinetic Energy Backscatter Scheme (SKEBS) where the amplitude of the backscatter is based on either the GM work rate and/or the viscous work rate (FrictWork). Each of these can be multiplied by a coefficient so that, e.g. the backscatter rate could be 50% of the GM work rate plus 20% of the viscous work rate. The vertical structure of the backscatter rate associated with GM is given by the EBT profile. The component that uses FrictWork has not been extensively tested, but it does run stably for at least few years.
This code was developed starting from Phil Pegion's branch. This PR may break the stochastic EPBL and SPPT schemes developed by @pjpegion.
If you don't use the external
stochastic_physics
package developed by Phil Pegion, Niraj Agarwal, and collaborators, the code calls a stub inconfig_src/externals/stochastic_physics/
that does nothing. To use the external package, replace the contents ofconfig_src/externals/stochastic_physics/
with the code at the link above. You can control the time and length scales of the random pattern through a namelist included ininput.nml
. An example is included below.New MOM_input parameters:
New available diagnostics:
Namelist example to control the stochastic physics package:
The only relevant parameters for stochastic backscatter are
OCNSKEB
which should be 1. Backscatter amplitude should be controlled by the MOM_input parameters above.OCNSKEB_LSCALE
is the characteristic wavelength of the backscatter in metersOCNSKEB_TAU
is the decorrelation time of the noise in sectionsISEED_OCNSKEB
is a seed for the random pattern generator. If using ensembles, this seed value should be changed across ensemble members.