Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17361)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17362)
@@ -109,4 +109,5 @@
 
 	/*Intermediaries*/
+	bool     active_element;
 	int      meshtype;
 	Element* basalelement;
@@ -125,6 +126,9 @@
 	}
 
+	Input* active_element_input = basalelement->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input);
+	active_element_input->GetInputValue(&active_element);
+
 	/*Check that all nodes are active, else return empty matrix*/
-	if(!basalelement->AllActive()) { /*This could probably replaced by active elt mask*/
+	if(!active_element) {
 	if(meshtype!=Mesh2DhorizontalEnum){
 			basalelement->DeleteMaterials(); 
@@ -204,4 +208,5 @@
 
 	/*Intermediaries*/
+	bool     active_element;
 	int      meshtype;
 	Element* basalelement;
@@ -219,6 +224,10 @@
 		default: _error_("mesh "<<EnumToStringx(meshtype)<<" not supported yet");
 	}
+
+	Input* active_element_input = basalelement->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input);
+	active_element_input->GetInputValue(&active_element);
+
 	/*Check that all nodes are active, else return empty matrix*/
-	if(!basalelement->AllActive()) {
+	if(!active_element) {
 	if(meshtype!=Mesh2DhorizontalEnum){
 			basalelement->DeleteMaterials(); 
Index: /issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp	(revision 17361)
+++ /issm/trunk-jpl/src/c/analyses/L2ProjectionEPLAnalysis.cpp	(revision 17362)
@@ -86,5 +86,4 @@
 	bool     active_element;
 	Element* basalelement;
-	Input*   active_element_input=NULL;
 
 	/*Get basal element*/
@@ -105,11 +104,10 @@
 	}
 
-	/* active_element_input=basalelement->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input); */
-	/* active_element_input->GetInputValue(&active_element); */
+	Input* active_element_input=basalelement->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input);
+	active_element_input->GetInputValue(&active_element);
 
 	/* Check that all nodes are active, else return empty matrix */
-	if(!basalelement->AllActive()) {
-		//	if(!active_element){
-	if(meshtype!=Mesh2DhorizontalEnum){
+	if(!active_element){
+		if(meshtype!=Mesh2DhorizontalEnum){
 			basalelement->DeleteMaterials();
 			delete basalelement;
@@ -160,5 +158,4 @@
 	bool     active_element;
 	Element* basalelement;
-	Input*   active_element_input=NULL;
 
 	/*Get basal element*/
@@ -175,10 +172,9 @@
 	}
 
-	/* active_element_input=basalelement->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input); */
-	/* active_element_input->GetInputValue(&active_element); */
+	Input* active_element_input = basalelement->GetInput(HydrologydcMaskEplactiveEltEnum); _assert_(active_element_input);
+	active_element_input->GetInputValue(&active_element);
 
 	/*Check that all nodes are active, else return empty matrix*/
-	if(!basalelement->AllActive()) {
-		/* if(!active_element){ */
+	if(!active_element) {
 		if(meshtype!=Mesh2DhorizontalEnum){
 			basalelement->DeleteMaterials();
Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 17361)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 17362)
@@ -37,20 +37,20 @@
 	this->inputs->AddInput(input_in);
 }/*}}}*/
-bool Element::AllActive(void){/*{{{*/
-
-	int numnodes = this->GetNumberOfNodes();
-	for(int i=0;i<numnodes;i++){
-		if(!this->nodes[i]->IsActive()) return false;
-	}
-	return true;
-}/*}}}*/
-bool Element::AnyActive(void){/*{{{*/
-
-	int numnodes = this->GetNumberOfNodes();
-	for(int i=0;i<numnodes;i++){
-		if(this->nodes[i]->IsActive()) return true;
-	}
-	return false;
-}/*}}}*/
+/* bool Element::AllActive(void){/\*{{{*\/ */
+
+/* 	int numnodes = this->GetNumberOfNodes(); */
+/* 	for(int i=0;i<numnodes;i++){ */
+/* 		if(!this->nodes[i]->IsActive()) return false; */
+/* 	} */
+/* 	return true; */
+/* }/\*}}}*\/ */
+/* bool Element::AnyActive(void){/\*{{{*\/ */
+
+/* 	int numnodes = this->GetNumberOfNodes(); */
+/* 	for(int i=0;i<numnodes;i++){ */
+/* 		if(this->nodes[i]->IsActive()) return true; */
+/* 	} */
+/* 	return false; */
+/* }/\*}}}*\/ */
 void Element::CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes_list,int numnodes,int* cs_array){/*{{{*/
 
Index: /issm/trunk-jpl/src/c/classes/Elements/Element.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 17361)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 17362)
@@ -52,6 +52,6 @@
 		/*Functions*/
 		void       AddInput(Input* input_in);
-		bool       AllActive(void);
-		bool       AnyActive(void);
+		/* bool       AllActive(void); */
+		/* bool       AnyActive(void); */
 		void       CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes,int numnodes,int* cs_array);
 		void       DeleteMaterials(void);
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 17361)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp	(revision 17362)
@@ -117,5 +117,4 @@
 					if(num_unstable_constraints==0) sedconverged = true;
 					if (sedcount>=hydro_maxiter){
-						//sedconverged = true;
 						_error_("   maximum number of Sediment iterations (" << hydro_maxiter << ") exceeded");
 					}
@@ -174,4 +173,5 @@
 				/* {{{ *//*Retrieve the EPL head slopes and compute EPL Thickness*/
 				if(VerboseSolution()) _printf0_("computing EPL Head slope...\n");
+				inefanalysis->ElementizeEplMask(femmodel);
 				femmodel->SetCurrentConfiguration(L2ProjectionEPLAnalysisEnum);
 				femmodel->UpdateConstraintsL2ProjectionEPLx();