Context Navigation

ISMIPB.par@ 16167

Last change on this file since 16167 was 16167, checked in by Eric.Larour, 12 years ago
CHG: switched off Matdamageice. Damange D and rheology B are now equally treated in the material Matice. Adjusted all control methods and parameter files accordingly.
Property svn:executable set to ``*
File size: 898 bytes

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

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