Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetShelf.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetShelf.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetShelf.par	(revision 13137)
@@ -59,7 +59,6 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 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
 md.surfaceforcings.mass_balance=-10*ones(md.mesh.numberofvertices,1);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetEISMINT.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetEISMINT.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetEISMINT.par	(revision 13137)
@@ -19,7 +19,6 @@
 disp('      creating flow law paramter');
 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');
 smb_max=0.5; %m/yr
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/79North.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/79North.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/79North.par	(revision 13137)
@@ -12,7 +12,6 @@
 %Materials
 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);
 
 %Friction
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetShelf.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetShelf.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetShelf.par	(revision 13137)
@@ -24,7 +24,6 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 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
 md.surfaceforcings.mass_balance=10*ones(md.mesh.numberofvertices,1);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelfConstrained.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 13137)
@@ -21,7 +21,6 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 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
 md.surfaceforcings.mass_balance=10*ones(md.mesh.numberofvertices,1);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareThermal.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareThermal.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareThermal.par	(revision 13137)
@@ -27,7 +27,6 @@
 disp('      creating flow law paramter');
 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');
 md.surfaceforcings.mass_balance=ones(md.mesh.numberofvertices,1)/md.constants.yts; %1m/a
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/Pig.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/Pig.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/Pig.par	(revision 13137)
@@ -16,7 +16,6 @@
 %Materials
 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: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetStaticEISMINT.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/RoundSheetStaticEISMINT.par	(revision 13137)
@@ -24,7 +24,6 @@
 disp('      creating flow law paramter');
 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');
 smb_max=0.5; %m/yr
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPA.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPA.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPA.par	(revision 13137)
@@ -16,7 +16,6 @@
 disp('      creating flow law paramter');
 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');
 %Create node on boundary fist (because we cannot use mesh)
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelf.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelf.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelf.par	(revision 13137)
@@ -21,7 +21,6 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 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
 pos=find(md.mask.elementonfloatingice);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPB.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPB.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPB.par	(revision 13137)
@@ -16,7 +16,6 @@
 disp('      creating flow law paramter');
 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');
 %Create node on boundary fist (because we cannot use mesh)
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 13137)
@@ -21,7 +21,6 @@
 md.initialization.temperature=(273-20)*ones(md.mesh.numberofvertices,1);
 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
 pos=find(md.mask.elementonfloatingice);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPC.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPC.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPC.par	(revision 13137)
@@ -17,7 +17,6 @@
 disp('      creating flow law paramter');
 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: ');
 %Create node on boundary fist (because wi can not use mesh)
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPD.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPD.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPD.par	(revision 13137)
@@ -16,7 +16,6 @@
 disp('      creating flow law paramter');
 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: ');
 %Create node on boundary fist (because wi can not use mesh)
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPE.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPE.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPE.par	(revision 13137)
@@ -25,7 +25,6 @@
 disp('      creating flow law paramter');
 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: ');
 %Create node on boundary fist (because wi can not use mesh)
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPF.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPF.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPF.par	(revision 13137)
@@ -15,7 +15,6 @@
 disp('      creating flow law paramter');
 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';
 
 disp('      boundary conditions for diagnostic model');
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareEISMINT.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareEISMINT.par	(revision 13136)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareEISMINT.par	(revision 13137)
@@ -26,7 +26,6 @@
 disp('      creating flow law paramter');
 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');
 md.surfaceforcings.mass_balance=0.2*ones(md.mesh.numberofvertices,1); %0m/a