Index: ../trunk-jpl/test/NightlyRun/test275.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test275.m	(revision 18822)
+++ ../trunk-jpl/test/NightlyRun/test275.m	(revision 18823)
@@ -15,16 +15,14 @@
 md.damage.stabilization=1;
 
 md.damage.requested_outputs={'default','DamageF'};
-md.transient.isdamageevolution=1;
 
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',3);
-%md=solve(md,DamageEvolutionSolutionEnum());
+md=solve(md,DamageEvolutionSolutionEnum());
 
-%field_names={'D','F'};
-%field_tolerances={1e-13,1e-13};
-%field_values={...
-%		(md.results.DamageEvolutionSolution.DamageDbar),...
-%		(md.results.DamageEvolutionSolution.DamageF),...
-%	};
-md=solve(md,TransientSolutionEnum());
+field_names={'D','F'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+		(md.results.DamageEvolutionSolution.DamageDbar),...
+		(md.results.DamageEvolutionSolution.DamageF),...
+	};
Index: ../trunk-jpl/examples/ISMIP/CheatyRunme.m
===================================================================
--- ../trunk-jpl/examples/ISMIP/CheatyRunme.m	(revision 18822)
+++ ../trunk-jpl/examples/ISMIP/CheatyRunme.m	(revision 18823)
@@ -2,7 +2,7 @@
 %step 7 is specific to ISMIPA
 %step 8 is specific to ISMIPF
 
-steps=[1:4];
+steps=[1:7];
 
 % parameter file to be used, choose between CheatyIsmipA.par or CheatyIsmipF.par
 ParamFile='CheatyIsmipF.par'
Index: ../trunk-jpl/examples/SquareIceShelf/Square.par
===================================================================
--- ../trunk-jpl/examples/SquareIceShelf/Square.par	(revision 18822)
+++ ../trunk-jpl/examples/SquareIceShelf/Square.par	(revision 18823)
@@ -24,7 +24,7 @@
 disp('      creating flow law parameter');
 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=0.9*ones(md.mesh.numberofvertices,1);
+md.damage.D=zeros(md.mesh.numberofvertices,1);
 
 disp('      creating boundary conditions');
 md=SetIceShelfBC(md,'Front.exp');
Index: ../trunk-jpl/examples/SquareIceShelf/runme.m
===================================================================
--- ../trunk-jpl/examples/SquareIceShelf/runme.m	(revision 18822)
+++ ../trunk-jpl/examples/SquareIceShelf/runme.m	(revision 18823)
@@ -1,38 +1,7 @@
 md=model;
-
 md=triangle(md,'DomainOutline.exp',100000);
 md=setmask(md,'all','');
-
-md.materials=matdamageice();
 md=parameterize(md,'Square.par');
-md.initialization.vx = ones(md.mesh.numberofvertices,1);
-md.initialization.vy = 0*ones(md.mesh.numberofvertices,1);
-md.initialization.pressure = 0*ones(md.mesh.numberofvertices,1);
-md.damage.isdamage=1;
-md.damage.law=1;
-md.damage.max_damage=0.8;
-md.damage.D=0.5*ones(md.mesh.numberofvertices,1);
-%md.damage.spcdamage=NaN*ones(md.mesh.numberofvertices,1);
-md.damage.c1=1;
-md.damage.c2=1;
-md.damage.equiv_stress=0;
-md.damage.stabilization=2;
-
-%md=extrude(md,7,3);
-md.surfaceforcings.mass_balance = zeros(md.mesh.numberofvertices,1);
-md.basalforcings.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1);
-md.basalforcings.floatingice_melting_rate = zeros(md.mesh.numberofvertices,1);
 md=setflowequation(md,'SSA','all');
-
-md.timestepping.start_time = 0;
-md.timestepping.final_time = 1.;
-md.timestepping.time_step  = 0.1;
-
-md.transient.isdamageevolution=1;
-md.transient.isthermal=0;
-
 md.cluster=generic('name',oshostname,'np',2);
-md.stressbalance.requested_outputs={'default','DeviatoricStressxx','DeviatoricStressyy','DeviatoricStressxy'};
-%md=solve(md,StressbalanceSolutionEnum());
-md.verbose.solution = 1;
-md=solve(md,TransientSolutionEnum());
+md=solve(md,StressbalanceSolutionEnum());
