Skip to content

Fix ISOMIP+ initial state and top drag coefficient - #998

Merged
xylar merged 2 commits into
MPAS-Dev:mainfrom
xylar:fix-isomip-plus
Sep 28, 2026
Merged

xylar merged 2 commits into
MPAS-Dev:mainfrom
xylar:fix-isomip-plus

Conversation

@xylar

@xylar xylar commented Sep 28, 2026 •

Copy link
Copy Markdown
Collaborator

The ISOMIP+ tasks fail in initial_state with current xarray, which no longer accepts a Dataset in the Dataset constructor. This PR fixes that and sets the implicit top drag coefficient to the ISOMIP+ value of 2.5e-3, which had been left at the MPAS-Ocean default of 1.0e-3. The drag change alters answers for all ISOMIP+ tasks.

These tasks are the reference for the Polaris ISOMIP+ tasks in E3SM-Project/polaris#826.

Checklist

  • Document (in a comment titled Testing in this PR) any testing that was used to verify the changes

🤖 Generated with Claude Code

xylar and others added 2 commits September 28, 2026 08:34
Newer xarray no longer accepts a Dataset as the data_vars argument of
the Dataset constructor.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
The ISOMIP+ protocol prescribes a top drag coefficient of 2.5e-3, but
implicit top drag used the MPAS-Ocean default of 1.0e-3.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
@xylar

xylar commented Sep 28, 2026 •

Copy link
Copy Markdown
Collaborator Author

Testing

I ran ocean/isomip_plus/planar/2km/z-star/ Ocean0, Ocean1, Ocean2 and thin_film_Ocean0 on Chrysalis with gnu and OpenMPI, followed by one month of each simulation step. All passed:

/lcrc/group/e3sm/ac.xylar/polaris_1.1/chrysalis/test_20260928/compass_stage3_2km

Without the fix, Ocean0, Ocean1 and Ocean2 failed in initial_state with the xarray TypeError.


Posted by Claude Code on @xylar's behalf. The testing, analysis and wording above are AI-authored; please check them accordingly.

@xylar
xylar marked this pull request as ready for review September 28, 2026 20:27
@xylar xylar self-assigned this Sep 28, 2026
@xylar xylar added bug Something isn't working ocean labels Sep 28, 2026
@xylar
xylar merged commit a876e74 into MPAS-Dev:main Sep 28, 2026
5 checks passed
@xylar
xylar deleted the fix-isomip-plus branch September 28, 2026 20:27
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

bug Something isn't working ocean

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant