Skip to content

Run ISSM over multiple input files (multiple ice sheets or glaciers)#3

Merged
agstub merged 17 commits intofeature/agstub/newcouplingfrom
feature/agstub/multi-glac
Feb 24, 2026
Merged

Run ISSM over multiple input files (multiple ice sheets or glaciers)#3
agstub merged 17 commits intofeature/agstub/newcouplingfrom
feature/agstub/multi-glac

Conversation

@agstub
Copy link
Collaborator

@agstub agstub commented Feb 23, 2026

This PR adds the capability to run ISSM over multiple input files (e.g., Greenland + Antarctica) intended for coupling with the GEOS Earth system model. This solution is based on the fact that ISSM simulations are encapsulated in FemModel objects. The main changes to the workflow in esmfbinders.cpp are as follows:

  1. InitializeISSM: Given an experiment directory EXPDIR, count all of the input (.bin) files and initialize a FemModel object for each one. Return the number of elements and nodes across all models on each process.
  2. GetElementsISSM and GetNodesISSM: Create IDs for the mesh elements and nodes that are unique across all models (input files / glaciers).
  3. On the ESMF/GEOS side: Create a global ESMF_Mesh object that is the union of all ISSM meshes (i.e., Greenland + Antarctica + whatever else is in EXPDIR). This ESMF_Mesh object can be used for regridding imports/exports between the GCM and ISSM. For an example see examples/GEOSInput/main_esmf.f90
  4. RunISSM: Given the time step and SMB forcing from the GCM, run a single time step of each model. Import/export pointers (defined on the global mesh) are split up appropriately across models and processes.

Just FYI: @lcandre2 @gmao-rreichle @weiyuan-jiang

@agstub agstub self-assigned this Feb 23, 2026
@agstub
Copy link
Collaborator Author

agstub commented Feb 24, 2026

works as expected in GEOS
(on discover, see: /discover/nobackup/agstubbl/GEOS/notebooks/plotting.ipynb)
will merge branch, create new module, update gridcomp branch

globalsurf aisflow

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.

1 participant