Skip to content

Adjustments to ISMIP7 streams, namelist, and albany_input.yaml - #999

Merged
matthewhoffman merged 20 commits into
MPAS-Dev:mainfrom
trhille:landice/ismip7_streams
Nov 12, 2026
Merged

matthewhoffman merged 20 commits into
MPAS-Dev:mainfrom
trhille:landice/ismip7_streams

Conversation

@trhille

@trhille trhille commented Sep 29, 2026 •

Copy link
Copy Markdown
Collaborator

This merge finalizes the streams and namelist configurations for ISMIP7.

This addresses issues tracked in #990 and #991.

The GrIS configuration in this PR is production-ready. The AIS configuration will have further adjustments in a follow-on PR.

@trhille trhille mentioned this pull request Sep 29, 2026
2 of 6 tasks
@trhille trhille added this to the ISMIP7 milestone Sep 29, 2026
@trhille trhille changed the title Final adjustments to ISMIP7 streams and namelist Final adjustments to ISMIP7 streams, namelist, and albany_input.yaml Oct 28, 2026
Add necessary melt parameters for ISMIP7 sub-shelf melt parameterization.
Remove vestigial ismip6_melt_params stream.
Remove deprecated config_enable_timeAvgRestarts option from namelist files.
Use shorter restart intervals to reduce wasted compute time.
Fix type ("Cubrature" instead of "Cubature") in yaml file.
@trhille
trhille force-pushed the landice/ismip7_streams branch from e0614d0 to ac49687 Compare October 28, 2026 15:21
* OCX forcing path
* enable 3d TF
* ctrl forcing paths

Note that it looks like the setup function is not equipped to handle
ctrl 3d TF forcing that is broken into separate files, which is the case
because of large file sizes.  This may be why the 3d tf had been
disabled.

Also note that each of the two controls must be set up separately
because the cfg uses direct paths to the forcing files themselves.
After fixing the MOLHO and FCT bugs, it became necessary to increase
CD strainrate mult to 8 to match restore-calving VAF.  However, at that
level, Jakobhavn was exhibiting much too high and possibly borderline
unstable retreat.

I explored adding crevasse water depth to also increase calving.  Amaral
et al. (2020) recommend crevasse water depth=24 m for a GIS-wide
parameter value, based on observed velocity and terminus position maps.
Water depth=24m alone still greatly underestimates the VAF from
restore-calving.  So I then include strainrate multiplier on top of a
prescribed 24 m water depth, and find strain rate multiplier of 4 works
best.

This commit also changes the allowable CFL fraction to 0.75 after
testing shows similar accuracy to using 0.5.
@matthewhoffman matthewhoffman changed the title Final adjustments to ISMIP7 streams, namelist, and albany_input.yaml Adjustments to ISMIP7 streams, namelist, and albany_input.yaml Nov 4, 2026
@matthewhoffman

Copy link
Copy Markdown
Member

I've noted 2 additional issues. I don't think I want to bother addressing them because they are easy to work around, but I want to flag them here.

  1. For the ctrl run, the 3d TF files don't get set up right because the cfg assumes it's a single file, but they actually are broken into multiyear chunks. Workaround is to manually add the symlinks after the run sets up. (I used copilot because there's a lot of them.)
  2. For the ssp585 run, the 3d TF stream gets an incorrect filename_interval. I looked at the code that is meant to automate that calculation, and it's not immediately obvious what the problem is. Workaround: manually adjust the filename_interval after the run gets created.

matthewhoffman and others added 5 commits November 5, 2026 20:07
That was the case in ISMIP6, but for ISMIP7, the ctrl runs are a monthly
climatology, so the same stream settings apply.
The ctrl experiments were special-cased in ismip7_gris and ismip7_ais,
a holdover from the old constant-climate protocol. They used separate
ctrl_*_climatology_path options and assumed a single 3D thermal forcing
file, so the multi-year GrIS 3D TF chunks were never symlinked.

Both ctrl runs now use {forcing_basepath}/{model}_ctrl and the same
forcing-file discovery and intervals as the other experiments. In AIS,
mesh-name filtering and annual TF are kept, and ctrl still skips the
hydrofracture shelf-collapse mask. The ctrl_*_climatology_path options
are removed from the configs and the docs are updated.

Co-authored-by: Copilot <223556219+Copilot@users.noreply.github.com>
The 3D thermal forcing stream used a fixed reference_time of 2000-01-01,
which puts filename boundaries at 2010, 2020, ... and does not match
projection chunks that start in 2015. Set reference_time to the start of
the first chunk, and fail setup if multiple files share a chunk start
year (e.g. two versions of the series in one directory).

Co-authored-by: Copilot <223556219+Copilot@users.noreply.github.com>
Create graph.info.part.<ntasks> alongside graph.info for GrIS and AIS
so job scripts that run MALI directly find the decomposition.

Co-authored-by: Copilot <223556219+Copilot@users.noreply.github.com>
Add custom_job_script and custom_job_env_cmds options for GrIS and AIS.
With custom_job_script = custom, each experiment also gets a
job_script.custom.sh that runs the linked MALI executable directly.
Enable it in ismip7_gris_pm_production.cfg.

Co-authored-by: Copilot <223556219+Copilot@users.noreply.github.com>
Another round of testing indicates this is the optimal OCX.
Also set ismip7_gris_pm_production.cfg to set up OCX plus all CESM runs initially.
And set account to m4274.
@matthewhoffman
matthewhoffman merged commit ceae905 into MPAS-Dev:main Nov 12, 2026
5 checks passed
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.

2 participants