Skip to content

Conversation

@vanroekel
Copy link
Contributor

Adds the ability to calculate ocean heat content during a forward model run

It computes heat content over standard layers
0-700m
700-2000m
2000m - bot
full Depth

@vanroekel
Copy link
Contributor Author

I've tested this in a brief 60to30 stand alone test case. Results look reasonable. Here is the sfctobottom anomaly (units of 10^22 J)

Screen Shot 2020-02-25 at 3 23 15 PM

@vanroekel
Copy link
Contributor Author

vanroekel commented Feb 25, 2020

I've added four arrays (nCells dimension)

oceanHeatContentSFCtoBot
oceanHeatContentSfcto700m
oceanHeatContent700mto2000m
oceanHeatContent2000mtoBot

@milenaveneziani
Copy link
Contributor

@vanroekel: we normally compute the OHC from the layerWeigthedAvg AM (and then define the depth ranges over which to do the accumulated quantity through the MPAS-Analysis config file). So I am not sure why we need this AM? Does the layerWeigthedAvg AM have a problem?

@vanroekel
Copy link
Contributor Author

my read of that AM is it can do the globally averaged OHC, this returns a horizontal map. This PR is needed to return the map but more importantly for work on a LDRD-DR effort where we want to compute the 2-D maps at high frequency while the model runs and do statistics on the data.

@milenaveneziani
Copy link
Contributor

Indeed the layerWeightedAvg AM computes quantities over specified regions (including the global region). I was thinking that we could potentially modify it, to output the T and S content globally over specified depth ranges, but perhaps better to keep things separate.

So I assume this AM does what your postprocessing script does to get the global maps?

@mark-petersen mark-petersen self-assigned this Mar 11, 2020
@mark-petersen mark-petersen force-pushed the ocean/addOceanHeatContentAM branch from a251a7f to 66e5b41 Compare March 11, 2020 20:51
@mark-petersen
Copy link
Contributor

rebased. Added this AM and output variables to analysis test in nightly regression suite.

On above comments, it sounds like this AM is sufficiently different from others. If it is useful, it's fine to have some redundancy between AMs and MPAS-Analysis. We can also write these out at higher frequency.

@mark-petersen
Copy link
Contributor

This can can be merged in next week. It requires E3SM flag additions, so will be put in an E3SM PR for that.

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.

Turned this AM on in QU240 analysis test. Passes gnu/optimized and intel/debug.

vanroekel and others added 3 commits May 27, 2020 13:34
Adds the ability to calculate ocean heat content in the forward model

It computes heat content over standard layers
0-700m
700-2000m
2000m - bot
full Depth
@mark-petersen mark-petersen force-pushed the ocean/addOceanHeatContentAM branch from 66e5b41 to 3a5dcab Compare May 29, 2020 23:18
@mark-petersen
Copy link
Contributor

Rebased. Passes nightly regression with gnu and intel/debug, including this AM in the analysis test.

@mark-petersen
Copy link
Contributor

@xylar and @milenaveneziani I would like to merge this. Do you have any comments?

@xylar
Copy link
Collaborator

xylar commented May 30, 2020

Sorry, I didn't realize my review had been requested. I don't have any input so will take myself off as a reviewer.

@xylar xylar removed their request for review May 30, 2020 13:45
@mark-petersen mark-petersen removed the request for review from milenaveneziani June 1, 2020 15:39
@mark-petersen mark-petersen merged commit 8d14c0e into MPAS-Dev:ocean/develop Jun 1, 2020
mark-petersen added a commit that referenced this pull request Jun 1, 2020
Adds normalVelocity arrays to eddy products #544
Adds Ocean Heat Content AM #455
Ocean/remove unused modules #569
jonbob added a commit to E3SM-Project/E3SM that referenced this pull request Jun 2, 2020
Update mpas-source: analysis member update, remove unused modules

This PR brings in a new mpas-source submodule with changes only to the ocean,
including several updates to analysis members and removing modules that are not
used in E3SM simulations:
* Adds normalVelocity arrays to eddy products (MPAS-Dev/MPAS-Model#544)
* Adds Ocean Heat Content AM (MPAS-Dev/MPAS-Model#455)
* Ocean/remove unused modules (MPAS-Dev/MPAS-Model#569)
It includes corresponding changes to E3SM namelists and streams files for
MPAS-Ocean. COMPASS changes associated with these PRs are also included.

[NML]
[BFB]
jonbob added a commit to E3SM-Project/E3SM that referenced this pull request Jun 3, 2020
…3613)

Update mpas-source: analysis member update, remove unused modules

Re-merged after fixes for debug failures, excess output

This PR brings in a new mpas-source submodule with changes only to the ocean,
including several updates to analysis members and removing modules that are not
used in E3SM simulations:
* Adds normalVelocity arrays to eddy products (MPAS-Dev/MPAS-Model#544)
* Adds Ocean Heat Content AM (MPAS-Dev/MPAS-Model#455)
* Ocean/remove unused modules (MPAS-Dev/MPAS-Model#569)
It includes corresponding changes to E3SM namelists and streams files for
MPAS-Ocean. COMPASS changes associated with these PRs are also included.

[NML]
[BFB]
jonbob added a commit to E3SM-Project/E3SM that referenced this pull request Jun 4, 2020
Update mpas-source: analysis member update, remove unused modules

This PR brings in a new mpas-source submodule with changes only to the ocean,
including several updates to analysis members and removing modules that are not
used in E3SM simulations:
* Adds normalVelocity arrays to eddy products (MPAS-Dev/MPAS-Model#544)
* Adds Ocean Heat Content AM (MPAS-Dev/MPAS-Model#455)
* Ocean/remove unused modules (MPAS-Dev/MPAS-Model#569)
It includes corresponding changes to E3SM namelists and streams files for
MPAS-Ocean. COMPASS changes associated with these PRs are also included.

[NML]
[BFB]
mark-petersen added a commit to mark-petersen/MPAS-Model that referenced this pull request Jan 11, 2021
…ean/develop

Adds Ocean Heat Content AM MPAS-Dev#455

Adds the ability to calculate ocean heat content during a forward model
run

It computes heat content over standard layers
0-700m
700-2000m
2000m - bot
full Depth
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