Changeset 23090

Timestamp:

08/15/18 07:41:20 (7 years ago)

Author:

tsantos

Message:

BUG: old friction model in examples/shakti

Location:

issm/trunk-jpl/examples/shakti

Files:

• moulin.par (modified) (1 diff)
• runme.m (modified) (1 diff)

issm/trunk-jpl/examples/shakti/moulin.par

@@ -25 +25 @@
 
 % Friction law and coefficient
-    md.friction=frictionsommers(md.friction);
-    md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
+md.friction=frictionshakti(md.friction);
+md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
 
 
 %Numerical parameters
-md.stressbalance.viscosity_overshoot=0.0;
+%md.stressbalance.viscosity_overshoot=0.0;
 md.masstransport.stabilization=1.;
 md.thermal.stabilization=1.;

issm/trunk-jpl/examples/shakti/runme.m

@@ -20 +20 @@
         % HYDROLOGY SPECIFIC PARAMETERIZATION:
         % Change hydrology class to Sommers' SHaKTI model
-        md.hydrology=hydrologysommers();
-
+        md.hydrology=hydrologyshakti();
+        
         % Define initial water head such that water pressure is 50% of ice overburden pressure
         md.hydrology.head = 0.5*md.materials.rho_ice/md.materials.rho_freshwater*md.geometry.thickness + md.geometry.base;