Index: /issm/trunk-jpl/test/Par/79North.par
===================================================================
--- /issm/trunk-jpl/test/Par/79North.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/79North.par	(revision 17767)
@@ -20,4 +20,7 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+
+%Damage
+md.damage.isdamage=0;
 
 %Friction
Index: /issm/trunk-jpl/test/Par/79North.py
===================================================================
--- /issm/trunk-jpl/test/Par/79North.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/79North.py	(revision 17767)
@@ -33,4 +33,7 @@
 md.initialization.temperature=md.initialization.temperature
 
+#Damage
+md.damage.isdamage=0
+
 #Friction
 md.friction.coefficient=50.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par
===================================================================
--- /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par	(revision 17767)
@@ -44,4 +44,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Friction
 md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par
===================================================================
--- /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par	(revision 17767)
@@ -43,4 +43,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Friction
 md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/ISMIPA.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPA.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPA.par	(revision 17767)
@@ -16,4 +16,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      boundary conditions for stressbalance model');
 %Create node on boundary first (because we cannot use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPA.py
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPA.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPA.py	(revision 17767)
@@ -19,4 +19,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      boundary conditions for stressbalance model"
 #Create node on boundary first (because we cannot use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPB.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPB.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPB.par	(revision 17767)
@@ -16,4 +16,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      boundary conditions for stressbalance model');
 %Create node on boundary first (because we cannot use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPB.py
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPB.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPB.py	(revision 17767)
@@ -19,4 +19,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      boundary conditions for stressbalance model"
 #Create node on boundary first (because we cannot use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPC.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPC.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPC.par	(revision 17767)
@@ -17,4 +17,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      boundary conditions for stressbalance model:');
 %Create node on boundary first (because we can not use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPC.py
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPC.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPC.py	(revision 17767)
@@ -20,4 +20,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      boundary conditions for stressbalance model:"
 #Create node on boundary first (because we can not use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPD.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPD.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPD.par	(revision 17767)
@@ -16,4 +16,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      boundary conditions for stressbalance model:');
 %Create node on boundary first (because we can not use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPD.py
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPD.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPD.py	(revision 17767)
@@ -19,4 +19,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      boundary conditions for stressbalance model:"
 #Create node on boundary first (because we can not use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPE.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPE.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPE.par	(revision 17767)
@@ -26,4 +26,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      boundary conditions for stressbalance model:');
 %Create node on boundary first (because we can not use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPE.py
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPE.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPE.py	(revision 17767)
@@ -31,4 +31,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      boundary conditions for stressbalance model:"
 #Create node on boundary first (because we can not use mesh)
Index: /issm/trunk-jpl/test/Par/ISMIPF.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPF.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPF.par	(revision 17767)
@@ -17,4 +17,7 @@
 md.materials.rheology_n=1.*ones(md.mesh.numberofelements,1);
 md.materials.rheology_law='None';
+
+%Damage
+md.damage.isdamage=0;
 
 disp('      boundary conditions for stressbalance model');
Index: /issm/trunk-jpl/test/Par/ISMIPF.py
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPF.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ISMIPF.py	(revision 17767)
@@ -20,4 +20,7 @@
 md.materials.rheology_n=1.*numpy.ones((md.mesh.numberofelements,1))
 md.materials.rheology_law='None'
+
+#Damage
+md.damage.isdamage=0
 
 print "      boundary conditions for stressbalance model"
Index: /issm/trunk-jpl/test/Par/Pig.par
===================================================================
--- /issm/trunk-jpl/test/Par/Pig.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/Pig.par	(revision 17767)
@@ -26,4 +26,7 @@
 md.initialization.temperature=md.initialization.temperature;
 
+%Damage
+md.damage.isdamage=0;
+
 %Friction
 md.friction.coefficient=50*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/Pig.py
===================================================================
--- /issm/trunk-jpl/test/Par/Pig.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/Pig.py	(revision 17767)
@@ -37,4 +37,7 @@
 md.initialization.temperature=md.initialization.temperature
 
+#Damage
+md.damage.isdamage=0
+
 #Friction
 md.friction.coefficient=50.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetEISMINT.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/RoundSheetEISMINT.par	(revision 17767)
@@ -20,4 +20,7 @@
 md.materials.rheology_B=6.81*10^7*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution 
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
+
+%Damage
+md.damage.isdamage=0;
 
 disp('      creating surface mass balance');
Index: /issm/trunk-jpl/test/Par/RoundSheetEISMINT.py
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetEISMINT.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/RoundSheetEISMINT.py	(revision 17767)
@@ -23,4 +23,7 @@
 md.materials.rheology_B=6.81*10**7*numpy.ones((md.mesh.numberofvertices,1))    #to have the same B as the analytical solution 
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
+
+#Damage
+md.damage.isdamage=0
 
 print "      creating surface mass balance"
Index: /issm/trunk-jpl/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetShelf.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/RoundSheetShelf.par	(revision 17767)
@@ -61,4 +61,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Surface mass balance and basal melting
 md.surfaceforcings.mass_balance=-10.*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/RoundSheetShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetShelf.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/RoundSheetShelf.py	(revision 17767)
@@ -68,4 +68,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 #Surface mass balance and basal melting
 md.surfaceforcings.mass_balance=-10.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par	(revision 17767)
@@ -24,4 +24,7 @@
 md.materials.rheology_B=6.81*10^7*ones(md.mesh.numberofvertices,1); %to have the same B as the analytical solution 
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
+
+%Damage
+md.damage.isdamage=0;
 
 disp('      creating surface mass balance');
Index: /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.py
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.py	(revision 17767)
@@ -28,4 +28,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      creating surface mass balance"
 smb_max=0.5    #m/yr
Index: /issm/trunk-jpl/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareEISMINT.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareEISMINT.par	(revision 17767)
@@ -26,4 +26,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      creating surface mass balance');
 md.surfaceforcings.mass_balance=0.2*ones(md.mesh.numberofvertices,1); %0m/a
Index: /issm/trunk-jpl/test/Par/SquareEISMINT.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareEISMINT.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareEISMINT.py	(revision 17767)
@@ -29,4 +29,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      creating surface mass balance"
 md.surfaceforcings.mass_balance=0.2*numpy.ones((md.mesh.numberofvertices,1))    #0m/a
Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 17767)
@@ -29,4 +29,7 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
+
+%Damage
+md.damage.isdamage=0;
 
 md.masstransport.calvingrate=zeros(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.py	(revision 17767)
@@ -40,4 +40,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 #Masstransport
 md.masstransport.calvingrate=0.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetShelf.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareSheetShelf.par	(revision 17767)
@@ -32,4 +32,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Accumulation and melting
 md.surfaceforcings.mass_balance=10.*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/SquareSheetShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetShelf.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareSheetShelf.py	(revision 17767)
@@ -43,4 +43,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 #Accumulation and melting
 md.surfaceforcings.mass_balance=10.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/SquareShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareShelf.par	(revision 17767)
@@ -29,4 +29,7 @@
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Friction
 md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/SquareShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareShelf.py	(revision 17767)
@@ -59,4 +59,7 @@
 md.materials.rheology_n = 3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 #Friction
 md.friction.coefficient = 20.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/SquareShelf2.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf2.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareShelf2.par	(revision 17767)
@@ -29,4 +29,7 @@
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Friction
 md.friction.coefficient=20*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/SquareShelf2.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf2.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareShelf2.py	(revision 17767)
@@ -59,4 +59,7 @@
 md.materials.rheology_n = 3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 #Friction
 md.friction.coefficient = 20.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 17767)
@@ -29,4 +29,7 @@
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 %Surface mass balance and basal melting
 md.surfaceforcings.mass_balance=10*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/SquareShelfConstrained.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelfConstrained.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareShelfConstrained.py	(revision 17767)
@@ -40,4 +40,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 #Surface mass balance and basal melting
 md.surfaceforcings.mass_balance=10.*numpy.ones((md.mesh.numberofvertices,1))
Index: /issm/trunk-jpl/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareThermal.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareThermal.par	(revision 17767)
@@ -27,4 +27,7 @@
 md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 
+%Damage
+md.damage.isdamage=0;
+
 disp('      creating surface mass balance');
 md.surfaceforcings.mass_balance=ones(md.mesh.numberofvertices,1)/md.constants.yts; %1m/a
Index: /issm/trunk-jpl/test/Par/SquareThermal.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareThermal.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/SquareThermal.py	(revision 17767)
@@ -31,4 +31,7 @@
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
 
+#Damage
+md.damage.isdamage=0
+
 print "      creating surface mass balance"
 md.surfaceforcings.mass_balance=numpy.ones((md.mesh.numberofvertices,1))/md.constants.yts    #1m/a
Index: /issm/trunk-jpl/test/Par/ValleyGlacierShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/ValleyGlacierShelf.par	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ValleyGlacierShelf.par	(revision 17767)
@@ -40,4 +40,7 @@
 md.materials.rheology_B = paterson(md.initialization.temperature);
 md.materials.rheology_n = 3.*ones(md.mesh.numberofelements,1);
+
+%Damage
+md.damage.isdamage=0;
 
 %Thermal
Index: /issm/trunk-jpl/test/Par/ValleyGlacierShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/ValleyGlacierShelf.py	(revision 17766)
+++ /issm/trunk-jpl/test/Par/ValleyGlacierShelf.py	(revision 17767)
@@ -47,4 +47,7 @@
 md.materials.rheology_B=paterson(md.initialization.temperature)
 md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
+
+#Damage
+md.damage.isdamage=0
 
 #Thermal