Index: /issm/trunk/src/m/classes/@matpar/display.m
===================================================================
--- /issm/trunk/src/m/classes/@matpar/display.m	(revision 119)
+++ /issm/trunk/src/m/classes/@matpar/display.m	(revision 120)
@@ -14,4 +14,3 @@
 disp(sprintf('   mixed_layer_capacity: %i',matpar.mixed_layer_capacity));
 disp(sprintf('   thermal_exchange_velocity: %i',matpar.thermal_exchange_velocity));
-disp(sprintf('   viscosity_overshoot: %i',matpar.viscosity_overshoot));
 disp(sprintf('   g: %g',matpar.g));
Index: /issm/trunk/src/m/classes/@matpar/matpar.m
===================================================================
--- /issm/trunk/src/m/classes/@matpar/matpar.m	(revision 119)
+++ /issm/trunk/src/m/classes/@matpar/matpar.m	(revision 120)
@@ -18,5 +18,4 @@
 	matpar.thermal_exchange_velocity=NaN;
 	matpar.g=9.81;
-	matpar.viscosity_overshoot=0.3;
 	matpar=class(matpar,'matpar');
 				
Index: /issm/trunk/src/m/classes/@pentaelem/CreateKMatrix.m
===================================================================
--- /issm/trunk/src/m/classes/@pentaelem/CreateKMatrix.m	(revision 119)
+++ /issm/trunk/src/m/classes/@pentaelem/CreateKMatrix.m	(revision 120)
@@ -142,5 +142,4 @@
 %recover material parameters
 gravity=matpar.g;
-viscosity_overshoot=matpar.viscosity_overshoot;
 rho_ice=matpar.rho_ice;
 rho_water=matpar.rho_water;
@@ -732,5 +731,4 @@
 	%recover material parameters
 	gravity=matpar.g;
-	viscosity_overshoot=matpar.viscosity_overshoot;
 	rho_ice=matpar.rho_ice;
 	rho_water=matpar.rho_water;
@@ -841,5 +839,5 @@
 		newviscosity=GetViscosity3d(matice,epsilon);
 		oldviscosity=GetViscosity3d(matice,oldepsilon);
-		viscosity=newviscosity+viscosity_overshoot*(newviscosity-oldviscosity);
+		viscosity=newviscosity+pentaelem.viscosity_overshoot*(newviscosity-oldviscosity);
 		%disp(sprintf('Element id %i Viscosity: %g \n',pentaelem.id,viscosity));
 
Index: /issm/trunk/src/m/classes/@pentaelem/display.m
===================================================================
--- /issm/trunk/src/m/classes/@pentaelem/display.m	(revision 119)
+++ /issm/trunk/src/m/classes/@pentaelem/display.m	(revision 120)
@@ -13,4 +13,5 @@
 disp(sprintf('   bedrock: %g %g %g %g %g %g',pentaelem.b(1),pentaelem.b(2),pentaelem.b(3),pentaelem.b(4),pentaelem.b(5),pentaelem.b(6)));
 disp(sprintf('   friction type: %i',pentaelem.friction_type));
+disp(sprintf('   viscosity overshoot: %g',pentaelem.viscosity_overshoot));
 disp(sprintf('   p: %i',pentaelem.p));
 disp(sprintf('   q: %i',pentaelem.q));
Index: /issm/trunk/src/m/classes/@pentaelem/pentaelem.m
===================================================================
--- /issm/trunk/src/m/classes/@pentaelem/pentaelem.m	(revision 119)
+++ /issm/trunk/src/m/classes/@pentaelem/pentaelem.m	(revision 120)
@@ -16,4 +16,5 @@
 	pentaelem.b=zeros(6,1);
 	pentaelem.friction_type=2;
+	pentaelem.viscosity_overshoot=0;
 	pentaelem.k=zeros(6,1);
 	pentaelem.p=1;
Index: /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m
===================================================================
--- /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m	(revision 119)
+++ /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m	(revision 120)
@@ -218,5 +218,4 @@
 %recover material parameters
 gravity=matpar.g;
-viscosity_overshoot=matpar.viscosity_overshoot;
 rho_ice=matpar.rho_ice;
 rho_water=matpar.rho_water;
@@ -243,8 +242,6 @@
 temperature_average_list=zeros(numgrids,1);
 
-
 %Get all element grid data:
 xyz_list=getgriddata(triaelem,grids);
-
 
 %Build linear indices for elementary stiffness matrix.
@@ -298,6 +295,4 @@
 	if(temperature_average_is_present)disp(sprintf('   temperature_average [%g %g %g]',temperature_average_list(1),temperature_average_list(2),temperature_average_list(3)));end;
 end
-
-
 
 alpha2_list=zeros(3,1);
@@ -402,5 +397,5 @@
  	newviscosity=GetViscosity2d(matice,epsilon);
  	oldviscosity=GetViscosity2d(matice,oldepsilon);
-	viscosity=newviscosity+viscosity_overshoot*(newviscosity-oldviscosity);
+	viscosity=newviscosity+triaelem.viscosity_overshoot*(newviscosity-oldviscosity);
 
 	%Get Jacobian determinant:
Index: /issm/trunk/src/m/classes/@triaelem/display.m
===================================================================
--- /issm/trunk/src/m/classes/@triaelem/display.m	(revision 119)
+++ /issm/trunk/src/m/classes/@triaelem/display.m	(revision 120)
@@ -13,4 +13,5 @@
 disp(sprintf('   grids bedrock: %g %g %g',triaelem.b(1),triaelem.b(2),triaelem.b(3)));
 disp(sprintf('   friction type: %i',triaelem.friction_type));
+disp(sprintf('   viscosity overshoot: %g',triaelem.viscosity_overshoot));
 disp(sprintf('   p: %i',triaelem.p));
 disp(sprintf('   q: %i',triaelem.q));
Index: /issm/trunk/src/m/classes/@triaelem/triaelem.m
===================================================================
--- /issm/trunk/src/m/classes/@triaelem/triaelem.m	(revision 119)
+++ /issm/trunk/src/m/classes/@triaelem/triaelem.m	(revision 120)
@@ -16,4 +16,5 @@
 	triaelem.b=zeros(3,1);
 	triaelem.friction_type=2;
+	triaelem.viscosity_overshoot=0;
 	triaelem.k=zeros(3,1);
 	triaelem.p=1;