Index: ../trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 19056)
+++ ../trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 19057)
@@ -702,7 +702,7 @@
 			vec_mask->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);
 		}
 
-		/*If mask was already one, keep one*/
+		/*If mask was already one, keep one or colapse*/
 		else if(old_active[i]>0.){
 			
 			vec_mask->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);
@@ -714,18 +714,12 @@
 					vec_mask->SetValue(basalelement->nodes[i]->Sid(),0.,INS_VAL);
 					epl_thickness[i]=init_thick;
 				}
-				else{
-					vec_mask->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);
-				}
 			}
 		}
 		/*Increase of the efficient system is needed if the epl head reach the maximum value (sediment max value for now)*/
 		GetHydrologyDCInefficientHmax(&h_max,basalelement,basalelement->nodes[i]);
 		if(eplhead[i]>=h_max && active_element){
 			for(j=0;j<numnodes;j++){
-				if(old_active[j]>0.){
-					vec_mask->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);
-				}
 				/*Increase of the domain is on the downstream node in term of sediment head*/
 				if(sedhead[j] == sedheadmin){
 					if(old_active[j]==0) eplzigzag_counter[basalelement->nodes[j]->Lid()] ++;
