md=mesh(model,'../Exp/Square.exp',200000); md=geography(md,'all',''); md=parameterize(md,'../Par/SquareShelf.par'); md=extrude(md,3,1); md=setelementstype(md,'pattyn','all'); %control parameters md.control_analysis=1; md.control_type='rheology_B'; md.cm_min=10^6; md.cm_max=2*10^9; md.nsteps=2; md.fit=0*ones(md.nsteps,1); md.weights=ones(md.numberofgrids,1); md.optscal=10^7*ones(md.nsteps,1); md.maxiter=2*ones(md.nsteps,1); md.cm_jump=0.3*ones(md.nsteps,1); md.vx_obs=md.vx; md.vy_obs=md.vy; md.cluster=oshostname(); md=solve(md,'analysis_type',DiagnosticSolutionEnum);