Changeset 17069

Timestamp:

01/08/14 08:46:57 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing tutorials with new damage field

Location:

issm/trunk-jpl/examples

Files:

• Jakobshavn/Jks.par (modified) (1 diff)
• Jakobshavn/runme.m (modified) (2 diffs)
• SquareIceShelf/Square.par (modified) (1 diff)

issm/trunk-jpl/examples/Jakobshavn/Jks.par

@@ -51 +51 @@
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 md.materials.rheology_B=paterson(md.initialization.temperature);
+md.damage.D=zeros(md.mesh.numberofvertices,1);
 %Reduce viscosity along the shear margins
 weakb=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,'WeakB.exp','node',2);

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

@@ -1 @@
-steps=[1];
+steps=[2:3];
 
 if any(steps==1)
@@ -52 +52 @@
 
         %Cost functions
-        md.inversion.cost_functions=[101*ones(md.inversion.nsteps,1) 103*ones(md.inversion.nsteps,1)];
+        md.inversion.cost_functions=[101 103];
         md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2);
         md.inversion.cost_functions_coefficients(:,1)=40;

issm/trunk-jpl/examples/SquareIceShelf/Square.par

@@ -25 +25 @@
 md.materials.rheology_B=paterson((273-20)*ones(md.mesh.numberofvertices,1));
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.damage.D=zeros(md.mesh.numberofvertices,1);
 
 disp('      creating boundary conditions');