1 | %Start defining model parameters here
|
---|
2 |
|
---|
3 | %Geometry and observation
|
---|
4 | load('./../Data/Pig.data','-mat');
|
---|
5 | md.vx_obs =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y);
|
---|
6 | md.vy_obs =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y);
|
---|
7 | md.surface =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y);
|
---|
8 | md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y);
|
---|
9 | md.bed=md.surface-md.thickness;
|
---|
10 | clear surface thickness vx_obs vy_obs x y index;
|
---|
11 | md.vx=md.vx_obs;
|
---|
12 | md.vy=md.vy_obs;
|
---|
13 | md.vz=zeros(md.numberofgrids,1);
|
---|
14 | md.pressure=zeros(md.numberofgrids,1);
|
---|
15 |
|
---|
16 | %Materials
|
---|
17 | md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
|
---|
18 | md.rheology_B=paterson(md.observed_temperature);
|
---|
19 | md.rheology_n=3*ones(md.numberofelements,1);
|
---|
20 | md.temperature=md.observed_temperature;
|
---|
21 |
|
---|
22 | %Friction
|
---|
23 | pos=find(md.elementoniceshelf);
|
---|
24 | md.drag_type=2;
|
---|
25 | md.drag_coefficient=50*ones(md.numberofgrids,1);
|
---|
26 | md.drag_coefficient(md.elements(pos,:))=0;
|
---|
27 | md.drag_p=ones(md.numberofelements,1);
|
---|
28 | md.drag_q=ones(md.numberofelements,1);
|
---|
29 |
|
---|
30 | %Numerical parameters
|
---|
31 | md.viscosity_overshoot=0.3;
|
---|
32 | md.artificial_diffusivity=1;
|
---|
33 | md.verbose=0;
|
---|
34 | md.waitonlock=30;
|
---|
35 | md.np=3;
|
---|
36 | md.dt=1;
|
---|
37 | md.ndt=2;
|
---|
38 | md.eps_res=0.05;
|
---|
39 | md.eps_rel=1;
|
---|
40 | md.eps_abs=NaN;
|
---|
41 |
|
---|
42 | %Boundary conditions:
|
---|
43 | md=SetMarineIceSheetBC(md);
|
---|
44 |
|
---|
45 | %Change name so that no test have the same name
|
---|
46 | A=dbstack;
|
---|
47 | if (length(A)>2), md.name=A(3).file(1:end-2); end
|
---|