Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix surface boundary conditions for TEOS-10 #812

Draft
wants to merge 14 commits into
base: master
Choose a base branch
from

Conversation

mjlosch
Copy link
Member

@mjlosch mjlosch commented Feb 28, 2024

What changes does this PR introduce?

This a bug fix towards completing the TEOS-10 implementation

What is the current behaviour?

When using TEOS-10 ('eosType='TEOS10', in data&PARMS01), the variable THETA is interpreted as conservative temperatures (CT), which introduces a bias at in bulk formulae at the surface. Therefore CT needs to be converted to potential (or at the surface: in-situ) temperature before using it in bulk formulae.

What is the new behaviour

add a new field seaSurfTemp that is computed once to be used in all bulk formulae. Convert this from CT to potential/in-situ temperature if 'eosType='TEOS10'`.

Does this PR introduce a breaking change?

  • changes results when using TEOS10
  • moved sstExtrapol from pkg/exf to the main namelist data&PARMS01, I doubt that anyone uses this.

Other information:

This is a draft to leave room for discussion as I may have missed a few points.
Also, I moved the extrapolation parameter sstExtrapol from pkg/exf to the main namelist data&PARMS01, so that it is available to all code. We might want to discuss a new name for this or find a good way of retiring it from the namelist without changing the name.

Suggested addition to tag-index

T.B.D.

mjlosch and others added 14 commits December 21, 2023 18:49
- convert conservative temperature to potential temperature to use as
SST, when eosType='TEOS10'
- this is a very rough first pitch and needs to be polished and
improved for performance especially for the new functions for
temperature conversions which are a poorly disguised copy of MOM6.
- some formatting of the new functions (incomplete)
replace functions by subroutine, documenation
- convert conservative temperature to potential temperature to use as
SST, when eosType='TEOS10'
- this is a very rough first pitch and needs to be polished and
improved for performance especially for the new functions for
temperature conversions which are a poorly disguised copy of MOM6.
- some formatting of the new functions (incomplete)
replace functions by subroutine, documenation
In the next step, the runtime parameter sstExtrapol will be moved from
data.exf to the global paramter list
without properly retiring it from the EXF_NML_01
@jrscott jrscott mentioned this pull request May 7, 2024
@mjlosch mjlosch requested a review from jrscott May 8, 2024 15:25
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.

None yet

1 participant