source: issm/trunk/test/Miscellaneous/GJM_test1/SquareShelf.par@ 9734

%Start defining model parameters here

%Geometry
hmin=300;
hmax=1000;
ymin=min(md.mesh.y);
ymax=max(md.mesh.y);
md.geometry.thickness=hmax+(hmin-hmax)*(md.mesh.y-ymin)/(ymax-ymin);
md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness;
md.geometry.surface=md.geometry.bed+md.geometry.thickness;

%Materials
md.observed_temperature=(273-20)*ones(md.mesh.numberofvertices,1);
md.materials.rheology_B=paterson(md.observed_temperature);
md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
md.initialization.temperature=md.observed_temperature;

%Friction
pos=find(md.mask.elementonfloatingice);
md.drag_type=2;
md.drag_coefficient=20*ones(md.mesh.numberofvertices,1);
md.drag_coefficient(md.mesh.elements(pos,:))=0;
md.drag_p=ones(md.mesh.numberofelements,1);
md.drag_q=ones(md.mesh.numberofelements,1);

%Numerical parameters
md.viscosity_overshoot=0.3;
md.prognostic.stabilization=1;
md.thermal.stabilization=1;
md.settings.waitonlock=30;
md.verbose=verbose(0);
md.eps_res=0.10;
md.eps_rel=0.02;
md.eps_abs=NaN;
md.dt=1;
md.ndt=3;

%Boundary conditions:
md=SetIceShelfBC(md,'SquareFront.exp');

%Change name so that no test have the same name
A=dbstack;
if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end