Index: /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticBaseVert.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticBaseVert.m	(revision 120)
+++ /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticBaseVert.m	(revision 121)
@@ -56,5 +56,4 @@
 [materials(end).constants]=deal(matpar);
 materials(end)=SetStructureField(materials(end),'constants','g',md.g);
-materials(end)=SetStructureField(materials(end),'constants','viscosity_overshoot',md.viscosity_overshoot);
 materials(end)=SetStructureField(materials(end),'constants','rho_ice',md.rho_ice);
 materials(end)=SetStructureField(materials(end),'constants','rho_water',md.rho_water);
Index: /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticHoriz.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticHoriz.m	(revision 120)
+++ /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticHoriz.m	(revision 121)
@@ -46,4 +46,5 @@
 		elements(pos)=SetStructureField(elements(pos),'element','b',md.bed(md.elements(el3pos,1:3)));
 		elements(pos)=SetStructureField(elements(pos),'element','friction_type',md.drag_type);
+		elements(pos)=SetStructureField(elements(pos),'element','viscosity_overshoot',md.viscosity_overshoot);
 		elements(pos)=SetStructureField(elements(pos),'element','k',md.drag(md.elements(el3pos,1:3)));
 		elements(pos)=SetStructureField(elements(pos),'element','p',md.p(el3pos));
@@ -78,4 +79,5 @@
 		elements(pos)=SetStructureField(elements(pos),'element','b',md.bed(md.elements(el6pos,1:6)));
 		elements(pos)=SetStructureField(elements(pos),'element','friction_type',md.drag_type);
+		elements(pos)=SetStructureField(elements(pos),'element','viscosity_overshoot',md.viscosity_overshoot);
 		elements(pos)=SetStructureField(elements(pos),'element','k',md.drag(md.elements(el6pos,1:6)));
 		elements(pos)=SetStructureField(elements(pos),'element','p',md.p(el6pos));
@@ -108,5 +110,4 @@
 [materials(end).constants]=deal(matpar);
 materials(end)=SetStructureField(materials(end),'constants','g',md.g);
-materials(end)=SetStructureField(materials(end),'constants','viscosity_overshoot',md.viscosity_overshoot);
 materials(end)=SetStructureField(materials(end),'constants','rho_ice',md.rho_ice);
 materials(end)=SetStructureField(materials(end),'constants','rho_water',md.rho_water);
Index: /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticStokes.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticStokes.m	(revision 120)
+++ /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticStokes.m	(revision 121)
@@ -41,4 +41,5 @@
 elements(stokesnewnumber)=SetStructureField(elements(stokesnewnumber),'element','b',md.bed(md.elements(pos,1:6)));
 elements(stokesnewnumber)=SetStructureField(elements(stokesnewnumber),'element','friction_type',md.drag_type);
+elements(stokesnewnumber)=SetStructureField(elements(stokesnewnumber),'element','viscosity_overshoot',md.viscosity_overshoot);
 elements(stokesnewnumber)=SetStructureField(elements(stokesnewnumber),'element','k',md.drag(md.elements(pos,1:6)));
 elements(stokesnewnumber)=SetStructureField(elements(stokesnewnumber),'element','p',md.p(pos));
@@ -63,5 +64,4 @@
 [materials(end).constants]=deal(matpar);
 materials(end)=SetStructureField(materials(end),'constants','g',md.g);
-materials(end)=SetStructureField(materials(end),'constants','viscosity_overshoot',md.viscosity_overshoot);
 materials(end)=SetStructureField(materials(end),'constants','rho_ice',md.rho_ice);
 materials(end)=SetStructureField(materials(end),'constants','rho_water',md.rho_water);
Index: /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticVert.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticVert.m	(revision 120)
+++ /issm/trunk/src/m/solutions/ice/ModelProcessorDiagnosticVert.m	(revision 121)
@@ -38,4 +38,5 @@
 elements(pos)=SetStructureField(elements(pos),'element','b',md.bed(md.elements(pos,1:6)));
 elements(pos)=SetStructureField(elements(pos),'element','friction_type',md.drag_type);
+elements(pos)=SetStructureField(elements(pos),'element','viscosity_overshoot',md.viscosity_overshoot);
 elements(pos)=SetStructureField(elements(pos),'element','k',md.drag(md.elements(pos,1:6)));
 elements(pos)=SetStructureField(elements(pos),'element','p',md.p(pos));
@@ -59,5 +60,4 @@
 [materials(end).constants]=deal(matpar);
 materials(end)=SetStructureField(materials(end),'constants','g',md.g);
-materials(end)=SetStructureField(materials(end),'constants','viscosity_overshoot',md.viscosity_overshoot);
 materials(end)=SetStructureField(materials(end),'constants','rho_ice',md.rho_ice);
 materials(end)=SetStructureField(materials(end),'constants','rho_water',md.rho_water);
Index: /issm/trunk/src/m/solutions/ice/SetupParams.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/SetupParams.m	(revision 120)
+++ /issm/trunk/src/m/solutions/ice/SetupParams.m	(revision 121)
@@ -12,8 +12,6 @@
 params.eps_rel=md.eps_rel;
 params.eps_abs=md.eps_abs;
-params.viscosity_overshoot=md.viscosity_overshoot;
 params.debug=md.debug;
 params.min_thermal_constraints=md.min_thermal_constraints;
 params.element_debug=md.element_debug;
 params.element_debugid=md.element_debugid;
-