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 calculation of Pacific Equatorial Undercurrent transport #273

Merged
merged 2 commits into from
Jun 1, 2020

Conversation

ashao
Copy link
Collaborator

@ashao ashao commented May 31, 2019

Huan Gao and Rong Zhang noted that the transport in the equatorial
undercurrent should be around ~24 Sv in the annual mean, but in
the actual diagnostic after postprocessing, it was around -8Sv. This
was traced down to errors in the strait_transport_timeseries.

This commit fixes the logical and indexing errors in that script namely
in the way that the z-index to integrate to and the zeroing out of any
deeper layers. For the year 650 of the CM4 PI-Control experiment, the
routine now returns an annually averaged transport of 25 Sv.

Huan Gao and Rong Zhang noted that the transport in the equatorial
undercurrent should be around ~24 Sv in the annual mean, but in
the actual diagnostic after postprocessing, it was around -8Sv. This
was traced down to errors in the strait_transport_timeseries.

This commit fixes the logical and indexing errors in that script namely
in the way that the z-index to integrate to and the zeroing out of any
deeper layers. For the year 650 of the CM4 PI-Control experiment, the
routine now returns an annually averaged transport of 25 Sv.
@ashao
Copy link
Collaborator Author

ashao commented May 31, 2019

I tested the changes here manually by calling sum_transport_in_straits directly and pointing it at output in /archive/oar.gfdl.cmip6/CM4/warsaw_201710_om4_v1.0.1/CM4_piControl_C/gfdl.ncrc4-intel16-prod-openmp/history/tmp_650/.

However, I'm not sure if there's an alternate testing pipeline for the refineDiag part of the postprocessing in general to make sure that it didn't break anything else. Do you know @jkrasting?

@adcroft
Copy link
Member

adcroft commented Jul 12, 2019

@jkrasting Do you know if this has been tested?

@adcroft
Copy link
Member

adcroft commented Jul 23, 2019

@jkrasting Bump

@adcroft
Copy link
Member

adcroft commented Sep 13, 2019

This PR is labelled "draft" presumably because it wasn't tested in the workflow. I'm still waiting to hear whether this has been tested. @jkrasting

@adcroft adcroft marked this pull request as ready for review June 1, 2020 22:41
@adcroft adcroft merged commit 3c02903 into NOAA-GFDL:dev/gfdl Jun 1, 2020
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.

2 participants