Context Navigation

ISMIPB.par@ 9636

Last change on this file since 9636 was 9636, checked in by Mathieu Morlighem, 14 years ago
Added material class
Property svn:executable set to ``*
File size: 789 bytes

Rev	Line
[5049]	1	%Ok, start defining model parameters here
	2
	3	disp(' creating thickness');
	4	md.surface=-md.xtan(0.5pi/180);
	5	md.bed=md.surface-1000+500sin(md.x2*pi/max(md.x));
	6	md.thickness=md.surface-md.bed;
	7
	8	disp(' creating drag');
[9610]	9	md.friction.coefficient=200*ones(md.numberofnodes,1); %q=1.
[5049]	10	%Take care of iceshelves: no basal drag
	11	pos=find(md.elementoniceshelf);
[9610]	12	md.friction.coefficient(md.elements(pos,:))=0;
	13	md.friction.p=ones(md.numberofelements,1);
	14	md.friction.q=ones(md.numberofelements,1);
[5049]	15
	16	disp(' creating flow law paramter');
[9636]	17	md.materials.rheology_B=6.806710^7ones(md.numberofnodes,1);
	18	md.materials.rheology_n=3*ones(md.numberofelements,1);
[5049]	19
	20	disp(' boundary conditions for diagnostic model');
[8308]	21	%Create node on boundary fist (because we cannot use mesh)
[5049]	22	md=SetIceSheetBC(md);

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