md=model;
md=mesh(md,'DomainOutline.exp',50000);
md=geography(md,'Shelf.exp','');
md=parameterize(md,'Square.par');
md.acceleration=1;
md=solve(md,'diagnostic','ice');
