source: issm/trunk/test/Par/ISMIPB.par@ 9725

Last change on this file since 9725 was 9725, checked in by Mathieu Morlighem, 14 years ago

Some more objects in mesh
  • Property svn:executable set to *
File size: 882 bytes
RevLine 
[5049]1%Ok, start defining model parameters here
2
3disp(' creating thickness');
[9691]4md.geometry.surface=-md.x*tan(0.5*pi/180);
5md.geometry.bed=md.geometry.surface-1000+500*sin(md.x*2*pi/max(md.x));
6md.geometry.thickness=md.geometry.surface-md.geometry.bed;
[5049]7
8disp(' creating drag');
[9725]9md.friction.coefficient=200*ones(md.mesh.numberofvertices,1); %q=1.
[5049]10%Take care of iceshelves: no basal drag
[9641]11pos=find(md.mask.elementonfloatingice);
[9610]12md.friction.coefficient(md.elements(pos,:))=0;
[9725]13md.friction.p=ones(md.mesh.numberofelements,1);
14md.friction.q=ones(md.mesh.numberofelements,1);
[5049]15
16disp(' creating flow law paramter');
[9725]17md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
18md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
[5049]19
20disp(' boundary conditions for diagnostic model');
[8308]21%Create node on boundary fist (because we cannot use mesh)
[5049]22md=SetIceSheetBC(md);
Note: See TracBrowser for help on using the repository browser.