Last change
on this file since 9702 was 9702, checked in by Mathieu Morlighem, 14 years ago |
Added autofiff.m class and moved some other fields to settings
|
-
Property svn:executable
set to
*
|
File size:
1.2 KB
|
Rev | Line | |
---|
[8011] | 1 | %Start defining model parameters here
|
---|
| 2 |
|
---|
| 3 | %Geometry
|
---|
| 4 | hmin=300;
|
---|
| 5 | hmax=1000;
|
---|
| 6 | ymin=min(md.y);
|
---|
| 7 | ymax=max(md.y);
|
---|
[9691] | 8 | md.geometry.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
|
---|
| 9 | md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness;
|
---|
| 10 | md.geometry.surface=md.geometry.bed+md.geometry.thickness;
|
---|
[8011] | 11 |
|
---|
| 12 | %Materials
|
---|
[8307] | 13 | md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
|
---|
[9636] | 14 | md.materials.rheology_B=paterson(md.observed_temperature);
|
---|
| 15 | md.materials.rheology_n=3*ones(md.numberofelements,1);
|
---|
[9684] | 16 | md.initialization.temperature=md.observed_temperature;
|
---|
[8011] | 17 |
|
---|
| 18 | %Friction
|
---|
[9641] | 19 | pos=find(md.mask.elementonfloatingice);
|
---|
[8011] | 20 | md.drag_type=2;
|
---|
[8307] | 21 | md.drag_coefficient=20*ones(md.numberofnodes,1);
|
---|
[8011] | 22 | md.drag_coefficient(md.elements(pos,:))=0;
|
---|
| 23 | md.drag_p=ones(md.numberofelements,1);
|
---|
| 24 | md.drag_q=ones(md.numberofelements,1);
|
---|
| 25 |
|
---|
| 26 | %Numerical parameters
|
---|
| 27 | md.viscosity_overshoot=0.3;
|
---|
[9642] | 28 | md.prognostic.stabilization=1;
|
---|
[9632] | 29 | md.thermal.stabilization=1;
|
---|
[9702] | 30 | md.settings.waitonlock=30;
|
---|
[8011] | 31 | md.verbose=verbose(0);
|
---|
| 32 | md.eps_res=0.10;
|
---|
| 33 | md.eps_rel=0.02;
|
---|
| 34 | md.eps_abs=NaN;
|
---|
| 35 | md.dt=1;
|
---|
| 36 | md.ndt=3;
|
---|
| 37 |
|
---|
| 38 | %Boundary conditions:
|
---|
| 39 | md=SetIceShelfBC(md,'SquareFront.exp');
|
---|
| 40 |
|
---|
| 41 | %Change name so that no test have the same name
|
---|
| 42 | A=dbstack;
|
---|
[9625] | 43 | if (length(A)>2), md.miscellaneous.name=A(3).file(1:end-2); end
|
---|
Note:
See
TracBrowser
for help on using the repository browser.