md=triangle(model(),'Square.exp',80000.);
md=setmask(md,'all','');
md=parameterize(md,'SquareShelf.par');
md=extrude(md,3,1);
%Set flow equation
...
%Solve