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

Add conservative remapping of observational datasets onto Greenland Ice Sheet mesh #803

Conversation

trhille
Copy link
Collaborator

@trhille trhille commented Mar 28, 2024

This merge adds conservative remapping of observational data sets (thickness, bed topography, velocity, and uncertainties) onto the Greenland Ice Sheet mesh, and handles some of the peculiarities of the other gridded datasets (unit conversions, sign conventions, etc).

Checklist

  • User's Guide has been updated
  • Developer's Guide has been updated
  • Documentation has been built locally and changes look as expected
  • Document (in a comment titled Testing in this PR) any testing that was used to verify the changes

Add conservative remapping directly from GrIS datasets for bed, thickness, velocities, and uncertainties.
Add forward slash after data_path in case path is defined incorrectly, i.e.,
in case user forgets to finish path with forward slash.
Run E3SM_RegridWeightGen on single processor instead of using srun
and multiple nodes. srun was causing cryptic errors, usually claiming to
be out-of-memory errors, but the 3-30km mesh creation worked fine on
one processor. 1-10km case has not yet been tested, and may require the
issue with srun to be resolved.
Also do not create new scrip files for datasets if they already exist in
the data directory.
Take care of large negative dHdt value in a few cells and make sure
dHdt uncertainty is positive
Also add back ESMF_RegridWeightGen on multiple procs.
@trhille trhille added land ice in progress This PR is not ready for review or merging labels Mar 28, 2024
Change data path to default to Perlmutter rather than Chicoma.
Fix file names to reflect updated gridded data set at LCRC.
Make changes required by linter.
@trhille
Copy link
Collaborator Author

trhille commented Mar 28, 2024

@andrewdnolan, as of 62155ab, this runs successfully on one compute node on Perlmutter in ~11.5 minutes. We ultimately want this to follow generally the same workflow as the antarctica/mesh_gen case, so all (or almost all?) of the code that I've added to compass/landice/tests/greenland/mesh.py should be moved to compass/landice/mesh.py.

@trhille trhille marked this pull request as ready for review April 2, 2024 22:17
@trhille trhille removed the in progress This PR is not ready for review or merging label Apr 2, 2024
@andrewdnolan
Copy link
Contributor

Testing

  • Machine: perlmutter
  • Compass version: 1.3.0-alpha.2

Command run:

compass setup -n 45 -w wrkdir -f custom.cfg

where custom.cfg

[mesh]

levels = 20

min_spac = 6.e3
max_spac = 60.e3

[parallel]
account = m4288

[job]

qos=debug
wall_time = 0:30:00

Results:

landice/greenland/mesh_gen
  * step: mesh
  test execution:      SUCCESS
  test runtime:        10:23
Test Runtimes:
10:23 PASS landice_greenland_mesh_gen
Total runtime 10:27
PASS: All passed successfully!

Add explanation of GrIS data interpolation to users guide. Also
update the path to data sets in the Antarctica users guide.
Update developers guide to explain data interpolation in Greenland
mesh_gen test case.
Copy link
Member

@matthewhoffman matthewhoffman left a comment

Choose a reason for hiding this comment

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

@trhille , great to have this finalized! Looks good to me.

@matthewhoffman matthewhoffman merged commit dc641b6 into MPAS-Dev:main Apr 24, 2024
4 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants