source: issm/trunk-jpl/test/Par/ISMIPA.par@ 17769

Last change on this file since 17769 was 17769, checked in by cborstad, 11 years ago

CHG: reverting previous change, isdamage=0 is specified by default in the damage class, problem was a typo elsewhere...
  • Property svn:executable set to *
File size: 881 bytes
Line 
1%Ok, start defining model parameters here
2
3disp(' creating thickness');
4md.geometry.surface=-md.mesh.x*tan(0.5*pi/180.);
5md.geometry.base=md.geometry.surface-1000.+500.*sin(md.mesh.x*2.*pi/max(md.mesh.x)).*sin(md.mesh.y*2.*pi/max(md.mesh.x));
6md.geometry.thickness=md.geometry.surface-md.geometry.base;
7
8disp(' creating drag');
9md.friction.coefficient=200.*ones(md.mesh.numberofvertices,1);
10md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
11md.friction.p=ones(md.mesh.numberofelements,1);
12md.friction.q=ones(md.mesh.numberofelements,1);
13
14disp(' creating flow law parameter');
15md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
16md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
17
18disp(' boundary conditions for stressbalance model');
19%Create node on boundary first (because we cannot use mesh)
20md=SetIceSheetBC(md);
Note: See TracBrowser for help on using the repository browser.