Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 15249)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 15250)
@@ -1589,5 +1589,5 @@
 
 	/*Step 1: update maks, the mask might be extended by residual and/or using downstream sediment head*/
-	mask=new Vector<IssmDouble>(nodes->NumberOfNodes());
+	mask=new Vector<IssmDouble>(nodes->NumberOfNodes(HydrologyDCEfficientAnalysisEnum));
 	for (int i=0;i<elements->Size();i++){
 		Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
@@ -1605,5 +1605,5 @@
 
 	/*Step 2: update node activity. If one element is connected to mask=1, all nodes are active*/
-	active=new Vector<IssmDouble>(nodes->NumberOfNodes());
+	active=new Vector<IssmDouble>(nodes->NumberOfNodes(HydrologyDCEfficientAnalysisEnum));
 	for (int i=0;i<elements->Size();i++){
 		Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
@@ -1617,14 +1617,19 @@
 
 	/*Update node activation accordingly*/
+	int counter =0;
 	for (int i=0;i<nodes->Size();i++){
 		Node* node=dynamic_cast<Node*>(nodes->GetObjectByOffset(i));
-		if(serial_active[node->Sid()]==1.){
-			node->Activate();
-		}
-		else{
-			node->Deactivate();
+		if(node->InAnalysis(HydrologyDCEfficientAnalysisEnum)){
+			if(serial_active[node->Sid()]==1.){
+				node->Activate();
+				counter++;
+			}
+			else{
+				node->Deactivate();
+			}
 		}
 	}
 	xDelete<IssmDouble>(serial_active);
+	if(VerboseSolution()) _printf0_("   Number of active nodes in EPL layer: "<< counter <<"\n");
 
 }