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

Last change on this file since 5966 was 5966, checked in by Mathieu Morlighem, 15 years ago

set np to 8 for EISMINT and ISMIP
  • Property svn:executable set to *
File size: 861 bytes
Line 
1%Ok, start defining model parameters here
2
3disp(' creating thickness');
4md.surface=-md.x*tan(0.5*pi/180);
5md.bed=md.surface-1000+500*sin(md.x*2*pi/max(md.x));
6md.thickness=md.surface-md.bed;
7md.firn_layer=0*ones(md.numberofgrids,1);
8
9disp(' creating drag');
10md.drag_type=2; %0 none 1 plastic 2 viscous
11md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
12%Take care of iceshelves: no basal drag
13pos=find(md.elementoniceshelf);
14md.drag_coefficient(md.elements(pos,:))=0;
15md.drag_p=ones(md.numberofelements,1);
16md.drag_q=ones(md.numberofelements,1);
17
18disp(' creating flow law paramter');
19md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
20md.rheology_n=3*ones(md.numberofelements,1);
21
22disp(' boundary conditions for diagnostic model');
23%Create grid on boundary fist (because we cannot use mesh)
24md=SetIceSheetBC(md);
25md=SetParallel(md,8);
Note: See TracBrowser for help on using the repository browser.