Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 26205)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 26206)
@@ -2732,5 +2732,5 @@
 ElementMatrix* StressbalanceAnalysis::CreateKMatrixMLHO(Element* element){/*{{{*/
 	
-	_error_("Mono Layer Higher-Order called, not fully tested. If you are sure of using it, comment this line.");
+	//_error_("Mono Layer Higher-Order called, not fully tested. If you are sure of using it, comment this line.");
 	
 	/* Check if ice in element */
@@ -2783,5 +2783,5 @@
 
 	/*build friction object, used later on: */
-	Friction* friction=new Friction(element,dim);
+	Friction* friction=new Friction(element,dim,true);
 
 	/*Recover portion of element that is grounded*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 26205)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 26206)
@@ -20,11 +20,13 @@
 	this->dim=0;
 	this->law=0;
+	this->isMLHO=false;
 
 }
 /*}}}*/
-Friction::Friction(Element* element_in,int dim_in){/*{{{*/
+Friction::Friction(Element* element_in,int dim_in,bool isMLHO_in){/*{{{*/
 
 	this->element=element_in;
 	this->dim=dim_in;
+	this->isMLHO=isMLHO_in;
 	element_in->FindParam(&this->law,FrictionLawEnum);
 }
@@ -855,9 +857,7 @@
 
 	/*diverse*/
-	bool isMLHO=false; //employing mono layer higher order model? Default is no 
 	IssmDouble vx,vy,vz,vmag;
 	
-	element->parameters->FindParam(&isMLHO,FlowequationIsMLHOEnum);
-	if(isMLHO) _assert_(dim==2);
+	if(this->isMLHO) _assert_(dim==2);
 
 	switch(dim){
@@ -867,5 +867,5 @@
 			break;
 		case 2:
-			if(isMLHO){
+			if(this->isMLHO){
 				element->GetInputValue(&vx,gauss,VxBaseEnum);
 				element->GetInputValue(&vy,gauss,VyBaseEnum);
Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.h	(revision 26205)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.h	(revision 26206)
@@ -17,8 +17,9 @@
 		int      dim;
 		int      law;
+		bool		isMLHO=false;
 
 		/*methods: */
 		Friction();
-		Friction(Element* element_in,int dim_in);
+		Friction(Element* element_in,int dim_in,bool isMLHO=0);
 		~Friction();