Skip to content

Conversation

@vanroekel
Copy link
Contributor

In the GM bolus calculation, there is a specified diffusivity and phase
speed (see Ferrari et al. 2010). Currently MPAS-O assumes these values
are fixed in time and space. This PR allows separate flags to enable
2D varying (+time) phase speed (based on the first baroclinic mode phase
speed) and a separate flag for 3D (+time) varying diffusivity.

In the GM bolus calculation, there is a specified diffusivity and phase
speed (see Ferrari et al. 2010).  Currently MPAS-O assumes these values
are fixed in time and space.  This PR allows separate flags to enable
2D varying (+time) phase speed (based on the first baroclinic mode phase
speed) and a separate flag for 3D (+time) varying diffusivity.
@vanroekel
Copy link
Contributor Author

This has currently been tested in 5 year normal forced G-cases. Results are here

https://portal.nersc.gov/project/acme/lvroekel/GM_20190605_GMtest_Kappa3d_vary_c2.anvil_yr5/

and a baseline is here

https://portal.nersc.gov/project/acme/lvroekel/20190420.GMPAS-IAF.T62_oEC60to30v3_gradientRI_0.6.edison_yr5/

A full core cycle test is in progress, currently at year 25

@vanroekel
Copy link
Contributor Author

@mark-petersen here is the 3d varying bolus kappa pr.

@mark-petersen
Copy link
Contributor

Design document is here:
https://acme-climate.atlassian.net/wiki/spaces/OCNICE/pages/980713736/ocean+GM+Parameter+function+of+depth

From the July 10 meeting notes:
The core cycle of the new GM bolus and control are done. You are welcome to discuss them at the cryo meeting tomorrow, I'll be out of town. The links are:
CONTROL
https://portal.nersc.gov/project/acme/lvroekel/GM600_20190624_GMtest_constant_intel.compy_yr_50_60
Test of variable GM Bolus
https://portal.nersc.gov/project/acme/lvroekel/GM600_20190617_GMtest_Kappa3d_vary_c2_intel.compy_yr_50_60
A few notes from @vanroekel

  1. Surface salinity biases are pretty nicely improved, continues to subsurface as well (pretty clear in sea floor maps near antarctica)
  2. MLD is better, perhaps a bit too deep.
  3. Sea ice is better as well.

From @milenaveneziani:
4) Also note that the vertical structure of temperature, salinity and PD in the SO is much improved (look at the transects), especially for transects at 187, 253, 280 and 318 longitude.
5) SO bottom temperature is much better, while bottom salinity is only marginally improved.

@vanroekel
Copy link
Contributor Author

removing the don't merge label, this is tested in a full core cycle and is flag protected -- default off.

@xylar
Copy link
Collaborator

xylar commented Jul 10, 2019

Design document is here: (Confluence link)

Having the design doc on Confluence isn't a good idea. It defeats the purpose of open development. Could this be copied as pdf or something to the local directory with design docs?

@pwolfram
Copy link
Contributor

Ideally we would have design documents with markdown or rst text so that we can readily render them via sphinx to a webpage documentation for the code some day, but I know this is extra work so it may not be worth it in the short term.

@xylar
Copy link
Collaborator

xylar commented Jul 10, 2019

@pwolfram, I agree. I think new design docs should be written that way (like @sbrus89 did for the one design doc we currently have) but, like you, I didn't think that was a reasonable thing to ask for in this case.

Removes the vertical averaging of the GM kappa bolus diffusivity.
Testing shows this functionality is not beneficial to simulations
Copy link
Contributor

@mark-petersen mark-petersen left a comment

Choose a reason for hiding this comment

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

Looks good. I merged locally tested regression suite, is bfb with previous on grizzly with gnu.

@mark-petersen mark-petersen merged commit 931ab3f into MPAS-Dev:ocean/develop Aug 16, 2019
mark-petersen added a commit that referenced this pull request Aug 16, 2019
In the GM bolus calculation, there is a specified diffusivity and phase
speed (see Ferrari et al. 2010). Currently MPAS-O assumes these values
are fixed in time and space. This PR allows separate flags to enable
2D varying (+time) phase speed (based on the first baroclinic mode phase
speed) and a separate flag for 3D (+time) varying diffusivity.
mark-petersen added a commit to mark-petersen/MPAS-Model that referenced this pull request Mar 24, 2020
…te_maint1.2

Adds options for 3d varying GM bolus and 2d varying phase speed MPAS-Dev#288

In the GM bolus calculation, there is a specified diffusivity and phase
speed (see Ferrari et al. 2010). Currently MPAS-O assumes these values
are fixed in time and space. This PR allows separate flags to enable
2D varying (+time) phase speed (based on the first baroclinic mode phase
speed) and a separate flag for 3D (+time) varying diffusivity.

* pr_test2:dd
  Removes 2D GM code and debug statements
  Adds 2d variable GM and fixes phase speed bug
  Adds options for 3d varying GM bolus and 2d varying phase speed
mark-petersen added a commit that referenced this pull request Mar 24, 2020
Update maint1.2: 3D varying GM

This update adds the following PRs from ocean/develop:
- Adds options for 3d varying GM bolus and 2d varying phase speed #288
- Add GM bolus eddy stats #339
- Fix threading issue in MPAS-O GM routine #376
- compute landIceMask using geometric_features #447
- bug fixes for the nonlocal source term in KPP #305
jonbob added a commit to E3SM-Project/E3SM that referenced this pull request Mar 27, 2020
#3513)

Update maint1.2: add 3D varying GM

This update brings in a new mpas-source submodule with the following PRs from
MPAS ocean/develop:
* Adds options for 3d varying GM bolus and 2d varying phase speed
  (MPAS-Dev/MPAS-Model#288)
* Add GM bolus eddy stats (MPAS-Dev/MPAS-Model#339)
* Fix threading issue in MPAS-O GM routine (MPAS-Dev/MPAS-Model#376)
* Compute landIceMask using geometric_features (MPAS-Dev/MPAS-Model#447)
* Bug fixes for the nonlocal source term in KPP (MPAS-Dev/MPAS-Model#305)
and from E3SM master:
* Change all ocean output files to single precision (E3SM #3360)

[non-BFB]
ashwathsv pushed a commit to ashwathsv/MPAS-Model that referenced this pull request Jul 21, 2020
In the GM bolus calculation, there is a specified diffusivity and phase
speed (see Ferrari et al. 2010). Currently MPAS-O assumes these values
are fixed in time and space. This PR allows separate flags to enable
2D varying (+time) phase speed (based on the first baroclinic mode phase
speed) and a separate flag for 3D (+time) varying diffusivity.
mark-petersen added a commit to mark-petersen/MPAS-Model that referenced this pull request Jan 11, 2021
In the GM bolus calculation, there is a specified diffusivity and phase
speed (see Ferrari et al. 2010). Currently MPAS-O assumes these values
are fixed in time and space. This PR allows separate flags to enable
2D varying (+time) phase speed (based on the first baroclinic mode phase
speed) and a separate flag for 3D (+time) varying diffusivity.
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.

5 participants