Skip to content

Implement thermodynamics for melting of excess ice#13

Draft
rfiorella wants to merge 3 commits intolanl/IM1from
excess-ice-mass
Draft

Implement thermodynamics for melting of excess ice#13
rfiorella wants to merge 3 commits intolanl/IM1from
excess-ice-mass

Conversation

@rfiorella
Copy link
Collaborator

Replace the original subsidence routine in ActiveLayerMod by coupling excess ice to the energy balance equations now.

Major Changes:

  • Excess ice now stored as mass (kg/m2) internally, volumetric fraction for restart I/O only
  • Phase change (melting/freezing) controlled by energy balance in PhaseChange_beta
  • Sequential melting: pore ice melts first, then excess ice with surplus heat
  • Excess ice does NOT refreeze
  • Subsidence tracked cumulatively since 1989 based on actual mass change
  • Geometric melting code in ActiveLayerMod commented out but retained

Replace the original subsidence routine in ActiveLayerMod
by coupling excess ice to the energy balance equations now.

Major Changes:
- Excess ice now stored as mass (kg/m2) internally, volumetric fraction
  for restart I/O only
- Phase change (melting/freezing) controlled by energy balance in
  PhaseChange_beta
- Sequential melting: pore ice melts first, then excess ice with surplus heat
- Excess ice does NOT refreeze
- Subsidence tracked cumulatively since 1989 based on actual mass change
- Geometric melting code in ActiveLayerMod commented out but retained
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