Skip to content

Cherry pick to add an option to calculate omega in non-hydrostatic runs similar to hydrostatic method#344

Merged
bensonr merged 2 commits into
NOAA-GFDL:dev/emcfrom
XiaqiongZhou-NOAA:cherry_pick_omga
Jul 3, 2024
Merged

Cherry pick to add an option to calculate omega in non-hydrostatic runs similar to hydrostatic method#344
bensonr merged 2 commits into
NOAA-GFDL:dev/emcfrom
XiaqiongZhou-NOAA:cherry_pick_omga

Conversation

@XiaqiongZhou-NOAA
Copy link
Copy Markdown
Contributor

@XiaqiongZhou-NOAA XiaqiongZhou-NOAA commented Jun 14, 2024

Description
This update cherry-picks from Release 042022 (#184). A new parameter, pass_full_omega_to_physics_in_non_hydrostatic_mode, has been added. When set to true, the omega calculation will match that of hydrostatic runs, using the formula omega = dp/dt. The default value is false, where the calculation remains omega = w * dp/dz as before.

How Has This Been Tested?

This update is tested on Hercules with both intel and gnu compiler and all UFS regression tests produced identical results with default setting.

Checklist:

Please check all whether they apply or not

  • My code follows the style guidelines of this project
  • I have performed a self-review of my own code
  • I have commented my code, particularly in hard-to-understand areas
  • I have made corresponding changes to the documentation
  • My changes generate no new warnings
  • Any dependent changes have been merged and published in downstream modules

…ydrostatic mode when 'pass_full_omega_to_physics_in_non_hydrostatic_mode' is set to true
Copy link
Copy Markdown
Contributor

@lharris4 lharris4 left a comment

Choose a reason for hiding this comment

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

I think this looks good. Thank you for your careful implementation.

if (.not. Atm(n)%flagstruct%hydrostatic .and. .not. Atm(n)%flagstruct%pass_full_omega_to_physics_in_non_hydrostatic_mode) then
Atm(n)%omga(isc:iec,jsc:jec,1:npz) = Atm(n)%delp(isc:iec,jsc:jec,1:npz) / Atm(n)%delz(isc:iec,jsc:jec,1:npz) * Atm(n)%w(isc:iec,jsc:jec,1:npz)
if(Atm(n)%flagstruct%nf_omega>0) then
call del2_cubed(&
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.

Good use of existing API.

@zach1221
Copy link
Copy Markdown

zach1221 commented Jul 2, 2024

@laurenchilutti testing is complete on WM PR#2327. Please feel free to merge this cubed-sphere PR.
@FernandoAndrade-NOAA

@zach1221
Copy link
Copy Markdown

zach1221 commented Jul 2, 2024

Copying @bensonr as well. Rusty, if you're able to merge, this PR is ready.

@bensonr bensonr merged commit 577fd5e into NOAA-GFDL:dev/emc Jul 3, 2024
laurenchilutti added a commit to laurenchilutti/GFDL_atmos_cubed_sphere that referenced this pull request Jul 11, 2024
…lar to hydrostatic mode when 'pass_full_omega_to_physics_in_non_hydrostatic_mode' is set to true (NOAA-GFDL#344)"

This reverts commit 577fd5e.
@laurenchilutti laurenchilutti mentioned this pull request Jul 11, 2024
6 tasks
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.

4 participants