Last change
on this file since 4949 was 4949, checked in by seroussi, 15 years ago |
added some par files
|
-
Property svn:executable
set to
*
|
File size:
813 bytes
|
Rev | Line | |
---|
[4949] | 1 | %Start defining model parameters here
|
---|
| 2 |
|
---|
| 3 | load('Pigdata.mat');
|
---|
| 4 | md.surface=surface;
|
---|
| 5 | md.bed=bed;
|
---|
| 6 | md.thickness=thickness;
|
---|
| 7 | md.vx_obs=vx_obs;
|
---|
| 8 | md.vy_obs=vy_obs;
|
---|
| 9 |
|
---|
| 10 | %Create drag
|
---|
| 11 | md.drag_type=2;
|
---|
| 12 | md.drag_coefficient=50*ones(md.numberofgrids,1);
|
---|
| 13 | %Take care of iceshelves: no basal drag
|
---|
| 14 | pos=find(md.elementoniceshelf);
|
---|
| 15 | md.drag_coefficient(md.elements(pos,:))=0;
|
---|
| 16 | md.drag_p=ones(md.numberofelements,1);
|
---|
| 17 | md.drag_q=ones(md.numberofelements,1);
|
---|
| 18 | md.viscosity_overshoot=0.3;
|
---|
| 19 |
|
---|
| 20 | %Create temperature and rheology
|
---|
| 21 | md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
|
---|
| 22 | md.rheology_B=paterson(md.observed_temperature);
|
---|
| 23 | md.rheology_n=3*ones(md.numberofelements,1);
|
---|
| 24 |
|
---|
| 25 | %Deal with boundary conditions:
|
---|
| 26 | md=SetMarineIceSheetBC(md);
|
---|
| 27 |
|
---|
| 28 | %Parallel options
|
---|
| 29 | md.np=3;
|
---|
| 30 | md.time=50;
|
---|
| 31 | md.waitonlock=60;
|
---|
| 32 |
|
---|
| 33 | %Misc
|
---|
| 34 | md.verbose=0;
|
---|
| 35 | md.eps_abs=NaN;
|
---|
| 36 | md.eps_rel=0.005;
|
---|
Note:
See
TracBrowser
for help on using the repository browser.