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