Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 18726)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 18727)
@@ -1816,4 +1816,5 @@
 		
 	HydrologyDCEfficientAnalysis* effanalysis =  new HydrologyDCEfficientAnalysis();
+	HydrologyDCInefficientAnalysis* inefanalysis =  new HydrologyDCInefficientAnalysis();
 	/*Step 1: update mask, the mask might be extended by residual and/or using downstream sediment head*/
 	mask=new Vector<IssmDouble>(this->nodes->NumberOfNodes(HydrologyDCEfficientAnalysisEnum));
@@ -1834,5 +1835,5 @@
 	InputUpdateFromVectorx(this,serial_mask,HydrologydcMaskEplactiveNodeEnum,NodeSIdEnum);
 	xDelete<IssmDouble>(serial_mask);
-
+	inefanalysis->ElementizeEplMask(this);
 	/*Step 2: update node activity. If one element is connected to mask=1, all nodes are active*/
 	active=new Vector<IssmDouble>(nodes->NumberOfNodes(HydrologyDCEfficientAnalysisEnum));
@@ -1841,4 +1842,5 @@
 		effanalysis->HydrologyEPLGetActive(active,element);
 	}
+
 
 	/*Assemble and serialize*/
@@ -1854,5 +1856,5 @@
 			if(serial_active[node->Sid()]==1.){
 				node->Activate();
-				counter++;
+				if(!node->IsClone());counter++;
 			}
 			else{
@@ -1900,5 +1902,5 @@
 			if(serial_active[node->Sid()]==1.){
 				node->Activate();
-				counter++;
+				if(!node->IsClone());counter++;
 			}
 			else{