Index: /issm/trunk-jpl/test/Par/79North.par
===================================================================
--- /issm/trunk-jpl/test/Par/79North.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/79North.par	(revision 13137)
@@ -13,5 +13,4 @@
 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/trunk-jpl/test/Par/ISMIPA.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPA.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/ISMIPA.par	(revision 13137)
@@ -17,5 +17,4 @@
 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/trunk-jpl/test/Par/ISMIPB.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPB.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/ISMIPB.par	(revision 13137)
@@ -17,5 +17,4 @@
 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/trunk-jpl/test/Par/ISMIPC.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPC.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/ISMIPC.par	(revision 13137)
@@ -18,5 +18,4 @@
 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/trunk-jpl/test/Par/ISMIPD.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPD.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/ISMIPD.par	(revision 13137)
@@ -17,5 +17,4 @@
 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/trunk-jpl/test/Par/ISMIPE.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPE.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/ISMIPE.par	(revision 13137)
@@ -26,5 +26,4 @@
 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/trunk-jpl/test/Par/ISMIPF.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPF.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/ISMIPF.par	(revision 13137)
@@ -16,5 +16,4 @@
 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/trunk-jpl/test/Par/Pig.par
===================================================================
--- /issm/trunk-jpl/test/Par/Pig.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/Pig.par	(revision 13137)
@@ -17,5 +17,4 @@
 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/trunk-jpl/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetEISMINT.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/RoundSheetEISMINT.par	(revision 13137)
@@ -20,5 +20,4 @@
 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/trunk-jpl/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetShelf.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/RoundSheetShelf.par	(revision 13137)
@@ -60,5 +60,4 @@
 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/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/RoundSheetStaticEISMINT.par	(revision 13137)
@@ -25,5 +25,4 @@
 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/trunk-jpl/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareEISMINT.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/SquareEISMINT.par	(revision 13137)
@@ -27,5 +27,4 @@
 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/trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 13137)
@@ -22,5 +22,4 @@
 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/trunk-jpl/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetShelf.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/SquareSheetShelf.par	(revision 13137)
@@ -25,5 +25,4 @@
 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/trunk-jpl/test/Par/SquareShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/SquareShelf.par	(revision 13137)
@@ -22,5 +22,4 @@
 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/trunk-jpl/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 13137)
@@ -22,5 +22,4 @@
 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/trunk-jpl/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareThermal.par	(revision 13136)
+++ /issm/trunk-jpl/test/Par/SquareThermal.par	(revision 13137)
@@ -28,5 +28,4 @@
 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');