source: issm/trunk-jpl/test/Par/ISMIPC.par@ 13137

Last change on this file since 13137 was 13137, checked in by Mathieu Morlighem, 13 years ago

CHG: clean up as no damage is needed
  • Property svn:executable set to *
File size: 1.2 KB
RevLine 
[5049]1%Ok, start defining model parameters here
2
3disp(' creating thickness');
[9734]4md.geometry.surface=2000-md.mesh.x*tan(0.1*pi/180); %to have z>0
[9691]5md.geometry.bed=md.geometry.surface-1000;
6md.geometry.thickness=md.geometry.surface-md.geometry.bed;
[5049]7
8disp(' creating drag');
[9734]9%md.friction.coefficient=sqrt(md.constants.yts.*(1000+1000*sin(md.mesh.x*2*pi/max(md.mesh.x/2)).*sin(md.mesh.y*2*pi/max(md.mesh.x/2)))./(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.bed)));
10md.friction.coefficient=sqrt(md.constants.yts.*(1000+1000*sin(md.mesh.x*2*pi/max(md.mesh.x)).*sin(md.mesh.y*2*pi/max(md.mesh.x))));
[5049]11%Take care of iceshelves: no basal drag
[9641]12pos=find(md.mask.elementonfloatingice);
[9733]13md.friction.coefficient(md.mesh.elements(pos,:))=0;
[9725]14md.friction.p=ones(md.mesh.numberofelements,1);
15md.friction.q=zeros(md.mesh.numberofelements,1);
[5049]16
17disp(' creating flow law paramter');
[9725]18md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
19md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
[5049]20
21disp(' boundary conditions for diagnostic model: ');
[8308]22%Create node on boundary fist (because wi can not use mesh)
[5049]23md=SetIceSheetBC(md);
Note: See TracBrowser for help on using the repository browser.