source: issm/trunk/test/Par/Pig.par@ 9625

Last change on this file since 9625 was 9625, checked in by Mathieu Morlighem, 14 years ago

Added miscellaneous class

  • Property svn:executable set to *
File size: 1.3 KB
Line 
1%Start defining model parameters here
2
3%Geometry and observation
4load('../Data/Pig.data','-mat');
5md.vx_obs =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.x,md.y);
6md.vy_obs =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y);
7md.surface =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y);
8md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y);
9md.bed=md.surface-md.thickness;
10clear surface thickness vx_obs vy_obs x y index;
11md.vx=md.vx_obs;
12md.vy=md.vy_obs;
13md.vz=zeros(md.numberofnodes,1);
14md.pressure=zeros(md.numberofnodes,1);
15
16%Materials
17md.temperature=(273-20)*ones(md.numberofnodes,1);
18md.rheology_B=paterson(md.temperature);
19md.rheology_n=3*ones(md.numberofelements,1);
20md.temperature=md.temperature;
21
22%Friction
23pos=find(md.elementoniceshelf);
24md.friction.coefficient=50*ones(md.numberofnodes,1);
25md.friction.coefficient(md.elements(pos,:))=0;
26md.friction.p=ones(md.numberofelements,1);
27md.friction.q=ones(md.numberofelements,1);
28
29%Numerical parameters
30md.viscosity_overshoot=0.3;
31md.artificial_diffusivity=1;
32md.verbose=verbose(0);
33md.waitonlock=30;
34md.dt=1;
35md.ndt=2;
36md.eps_res=0.05;
37md.eps_rel=1;
38md.eps_abs=NaN;
39
40%Boundary conditions:
41md=SetMarineIceSheetBC(md);
42
43%Change name so that no test have the same name
44A=dbstack;
45if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
Note: See TracBrowser for help on using the repository browser.