Test on an idealized geometry 

Coupled model of ice shelf flow

For this test you can add a dirichlet value for the thickness of the upper boundary condition like this:
%Evolution of the ice shelf
pos=find(md.y==200000); %grids on the upper boundary condition
md.gridondirichlet_prog(pos)=1;
md.dirichletvalues_prog(pos)=500;
md.ndt=5000*md.yts;
md.dt=500*md.yts;

The ice shelf seems to stabilize around 5000 years.
