source: issm/trunk/test/Par/SquareSheetShelf.par@ 8392

%Start defining model parameters here

%Geometry
hmin=300;
hmax=1000;
ymin=min(md.y);
ymax=max(md.y);
md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
md.bed=-md.rho_ice/md.rho_water*md.thickness;
bed_sheet=-md.rho_ice/md.rho_water*(hmax+(hmin-hmax)*(ymax/2-ymin)/(ymax-ymin));
pos=find(md.y<=ymax/2);
md.bed(pos)=bed_sheet;
md.surface=md.bed+md.thickness;

%Initial velocity 
load('./../Data/SquareSheetShelf.data','-mat');
md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y);
md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y);
clear vx vy x y index;
md.vz=zeros(md.numberofnodes,1);
md.pressure=zeros(md.numberofnodes,1);

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

%Accumulation and melting
md.accumulation_rate=10*ones(md.numberofnodes,1);
md.forcings.accumulation_rate=[10*ones(md.numberofnodes,1);1];
md.basal_melting_rate=5*ones(md.numberofnodes,1);

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

%Numerical parameters
md.viscosity_overshoot=0.0;
md.artificial_diffusivity=1;
md.verbose=verbose(0);
md.waitonlock=30;
md.eps_res=0.05;
md.eps_rel=0.05;
md.eps_abs=NaN;
md.dt=1;
md.ndt=3;

%Deal with boundary conditions:
md=SetMarineIceSheetBC(md,'./../Exp/SquareFront.exp');

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