Skip to content

*+Revise non-Boussinesq find_coupling_coef calcs#373

Merged
marshallward merged 2 commits intoNOAA-GFDL:dev/gfdlfrom
Hallberg-NOAA:nonBous_coupling_coef
Jun 8, 2023
Merged

*+Revise non-Boussinesq find_coupling_coef calcs#373
marshallward merged 2 commits intoNOAA-GFDL:dev/gfdlfrom
Hallberg-NOAA:nonBous_coupling_coef

Conversation

@Hallberg-NOAA
Copy link
Copy Markdown
Member

Restructure one of the find_coupling_coef calculations to draw out the stress-magnitude terms, in preparation for future steps to reduce the dependency on the Boussinesq reference density. Using a value of VERT_FRICTION_ANSWER_DATE that is below 20230601 recovers the previous answers with non-Boussinesq test cases, but this is irrelevant for Boussinesq test cases. This updated code is mathematically equivalent to the previous expressions but it does change answers at roundoff in non-Boussinesq cases for recent answer dates. There are modifications to some comments in MOM_parameter_doc files.

  Restructure one of the find_coupling_coef calculations to draw out the
stress-magnitude terms, in preparation for future steps to reduce the dependency
on the Boussinesq reference density.  Using a value of VERT_FRICTION_ANSWER_DATE
that is below 20230601 recovers the previous answers with non-Boussinesq test
cases, but this is irrelevant for Boussinesq test cases.  This updated code is
mathematically equivalent to the previous expressions but it does change answers
at roundoff in non-Boussinesq cases for recent answer dates.  There are
modifications to some comments in MOM_parameter_doc files.
@codecov
Copy link
Copy Markdown

codecov Bot commented May 26, 2023

Codecov Report

Merging #373 (f4534c0) into dev/gfdl (1faa9ab) will increase coverage by 0.00%.
The diff coverage is 40.00%.

❗ Current head f4534c0 differs from pull request most recent head 206a7af. Consider uploading reports for the commit 206a7af to get more accurate results

@@            Coverage Diff            @@
##           dev/gfdl     #373   +/-   ##
=========================================
  Coverage     38.28%   38.29%           
=========================================
  Files           269      269           
  Lines         76018    76022    +4     
  Branches      13982    13983    +1     
=========================================
+ Hits          29107    29109    +2     
- Misses        41682    41684    +2     
  Partials       5229     5229           
Impacted Files Coverage Δ
...c/parameterizations/vertical/MOM_vert_friction.F90 58.59% <40.00%> (-0.13%) ⬇️

... and 1 file with indirect coverage changes

📣 We’re building smart automated test selection to slash your CI/CD build times. Learn more

@Hallberg-NOAA Hallberg-NOAA added the answer-changing A change in results (actual or potential) label May 26, 2023
@marshallward
Copy link
Copy Markdown
Member

marshallward commented Jun 8, 2023

Gaea regression: https://gitlab.gfdl.noaa.gov/ogrp/MOM6/-/pipelines/19434 ❌ (✔️) 🟡

As expected, answers have changed in the nonBous_global test.

@marshallward marshallward merged commit e672b98 into NOAA-GFDL:dev/gfdl Jun 8, 2023
@Hallberg-NOAA Hallberg-NOAA deleted the nonBous_coupling_coef branch November 8, 2023 23:15
dhruvbalwada pushed a commit to dhruvbalwada/MOM6 that referenced this pull request Jan 31, 2026
Updates the NUOPC cap to write the stochastic physics package's restart file in the run directory whenever using the Stochastic GM+E parameterization.

When restarting, enables the stochastic physics package to read a restart file from a name/location other than the default. If the new MOM6 namelist parameter CESM_RESTFILE is True then it uses the mom6 restart filename but with .r_stoch. in place of .r.. If CESM_RESTFILE is False (default) it simply uses whatever the default restart file name is in the stochastic physics package. (Edit: This PR no longer uses the CESM_RESTFILE parameter. It just reads the restart file directly from the run directory.) A separate PR to ESCOMP/stochastic_physics updates that package to allow the restart file name to change from its default value, which was originally hard-coded.

This partially addresses NOAA-GFDL#360. Changes will be needed in MOM_interface to make the whole thing work; specifically, we should always set CESM_RESTFILE = True in user_nl_mom, and when CONTINUE_RUN is true we need to modify the stochini parameter in the &nam_stochy section of input.nml to take the value .true..

This PR takes an alternate approach to NOAA-GFDL#373, which required the restart file name to be set in the stochastic physics namelist. I've closed NOAA-GFDL#373 in preference to this approach.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

answer-changing A change in results (actual or potential)

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants