Skip to content

Interpolate subglacial hydrology forcing fields in time #163

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.

Sign up for GitHub

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

Open
matthewhoffman wants to merge 8 commits into
base: develop
Choose a base branch
from
Open

Interpolate subglacial hydrology forcing fields in time #163

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
dd66ee6
Use hydro-specific domain masks
matthewhoffman Feb 2, 2026
0fe5852
Remove call and reference to li_mask module
matthewhoffman Feb 2, 2026
9089670
Add vars and routines for sgh forcing interp
matthewhoffman Feb 3, 2026
1f35ec7
Activate forcing interpolation
matthewhoffman Feb 4, 2026
0def959
Switch all model calculations to use interpolated-in-time forcing
matthewhoffman Feb 4, 2026
5ff0f08
Move mask calculation into forcing interpolation routine
matthewhoffman Feb 4, 2026
cc2c6ae
Fix issues with pools
matthewhoffman Feb 4, 2026
5263cb7
Fix compile errors
matthewhoffman Feb 4, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
32 changes: 31 additions & 1 deletion components/mpas-albany-landice/src/Registry_subglacial_hydro.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -167,8 +167,20 @@
<!-- state vars -->
<var name="iceThicknessHydro" type="real" dimensions="nCells Time" units="m"
description="ice thickness used by the hydrology model. Same as 'thickness' but with potential differences along domain boundaries that inhibit the formation of local hydropotential minima on boundaries." />
<var name="iceThicknessHydroPrev" type="real" dimensions="nCells Time" units="m"
description="ice thickness used by the hydrology model from previous master time step" />
<var name="bedTopographyPrev" type="real" dimensions="nCells Time" units="m"
description="bed topography used by the hydrology model from previous master time step" />
<var name="basalSpeedPrev" type="real" dimensions="nCells Time" units="m s^-1"
description="basal speed used by the hydrology model from previous master time step" />
<var name="hydroDomainCell" type="integer" dimensions="nCells Time"
description="logical mask of cells where subglacial hydrology model is active" />
<var name="hydroDomainEdge" type="integer" dimensions="nCells Time"
description="logical mask of edges where subglacial hydrology model is active" />
<var name="flowAHydro" type="real" dimensions="nCells Time" units="s^-1 Pa^-n" default_value="3.1709792e-24"
description="flow parameter A value used by subglacial hydrology model" />
<var name="flowAHydroPrev" type="real" dimensions="nCells Time" units="s^-1 Pa^-n" default_value="3.1709792e-24"
description="flow parameter A value used by subglacial hydrology model from previous master time step" />
<var name="waterThickness" type="real" dimensions="nCells Time" units="m" default_value="0.01"
description="water layer thickness in subglacial hydrology system" />
<var name="waterThicknessOld" type="real" dimensions="nCells Time" units="m"
Expand All @@ -190,8 +202,12 @@
<!-- inputs -->
<var name="basalMeltInput" type="real" dimensions="nCells Time" units="kg m^-2 s^-1"
description="basal meltwater input to subglacial hydrology system" />
<var name="basalMeltInputPrev" type="real" dimensions="nCells Time" units="kg m^-2 s^-1"
description="basal meltwater input to subglacial hydrology system from previous master time step" />
<var name="externalWaterInput" type="real" dimensions="nCells Time" units="kg m^-2 s^-1"
description="external water input to subglacial hydrology system" />
<var name="externalWaterInputPrev" type="real" dimensions="nCells Time" units="kg m^-2 s^-1"
description="external water input to subglacial hydrology system from previous master time step" />
<var name="frictionAngle" type="real" dimensions="nCells Time"
description="subglacial till friction angle" />
<!-- convenience variables -->
Expand Down Expand Up @@ -303,7 +319,21 @@
<var name="channelDiffusivity" type="real" dimensions="nEdges Time" units="m^2 s^-1"
description="diffusivity in channel in subglacial hydrology system" />
</var_struct>

<!-- Hydro Scratch Variables -->
<var_struct name="hydro_scratch" time_levs="1">
<var name="iceThicknessHydroNow" type="real" dimensions="nCells Time" units="m" persistence="scratch"
description="ice thickness used by the hydrology model for current subcycle" />
<var name="bedTopographyNow" type="real" dimensions="nCells Time" units="m" persistence="scratch"
description="bed topography used by the hydrology model for current subcycle" />
<var name="basalSpeedNow" type="real" dimensions="nCells Time" units="m s^-1" persistence="scratch"
description="basal speed used by the hydrology model for current subcycle" />
<var name="flowAHydroNow" type="real" dimensions="nCells Time" units="s^-1 Pa^-n" default_value="3.1709792e-24" persistence="scratch"
description="flow parameter A value used by subglacial hydrology model for current subcycle" />
<var name="basalMeltInputNow" type="real" dimensions="nCells Time" units="kg m^-2 s^-1" persistence="scratch"
description="basal meltwater input to subglacial hydrology system at current subcycle" />
<var name="externalWaterInputNow" type="real" dimensions="nCells Time" units="kg m^-2 s^-1" persistence="scratch"
description="external water input to subglacial hydrology system at current subcycle" />
</var_struct>
<!-- ======================================================================= -->


Loading