search

OGGM / oggm

Open Global Glacier Model
http://oggm.org
BSD 3-Clause "New" or "Revised" License
219 stars 105 forks source link

How to deal with mu=0 after calving inversion #1000

Open fmaussion opened 4 years ago

fmaussion commented 4 years ago

In @bearecinos 's paper, we noticed that certain glaciers don't have a solution to the equation: "ice deformation flux = calving law flux". In these cases, we set the calving flux to the maximum possible one while conserving mass, i.e. no melt on the glacier and all precipitation is calving under equilibrium conditions.

This was kind of OK for the inversion (it gave a nice "physically based" upper bound on calving), but running a simulation with these glaciers makes no sense since mu equal 0 and no melt can occur.

How to deal with these glaciers?

One solution would be to reduce k_calving for these glaciers until mu reaches a minimum threshold (which is not 0). Not great, but....

fmaussion commented 4 years ago

At very least we should maybe error when trying to run a simulation on a glacier than cannot melt...