source: issm/trunk/test/Validation/PattynStokes/test2_icesheet/Square.par@ 4875

Last change on this file since 4875 was 4875, checked in by seroussi, 15 years ago

Validation tests PattynStokes

  • Property svn:executable set to *
File size: 1.0 KB
Line 
1%Ok, start defining model parameters here
2
3disp(' creating thickness');
4hmin=300;
5hmax=1000;
6ymin=min(md.y);
7ymax=max(md.y);
8md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
9md.firn_layer=10*ones(md.numberofgrids,1);
10md.bed=-md.rho_ice/md.rho_water*md.thickness+10;
11md.surface=md.bed+md.thickness;
12
13disp(' creating drag');
14md.drag_type=2; %0 none 1 plastic 2 viscous
15md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
16%Take care of iceshelves: no basal drag
17pos=find(md.elementoniceshelf);
18md.drag_coefficient(md.elements(pos,:))=0;
19md.drag_p=ones(md.numberofelements,1);
20md.drag_q=ones(md.numberofelements,1);
21
22disp(' creating temperature');
23md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
24
25disp(' creating flow law paramter');
26md.rheology_B=paterson(md.observed_temperature);
27md.rheology_n=3*ones(md.numberofelements,1);
28
29%Deal with boundary conditions:
30disp(' boundary conditions for diagnostic model: ');
31md=SetMarineIceSheetBC(md,'Front.exp');
32
33%Parallel options
34md.np=3;
35md.time=50;
36md.waitonlock=10;
Note: See TracBrowser for help on using the repository browser.