[4949] | 1 | %Start defining model parameters here
|
---|
| 2 |
|
---|
[4953] | 3 | %Geometry and observation
|
---|
[8824] | 4 | load('../Data/Pig.data','-mat');
|
---|
[9681] | 5 | md.inversion.vx_obs =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y);
|
---|
| 6 | md.inversion.vy_obs =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y);
|
---|
[5035] | 7 | md.surface =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y);
|
---|
| 8 | md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y);
|
---|
[4953] | 9 | md.bed=md.surface-md.thickness;
|
---|
| 10 | clear surface thickness vx_obs vy_obs x y index;
|
---|
[9681] | 11 | md.vx=md.inversion.vx_obs;
|
---|
| 12 | md.vy=md.inversion.vy_obs;
|
---|
[8305] | 13 | md.vz=zeros(md.numberofnodes,1);
|
---|
| 14 | md.pressure=zeros(md.numberofnodes,1);
|
---|
[4949] | 15 |
|
---|
[4953] | 16 | %Materials
|
---|
[9425] | 17 | md.temperature=(273-20)*ones(md.numberofnodes,1);
|
---|
[9636] | 18 | md.materials.rheology_B=paterson(md.temperature);
|
---|
| 19 | md.materials.rheology_n=3*ones(md.numberofelements,1);
|
---|
[9425] | 20 | md.temperature=md.temperature;
|
---|
[4953] | 21 |
|
---|
| 22 | %Friction
|
---|
[9641] | 23 | pos=find(md.mask.elementonfloatingice);
|
---|
[9610] | 24 | md.friction.coefficient=50*ones(md.numberofnodes,1);
|
---|
| 25 | md.friction.coefficient(md.elements(pos,:))=0;
|
---|
| 26 | md.friction.p=ones(md.numberofelements,1);
|
---|
| 27 | md.friction.q=ones(md.numberofelements,1);
|
---|
[4953] | 28 |
|
---|
| 29 | %Numerical parameters
|
---|
[9679] | 30 | md.diagnostic.viscosity_overshoot=0.3;
|
---|
[9642] | 31 | md.prognostic.stabilization=1;
|
---|
[6307] | 32 | md.verbose=verbose(0);
|
---|
[4953] | 33 | md.waitonlock=30;
|
---|
[9628] | 34 | md.timestepping.time_step=1;
|
---|
| 35 | md.timestepping.final_time=2;
|
---|
[9679] | 36 | md.diagnostic.restol=0.05;
|
---|
| 37 | md.diagnostic.reltol=1;
|
---|
[9677] | 38 | md.steadystate.reltol=1;
|
---|
[9679] | 39 | md.diagnostic.abstol=NaN;
|
---|
[4949] | 40 |
|
---|
[4953] | 41 | %Boundary conditions:
|
---|
[4949] | 42 | md=SetMarineIceSheetBC(md);
|
---|
[4999] | 43 |
|
---|
| 44 | %Change name so that no test have the same name
|
---|
| 45 | A=dbstack;
|
---|
[9625] | 46 | if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
|
---|