Skip to content

Conditional registration of FrictWorkMax with MEKE#946

Merged
adcroft merged 2 commits into
mom-ocean:dev/gfdlfrom
marshallward:frictworkmax
Jul 10, 2019
Merged

Conditional registration of FrictWorkMax with MEKE#946
adcroft merged 2 commits into
mom-ocean:dev/gfdlfrom
marshallward:frictworkmax

Conversation

@marshallward
Copy link
Copy Markdown
Collaborator

The FrictWorkMax diagnostic requires that MEKE be enabled and that MEKE%MEKE field be allocated and defined. It is currently possible to register this diagnostic, even when MEKE is disabled.

This patch prevents the registration of FrictWorkMax when MEKE is turned off.

The FrictWorkMax diagnostic requires that MEKE be enabled and that
MEKE%MEKE field be allocated and defined.  It is currently possible to
register this diagnostic, even when MEKE is disabled.

This patch prevents the registration of FrictWorkMax when MEKE is turned
off.
@marshallward
Copy link
Copy Markdown
Collaborator Author

GitLab testing here: https://gitlab.gfdl.noaa.gov/ogrp/MOM6/pipelines/8446

@marshallward
Copy link
Copy Markdown
Collaborator Author

GitLab tests have passed.

@adcroft adcroft merged commit 97afc1d into mom-ocean:dev/gfdl Jul 10, 2019
@marshallward marshallward deleted the frictworkmax branch September 3, 2019 15:42
Hallberg-NOAA added a commit to Hallberg-NOAA/MOM6 that referenced this pull request Sep 6, 2025
…n#946)

* Added a 7th-, 5th- and 3rd-order WENO schemes to reconstruct the PV and added the Koren scheme to reconstruct the kinetic energy in MOM_CoriolisAdv:

- Set CORIOLIS_SCHEME = "WENOVI7TH_PV_ENSTRO", "WENOVI5TH_PV_ENSTRO", :WENOVI3RD_PV_ENSTRO" to use a 7th-, 5th-, and 3rd-order WENO scheme, respectively, to reconstruct PV.
- Set KE_SCHEME = "KE_UP3" to use a third-order upwind scheme with the Koren flux limiter to reconstruct the kinetic energy. Set KE_USE_LIMITER = False (default is True) to not use the flux limiter. The KE_UP3 scheme is recommended to be used along with the WENO Coriolis schemes.
- Set WENO_VELOCITY_SMOOTH = True to use the velocity field to estimate the smoothness indicator for WENO. This scheme tends to yield higher-order of accuracy and thus more energetic eddy fields.
- The WENO schemes conduct a thickness-weighted reconstruction for PV to avoid PV singularities.
- The order of WENO is reduced to 5th, 3rd, and 1st subsequently when it approaches lateral boundaries
- The halo size is 4 and 5 for the 7th- and 5th-order WENO schemes. When they are used, parameters USE_WENO and WENO_stencil are passed to MOM_dynamics_split_RK2 to increase the halo update size. Only MOM_dynamics_split_RK2 is changed now. The other RK2 modules will need to be updated later.

* Replaced the WENO_stencil by the CorAdv_stencil function that is passed to MOM_dynamics_split_RK2

* Debug

* Replaced a few ratios to variables. Corrected the unit of some variables

* Changed adding area_tiny and h_tiny to taking the max

* Removed the dimensional constant in WENO weighting

* Added stencil in $OMP(shared). Added additional comments

* Removed trailing space

* Changed the zero check if statement for fac_fn. Changed Jsq, Jeq, Isq, Ieq to private variables

* Addressed couple of rotational symmetry issues

* Removed commented lines in RK2 and added WENO stencil information to RK2B

* Added Is_q, Ie_q, Js_q, Je_q to set bounds for Coriolis schemes

* Enlarged halo update size for WENO schemes in dynamic_unsplit modules

* Removed trailing space

* Moved the update of h to a group pass in RK2 module

* Removed unnecessary lines

---------

Co-authored-by: Alistair Adcroft <Alistair.Adcroft@noaa.gov> and Robert Hallberg <Robert.Hallberg@noaa.gov>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants