source: issm/trunk/test/Miscellaneous/GJM_test1/SquareShelf.par@ 9734

Last change on this file since 9734 was 9734, checked in by seroussi, 14 years ago

finished mesh

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