Context Navigation

ISMIPC.par@ 17590

Last change on this file since 17590 was 17590, checked in by Mathieu Morlighem, 11 years ago
CHG: WARNING: md.geometry.bed is now md.geometry.base and bathymetry is bed
Property svn:executable set to ``*
File size: 1.2 KB

Rev	Line
[5049]	1	%Ok, start defining model parameters here
	2
	3	disp(' creating thickness');
[14134]	4	md.geometry.surface=2000.-md.mesh.xtan(0.1pi/180.); %to have z>0
[17590]	5	md.geometry.base=md.geometry.surface-1000.;
	6	md.geometry.thickness=md.geometry.surface-md.geometry.base;
[5049]	7
	8	disp(' creating drag');
[17590]	9	%md.friction.coefficient=sqrt(md.constants.yts.(1000.+1000.sin(md.mesh.x2.pi/max(md.mesh.x/2.)).sin(md.mesh.y2.pi/max(md.mesh.x/2.)))./(md.constants.g(md.materials.rho_icemd.geometry.thickness+md.materials.rho_watermd.geometry.base)));
[14134]	10	md.friction.coefficient=sqrt(md.constants.yts.(1000.+1000.sin(md.mesh.x2.pi/max(md.mesh.x)).sin(md.mesh.y2.*pi/max(md.mesh.x))));
[15988]	11	md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
[9725]	12	md.friction.p=ones(md.mesh.numberofelements,1);
	13	md.friction.q=zeros(md.mesh.numberofelements,1);
[5049]	14
[14134]	15	disp(' creating flow law parameter');
[9725]	16	md.materials.rheology_B=6.806710^7ones(md.mesh.numberofvertices,1);
[14134]	17	md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
[5049]	18
[16167]	19	%Damage
	20	md.damage.D=zeros(md.mesh.numberofvertices,1);
	21
[15771]	22	disp(' boundary conditions for stressbalance model:');
[14134]	23	%Create node on boundary first (because we can not use mesh)
[5049]	24	md=SetIceSheetBC(md);

Note: See TracBrowser for help on using the repository browser.