Index: /issm/trunk-jpl/examples/Jakobshavn/Jks.par
===================================================================
--- /issm/trunk-jpl/examples/Jakobshavn/Jks.par	(revision 17068)
+++ /issm/trunk-jpl/examples/Jakobshavn/Jks.par	(revision 17069)
@@ -51,4 +51,5 @@
 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);
Index: /issm/trunk-jpl/examples/Jakobshavn/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Jakobshavn/runme.m	(revision 17068)
+++ /issm/trunk-jpl/examples/Jakobshavn/runme.m	(revision 17069)
@@ -1,3 +1,3 @@
-steps=[1];
+steps=[2:3];
 
 if any(steps==1)
@@ -52,5 +52,5 @@
 
 	%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;
Index: /issm/trunk-jpl/examples/SquareIceShelf/Square.par
===================================================================
--- /issm/trunk-jpl/examples/SquareIceShelf/Square.par	(revision 17068)
+++ /issm/trunk-jpl/examples/SquareIceShelf/Square.par	(revision 17069)
@@ -25,4 +25,5 @@
 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');