Index: /issm/branches/trunk-jpl-damage/test/Par/79North.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/79North.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/79North.par	(revision 11467)
@@ -13,4 +13,5 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_B=paterson(md.initialization.temperature);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 
Index: /issm/branches/trunk-jpl-damage/test/Par/ISMIPA.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/ISMIPA.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/ISMIPA.par	(revision 11467)
@@ -17,4 +17,5 @@
 md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      boundary conditions for diagnostic model');
Index: /issm/branches/trunk-jpl-damage/test/Par/ISMIPB.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/ISMIPB.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/ISMIPB.par	(revision 11467)
@@ -17,4 +17,5 @@
 md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      boundary conditions for diagnostic model');
Index: /issm/branches/trunk-jpl-damage/test/Par/ISMIPC.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/ISMIPC.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/ISMIPC.par	(revision 11467)
@@ -18,4 +18,5 @@
 md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      boundary conditions for diagnostic model: ');
Index: /issm/branches/trunk-jpl-damage/test/Par/ISMIPD.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/ISMIPD.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/ISMIPD.par	(revision 11467)
@@ -17,4 +17,5 @@
 md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      boundary conditions for diagnostic model: ');
Index: /issm/branches/trunk-jpl-damage/test/Par/ISMIPE.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/ISMIPE.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/ISMIPE.par	(revision 11467)
@@ -26,4 +26,5 @@
 md.materials.rheology_B=6.8067*10^7*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      boundary conditions for diagnostic model: ');
Index: /issm/branches/trunk-jpl-damage/test/Par/ISMIPF.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/ISMIPF.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/ISMIPF.par	(revision 11467)
@@ -16,4 +16,5 @@
 md.materials.rheology_B=1.4734*10^14*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=1*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 md.materials.rheology_law='None';
 
Index: /issm/branches/trunk-jpl-damage/test/Par/Pig.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/Pig.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/Pig.par	(revision 11467)
@@ -17,4 +17,5 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_B=paterson(md.initialization.temperature);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 md.initialization.temperature=md.initialization.temperature;
Index: /issm/branches/trunk-jpl-damage/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/RoundSheetEISMINT.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/RoundSheetEISMINT.par	(revision 11467)
@@ -20,4 +20,5 @@
 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);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      creating surface mass balance');
Index: /issm/branches/trunk-jpl-damage/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/RoundSheetShelf.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/RoundSheetShelf.par	(revision 11467)
@@ -60,4 +60,5 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 %Surface mass balance and basal melting
Index: /issm/branches/trunk-jpl-damage/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/RoundSheetStaticEISMINT.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/RoundSheetStaticEISMINT.par	(revision 11467)
@@ -25,4 +25,5 @@
 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);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      creating surface mass balance');
Index: /issm/branches/trunk-jpl-damage/test/Par/SquareEISMINT.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/SquareEISMINT.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/SquareEISMINT.par	(revision 11467)
@@ -27,4 +27,5 @@
 md.materials.rheology_B=1.7687*10^8*ones(md.mesh.numberofvertices,1);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      creating surface mass balance');
Index: /issm/branches/trunk-jpl-damage/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/SquareSheetConstrained.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/SquareSheetConstrained.par	(revision 11467)
@@ -22,4 +22,5 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 %Friction
Index: /issm/branches/trunk-jpl-damage/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/SquareSheetShelf.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/SquareSheetShelf.par	(revision 11467)
@@ -25,4 +25,5 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 %Accumulation and melting
Index: /issm/branches/trunk-jpl-damage/test/Par/SquareShelf.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/SquareShelf.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/SquareShelf.par	(revision 11467)
@@ -22,4 +22,5 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 %Friction
Index: /issm/branches/trunk-jpl-damage/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/SquareShelfConstrained.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/SquareShelfConstrained.par	(revision 11467)
@@ -22,4 +22,5 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 %Surface mass balance and basal melting
Index: /issm/branches/trunk-jpl-damage/test/Par/SquareThermal.par
===================================================================
--- /issm/branches/trunk-jpl-damage/test/Par/SquareThermal.par	(revision 11466)
+++ /issm/branches/trunk-jpl-damage/test/Par/SquareThermal.par	(revision 11467)
@@ -28,4 +28,5 @@
 md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_Z=ones(md.mesh.numberofvertices,1);
 
 disp('      creating surface mass balance');