Last change
on this file since 9610 was 9610, checked in by Mathieu Morlighem, 14 years ago |
Added friction object
|
-
Property svn:executable
set to
*
|
File size:
769 bytes
|
Line | |
---|
1 | %Ok, start defining model parameters here
|
---|
2 |
|
---|
3 | disp(' creating thickness');
|
---|
4 | md.surface=-md.x*tan(0.5*pi/180);
|
---|
5 | md.bed=md.surface-1000+500*sin(md.x*2*pi/max(md.x));
|
---|
6 | md.thickness=md.surface-md.bed;
|
---|
7 |
|
---|
8 | disp(' creating drag');
|
---|
9 | md.friction.coefficient=200*ones(md.numberofnodes,1); %q=1.
|
---|
10 | %Take care of iceshelves: no basal drag
|
---|
11 | pos=find(md.elementoniceshelf);
|
---|
12 | md.friction.coefficient(md.elements(pos,:))=0;
|
---|
13 | md.friction.p=ones(md.numberofelements,1);
|
---|
14 | md.friction.q=ones(md.numberofelements,1);
|
---|
15 |
|
---|
16 | disp(' creating flow law paramter');
|
---|
17 | md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
|
---|
18 | md.rheology_n=3*ones(md.numberofelements,1);
|
---|
19 |
|
---|
20 | disp(' boundary conditions for diagnostic model');
|
---|
21 | %Create node on boundary fist (because we cannot use mesh)
|
---|
22 | md=SetIceSheetBC(md);
|
---|
Note:
See
TracBrowser
for help on using the repository browser.