Skip to content

Fix for LES tke in real-data cases - #1499

Merged
mgduda merged 1 commit into
MPAS-Dev:hotfix-v8.4.2from
dudhia:lesfix-8.4.1
Sep 3, 2026
Merged

Fix for LES tke in real-data cases#1499
mgduda merged 1 commit into
MPAS-Dev:hotfix-v8.4.2from
dudhia:lesfix-8.4.1

Conversation

@dudhia

@dudhia dudhia commented Aug 5, 2026

Copy link
Copy Markdown
Collaborator

For the prognostic tke option, initial tke was zero. This leads to no tke generation because source terms depend on non-zero tke through the K diffusion coefficients. The fix is to set a small minimum value (currently chosen as 0.01 m2/s) on the K coefficients. Two lines are added in the dissipation module.

@dudhia

dudhia commented Aug 5, 2026

Copy link
Copy Markdown
Collaborator Author

Not sure what branch to merge to.

@mgduda
mgduda changed the base branch from master to hotfix-v8.4.2 August 5, 2026 19:16
@mgduda
mgduda requested a review from skamaroc August 5, 2026 19:17

@skamaroc skamaroc left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This change is necessary as stated in the pull request and the commit message. The lower bound on the tke is consistent with practices in other models (e.g. WRF) and is in the correct location in the MPAS code base.

@mgduda
mgduda merged commit 6845943 into MPAS-Dev:hotfix-v8.4.2 Sep 3, 2026
mgduda added a commit that referenced this pull request Sep 4, 2026
This merge addresses several issues in the MPAS-Atmosphere model and in the MPAS
infrastructure. Specific changes include:

 * Fix an issue with the naming of MPAS's internal ESMF time-keeping library in
   CMake builds of MPAS. (PR #1506)

 * Fix an issue with restart runs in MPAS-Atmosphere on GPUs. (PR #1503)

 * Address a lack of TKE in MPAS-Atmosphere with the 1.5-order prognostic TKE
   option when there is no TKE in the initial conditions (for example, in
   real-data cases). (PR #1499)

 * Correct the ordering of the 'lbc_ni' and 'lbc_nr' scalars in the atmosphere
   core's Registry.xml file. Note that unless LBCs are provided for 'ni' and
   'nr', the previously incorrect order of 'lbc_ni' and 'lbc_nr' would not
   affect results, as both fields would be filled with zero values. (PR #1498)

 * Change the default value of the config_number_of_sub_steps namelist option
   from 2 to 4 in the atmosphere core's Registry.xml file to address an issue in
   convection-permitting simulations where noise can sometimes appear in the jet
   streams using the previous default value of 2. The noise, when present, also
   produces unphysical features in horizontal kinetic energy spectra at high
   wavenumbers. (PR #1496)

   Note: Increasing the number of sub-steps from 2 to 4 can be expected to
   increase the computational cost of the dynamics by ~10% and the
   computational cost of a full-physics simulation by ~5%.

 * Work around an issue with continuation lines in instances of the
   PACKAGE_LOGIC_PRINT macro that prevented successful compilation of
   MPAS-Atmosphere with the 'nvfortran' compiler. (PR #1495)

* hotfix-v8.4.2:
  Rename bundled ESMF time-manager library output to esmf_time
  Fix for GPU-enabled restart runs
  Correct order of lbc_ni and lbc_nr in the atmosphere core's Registry.xml file
  add minimum value for Kh and Kv to allow for zero initial tke
  Change default value of atmosphere core config_number_of_sub_steps from 2 to 4
  Work around nvfortran issue with continuation lines in PACKAGE_LOGIC_PRINT macro
  Update version number to 8.4.2
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants