source: issm/trunk/test/Validation/ISMIP/TestD/Square.par@ 904

Last change on this file since 904 was 904, checked in by seroussi, 16 years ago

simplified Square.par

  • Property svn:executable set to *
File size: 1.2 KB
RevLine 
[16]1%Ok, start defining model parameters here
2
[904]3disp(' creating thickness');
4md.surface=2000-md.x*tan(0.1*pi/180); %to have z>0
5md.bed=md.surface-1000;
6md.thickness=md.surface-md.bed;
7md.firn_layer=0*ones(md.numberofgrids,1);
[16]8
[904]9disp(' creating drag');
10md.drag_type=2; %0 none 1 plastic 2 viscous
11md.drag=sqrt(md.yts.*(1000+1000*sin(md.x*2*pi/max(md.x)))./(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)));
12%Take care of iceshelves: no basal drag
13pos=find(md.elementoniceshelf);
14md.drag(md.elements(pos,:))=0;
15md.p=ones(md.numberofelements,1);
16md.q=ones(md.numberofelements,1);
[16]17
[904]18disp(' creating flow law paramter');
19md.B=6.8067*10^7*ones(md.numberofgrids,1);
20md.n=3*ones(md.numberofelements,1);
[16]21
[904]22disp(' boundary conditions for diagnostic model: ');
23%Create grid on boundary fist (because wi can not use mesh)
24md.gridonboundary=zeros(md.numberofgrids,1);
25pos=find(md.x==0 | md.y==0 | md.x==max(md.x) | md.y==max(md.y));
26md.gridonboundary(pos)=1;
27md.gridondirichlet_diag=zeros(md.numberofgrids,1);
28pos=find(md.gridonboundary);md.gridondirichlet_diag(pos)=1;
29md.dirichletvalues_diag=zeros(md.numberofgrids,2);
30md.segmentonneumann_diag=zeros(0,3);
31md.gridondirichlet_thermal=ones(md.numberofgrids,1); %surface temperature
Note: See TracBrowser for help on using the repository browser.