md=mesh(model,'../Exp/Square.exp','../Exp/Rifts.exp',50000);
md=meshprocessrifts(md,'../Exp/Square.exp');
md=geography(md,'all','');
md=parameterize(md,'../Par/SquareShelf.par');
md=setelementstype(md,'macayeal','all');
md.cluster=oshostname();
md.solverstring=' -mat_type aijmumps -ksp_type preonly -pc_type lu -mat_mumps_icntl_14 80 ';
md=solversettomumps(md);
error('rifts not working for now')
md=solve(md,'analysis_type',DiagnosticSolutionEnum);
