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

Last change on this file since 4953 was 4953, checked in by Mathieu Morlighem, 15 years ago

updated par files

  • Property svn:executable set to *
File size: 1.0 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,0);
6md.vy_obs =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.x,md.y,0);
7md.surface =InterpFromMeshToMesh2d(index,x,y,surface,md.x,md.y,0);
8md.thickness=InterpFromMeshToMesh2d(index,x,y,thickness,md.x,md.y,0);
9md.bed=md.surface-md.thickness;
10clear surface thickness vx_obs vy_obs x y index;
11
12%Materials
13md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
14md.rheology_B=paterson(md.observed_temperature);
15md.rheology_n=3*ones(md.numberofelements,1);
16
17%Friction
18pos=find(md.elementoniceshelf);
19md.drag_type=2;
20md.drag_coefficient=50*ones(md.numberofgrids,1);
21md.drag_coefficient(md.elements(pos,:))=0;
22md.drag_p=ones(md.numberofelements,1);
23md.drag_q=ones(md.numberofelements,1);
24
25%Numerical parameters
26md.viscosity_overshoot=0.3;
27md.artificial_diffusivity=1;
28md.verbose=0;
29md.waitonlock=30;
30md.np=3;
31md.eps_res=0.05;
32md.eps_rel=0.005;
33md.eps_abs=NaN;
34
35%Boundary conditions:
36md=SetMarineIceSheetBC(md);
Note: See TracBrowser for help on using the repository browser.