source: issm/trunk-jpl/test/Par/ISMIPC.par@ 24863

Last change on this file since 24863 was 24863, checked in by Mathieu Morlighem, 5 years ago

CHG: renaming groundedice_levelset -> ocean_levelset (negative if ocean present, positive outside
  • Property svn:executable set to *
File size: 1.1 KB
Line 
1%Ok, start defining model parameters here
2
3disp(' creating thickness');
4md.geometry.surface=2000.-md.mesh.x*tan(0.1*pi/180.); %to have z>0
5md.geometry.base=md.geometry.surface-1000.;
6md.geometry.thickness=md.geometry.surface-md.geometry.base;
7
8disp(' creating drag');
9%md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x/2.)).*sin(md.mesh.y*2.*pi/max(md.mesh.x/2.)))./(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base)));
10md.friction.coefficient=sqrt(md.constants.yts.*(1000.+1000.*sin(md.mesh.x*2.*pi/max(md.mesh.x)).*sin(md.mesh.y*2.*pi/max(md.mesh.x))));
11md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
12md.friction.p=ones(md.mesh.numberofelements,1);
13md.friction.q=zeros(md.mesh.numberofelements,1);
14
15disp(' creating flow law parameter');
16md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
17md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
18
19disp(' boundary conditions for stressbalance model:');
20%Create node on boundary first (because we can not use mesh)
21md=SetIceSheetBC(md);
Note: See TracBrowser for help on using the repository browser.