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 issues in refactor #307

Merged
merged 3 commits into from
Sep 4, 2024

Conversation

jedwards4b
Copy link
Contributor

Description of changes

Fix issues in restart of multiple ice sheets.

Specific notes

Contributors other than yourself, if any:

CDEPS Issues Fixed (include github issue #):

Are there dependencies on other component PRs (if so list):

Are changes expected to change answers (bfb, different to roundoff, more substantial):

Any User Interface Changes (namelist or namelist defaults changes):

Testing performed (e.g. aux_cdeps, CESM prealpha, etc):

Hashes used for testing:

@jedwards4b jedwards4b self-assigned this Sep 4, 2024
@jedwards4b jedwards4b marked this pull request as draft September 4, 2024 04:09
@mvertens
Copy link
Collaborator

mvertens commented Sep 4, 2024

@jedwards4b @billsacks @Katetc - I'd really like the fix to the restarts to go in its own PR - particularly since I want to make sure that the results are bfb with the restart change. Can this please be done and have Kate's PR come in separately. I think that would be much cleaner and also important for me to use in my current NorESM simulations.

Copy link
Collaborator

@mvertens mvertens left a comment

Choose a reason for hiding this comment

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

I don't see how this worked for a multiple ice sheet test(see my comment in the code review). Could you please summarize the test that passed?

allocate(gindex(lsize))
call ESMF_DistGridGet(distGrid, localDe=0, seqIndexList=gindex, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return

call pio_initdecomp(pio_subsystem, pio_double, (/nx_global(ns),ny_global(ns)/), gindex, pio_iodesc)
call pio_write_darray(pioid, varid, pio_iodesc, Fgrg_rofi(ns)%ptr, rcode, fillval=shr_const_spval)
Copy link
Collaborator

Choose a reason for hiding this comment

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

This is incorrect in that varid should be an array and not a scalar. How does this ever work for multiple ice sheets?

Copy link
Collaborator

Choose a reason for hiding this comment

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

This refers to line 510 above which was not modified and I think should be.

Copy link
Collaborator

@mvertens mvertens left a comment

Choose a reason for hiding this comment

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

Please see my comment regarding making varid an array in the restart write.

@jedwards4b jedwards4b marked this pull request as ready for review September 4, 2024 13:56
Copy link
Collaborator

@mvertens mvertens left a comment

Choose a reason for hiding this comment

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

Thanks!

@jedwards4b jedwards4b merged commit 8b7314b into ESCOMP:main Sep 4, 2024
1 check passed
@jedwards4b jedwards4b deleted the dglc_fix_for_multiple_icesheets branch September 4, 2024 15:23
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.

3 participants