Page 1 of 1

Glacier width in FreeSurfaceSolver

Posted: 04 Jun 2015, 12:08
by vanGeffen

I have a question: is there an easy way to account for the glacier width W(x) in the FreeSurfaceSolver when running a flow-line model? I want to determine the free surface following the equation for the ice thickness as depicted in the attachment.


ChangeIceThickness.png (24.96 KiB) Viewed 2745 times
with H the ice thickness, u the depth-averaged horizontal velocity, W the glacier width and b the surface mass balance

Re: Glacier width in FreeSurfaceSolver

Posted: 09 Jun 2015, 15:58
by joeatodd
Hi Sharon,

It's possible to do this as an accumulation term in the freesurface equation, but I think that it's more accurate to directly account for the addition of mass in the stokes equations.

I wrote some code to do this in a flowline case by adding a "Pressure Load" following Peter's advice in this thread:


The plan is for this code to make its way into the elmer/ice repository before too long. If you like I can pass it on to you for the time being, but it comes without any warranty ;)