Index: /issm/trunk/src/c/objects/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Tria.cpp	(revision 76)
+++ /issm/trunk/src/c/objects/Tria.cpp	(revision 77)
@@ -336,5 +336,5 @@
 
 	double MAXSLOPE=.06; // 6 %
-	double MOUNTAINKEXPONENT=10;
+	double MOUNTAINKEXPONENT=12;
 
 
Index: /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m
===================================================================
--- /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m	(revision 76)
+++ /issm/trunk/src/m/classes/@triaelem/CreateKMatrix.m	(revision 77)
@@ -207,5 +207,5 @@
 %some parameters
 MAXSLOPE=.06;  %any element with slope>MAXSLOPE is considered a "rock" element, with infinite stiffness.
-MOUNTAINKEXPONENT=10; % "infinite" stiffness is going to be  10^MOUNTAINKEXPONENT
+MOUNTAINKEXPONENT=12; % "infinite" stiffness is going to be  10^MOUNTAINKEXPONENT
 
 %Create elementary stiffness matrix 
@@ -323,5 +323,5 @@
 	frictionparameters.k=K_list;
 
-	alpha2=Getalpha2(frictionparameters);
+	alpha2_list=Getalpha2(frictionparameters);
 	if (element_debug & triaelem.id==element_debugid),
 		disp(sprintf('   alpha2_list (%g %g %g )',alpha2_list(1),alpha2_list(2),alpha2_list(3)));
@@ -368,5 +368,5 @@
 
 		if (slope_magnitude>MAXSLOPE),
-			K=10^(-MOUNTAINKEXPONENT);
+			alpha2_list=10^(MOUNTAINKEXPONENT)*ones(3,1);
 		end
 	end
@@ -398,5 +398,4 @@
 		matice.B=B_param; clear B_param.
 	end
-
 
 	%Get viscosity at last two iterations: 
@@ -448,5 +447,5 @@
  
 		%compute alpha2 for the current gaussian point
-		alpha2_g=GetParameterValue(triaelem,alpha2,gauss_l1l2l3);
+		alpha2_g=GetParameterValue(triaelem,alpha2_list,gauss_l1l2l3);
         
 		if velocity_is_present
@@ -474,7 +473,6 @@
 		L
 	end
-
 end %for ig=1:num_gauss,
-	
+
 if (element_debug & triaelem.id==element_debugid),
 	disp(sprintf('Ke_gg->terms:\n'));