-
Notifications
You must be signed in to change notification settings - Fork 40
Add conservative remapping of observational datasets onto Greenland Ice Sheet mesh #803
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
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.
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.
@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 |
Testing
Command run:
where [mesh]
levels = 20
min_spac = 6.e3
max_spac = 60.e3
[parallel]
account = m4288
[job]
qos=debug
wall_time = 0:30:00 Results:
|
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.
There was a problem hiding this 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.
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
Testing
in this PR) any testing that was used to verify the changes