[4996] | 1 | %Start defining model parameters here
|
---|
| 2 |
|
---|
| 3 | %Geometry
|
---|
| 4 | hmin=300;
|
---|
| 5 | hmax=1000;
|
---|
| 6 | ymin=min(md.y);
|
---|
| 7 | ymax=max(md.y);
|
---|
| 8 | md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
|
---|
[9636] | 9 | md.bed=-md.materials.rho_ice/md.materials.rho_water*md.thickness;
|
---|
| 10 | bed_sheet=-md.materials.rho_ice/md.materials.rho_water*(hmax+(hmin-hmax)*(ymax/2-ymin)/(ymax-ymin));
|
---|
[4996] | 11 | pos=find(md.y<=ymax/2);
|
---|
| 12 | md.bed(pos)=bed_sheet;
|
---|
| 13 | md.surface=md.bed+md.thickness;
|
---|
| 14 |
|
---|
| 15 | %Initial velocity
|
---|
[8939] | 16 | load('../Data/SquareSheetShelf.data','-mat');
|
---|
[5035] | 17 | md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y);
|
---|
| 18 | md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y);
|
---|
[4996] | 19 | clear vx vy x y index;
|
---|
[8305] | 20 | md.vz=zeros(md.numberofnodes,1);
|
---|
| 21 | md.pressure=zeros(md.numberofnodes,1);
|
---|
[4996] | 22 |
|
---|
| 23 | %Materials
|
---|
[9425] | 24 | md.temperature=(273-20)*ones(md.numberofnodes,1);
|
---|
[9636] | 25 | md.materials.rheology_B=paterson(md.temperature);
|
---|
| 26 | md.materials.rheology_n=3*ones(md.numberofelements,1);
|
---|
[4996] | 27 |
|
---|
| 28 | %Accumulation and melting
|
---|
[9607] | 29 | md.surfaceforcings.mass_balance=10*ones(md.numberofnodes,1);
|
---|
[9611] | 30 | md.basalforcings.melting_rate=5*ones(md.numberofnodes,1);
|
---|
[4996] | 31 |
|
---|
| 32 | %Friction
|
---|
[9641] | 33 | pos=find(md.mask.elementonfloatingice);
|
---|
[9610] | 34 | md.friction.coefficient=20*ones(md.numberofnodes,1);
|
---|
| 35 | md.friction.coefficient(md.elements(pos,:))=0;
|
---|
| 36 | md.friction.p=ones(md.numberofelements,1);
|
---|
| 37 | md.friction.q=ones(md.numberofelements,1);
|
---|
[4996] | 38 |
|
---|
| 39 | %Numerical parameters
|
---|
| 40 | md.viscosity_overshoot=0.0;
|
---|
| 41 | md.artificial_diffusivity=1;
|
---|
[9632] | 42 | md.thermal.stabilization=1;
|
---|
[6307] | 43 | md.verbose=verbose(0);
|
---|
[4996] | 44 | md.waitonlock=30;
|
---|
| 45 | md.eps_res=0.05;
|
---|
| 46 | md.eps_rel=0.05;
|
---|
| 47 | md.eps_abs=NaN;
|
---|
[9628] | 48 | md.timestepping.time_step=1;
|
---|
| 49 | md.timestepping.final_time=3;
|
---|
[4996] | 50 |
|
---|
| 51 | %Deal with boundary conditions:
|
---|
[8939] | 52 | md=SetMarineIceSheetBC(md,'../Exp/SquareFront.exp');
|
---|
[4999] | 53 |
|
---|
| 54 | %Change name so that no test have the same name
|
---|
| 55 | A=dbstack;
|
---|
[9625] | 56 | if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
|
---|