Index: /issm/trunk/src/c/ModelProcessorx/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp	(revision 109)
+++ /issm/trunk/src/c/ModelProcessorx/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp	(revision 110)
@@ -66,4 +66,5 @@
 	double tria_meanvel;/*!scaling ratio for velocities*/
 	double tria_epsvel; /*!minimum velocity to avoid infinite velocity ratios*/
+	double tria_viscosity_overshoot; 
 
 	/*matice constructor input: */
@@ -274,6 +275,9 @@
 			tria_epsvel=model->epsvel;
 
+			/*viscosity_overshoot*/
+			tria_viscosity_overshoot=model->viscosity_overshoot;
+
 			/*Create tria element using its constructor:*/
-			tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel);
+			tria=new Tria(tria_id, tria_mid, tria_mparid, tria_g, tria_h, tria_s, tria_b, tria_k, tria_friction_type, tria_p, tria_q, tria_shelf, tria_meanvel, tria_epsvel, tria_viscosity_overshoot);
 
 			/*Add tria element to elements dataset: */
@@ -571,8 +575,8 @@
 		node_x[1]=model->y[i];
 		node_x[2]=model->z[i];
+
 		
 		node_onbed=(int)model->gridonbed[i];
-		node_onsurface=(int)model->gridonsurface[i];
-		
+		node_onsurface=(int)model->gridonsurface[i];	
 		if (strcmp(model->meshtype,"3d")==0){
 			if (isnan(model->uppernodes[i])){