Test on an idealized geometry 
Mass conservation
To launch this test, use runme

The velocity is fixed everywhere with:
u=0, v=-400m/yr
The thickness is initialised with 500m everywhere (in hydrostatic equilibrium)

The upper boundary (y=20000) is a dirichlet for prognostic model. 
Tis thickness has to evolve with time. So you have to add these two lines in the icetransient2d.m just before running the iceprognostic_core

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
pos=find(m_p.ys~=0);
m_p.ys(pos)=500+500/5*sin(2*pi*time/(200*md.yts))*ones(size(pos,1),1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
The propagation of the signal is attenuated with time, a little bit less if there is no artificial diffusivity.
One can usethe CrossLine to see precisely the thickness in the middle of the iceshelf.
