Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17963)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17964)
@@ -288,5 +288,5 @@
 			thickness_input ->GetInputValue(&epl_thickness,gauss);
 			
-			/*Dealing with the sediment part of the transfer term*/
+			/*Dealing with the epl part of the transfer term*/
 			transfer=GetHydrologyPVectorTransfer(basalelement,gauss,thickness_input,sed_head_input,epl_head_input,sed_trans_input,residual_input);
 			scalar = Jdet*gauss->weight*((water_head*epl_specificstoring*epl_thickness)+(transfer*dt));
@@ -410,18 +410,18 @@
 		residual_input->GetInputValue(&residual,gauss);
 		element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
-
-		if(epl_head>sed_head){
-			if(residual>0.0){	
-				transfer=0.0;
-			}
-			else{
-				transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
-				//				transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage);
-			}
-		}
-		else{
-			transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
-			//			transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage);
-		}
+		transfer=(sediment_transmitivity)/(sediment_thickness*leakage); 
+		/* if(epl_head>sed_head){ */
+		/* 	if(residual>0.0){	 */
+		/* 		transfer=0.0; */
+		/* 	} */
+		/* 	else{ */
+		/* 		transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* 		//transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* 	} */
+		/* } */
+		/* else{ */
+		/* 	transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* 	//transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* } */
 		break;
 	default:
@@ -463,18 +463,18 @@
 		residual_input->GetInputValue(&residual,gauss);
 		element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
-		
-		if(epl_head>sediment_head){
-			if(residual>0.0){	
-				transfer=0.0;
-			}
-			else{
-				transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
-				//				transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
-			}
-		}
-		else{
-			transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
-			//			transfer=(sediment_storing*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
-		}
+		transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage);
+		/* if(epl_head>sediment_head){ */
+		/* 	if(residual>0.0){	 */
+		/* 		transfer=0.0; */
+		/* 	} */
+		/* 	else{ */
+		/* 		transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
+		/* 		//transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
+		/* 	} */
+		/* } */
+		/* else{ */
+		/* 	transfer=(sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
+		/* 	//transfer=(sediment_storing*sediment_transmitivity*sediment_head)/(sediment_thickness*leakage); */
+		/* } */
 		break;
 	default:
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 17963)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 17964)
@@ -329,4 +329,5 @@
 
 		scalar = Jdet*gauss->weight*(water_load);
+
 		if(dt!=0.) scalar = scalar*dt;
 		for(int i=0;i<numnodes;i++){
@@ -601,19 +602,19 @@
 		sed_trans_input->GetInputValue(&sediment_transmitivity,gauss);
 		element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
-	
-		if(epl_head>sed_head){
-			h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input);
-			if(sed_head>=h_max){
-				transfer=0.0;
-			}
-			else{
-				transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
-				//				transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage);
-			}
-		}
-		else{
-			transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
-			//			transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage);
-		}
+		transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
+		/* if(epl_head>sed_head){ */
+		/* 	h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input); */
+		/* 	if(sed_head>=h_max){ */
+		/* 		transfer=0.0; */
+		/* 	} */
+		/* 	else{ */
+		/* 		transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* 		//transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* 	} */
+		/* } */
+		/* else{ */
+		/* 	transfer=(sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* 	//transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage); */
+		/* } */
 		break;
 	default:
@@ -658,19 +659,19 @@
 
 		element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
-
-		if(epl_head>sediment_head){
-			h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input);
-			if(sediment_head>=h_max){
-				transfer=0.0;
-			}
-			else{
-				transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
-				//				transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
-			}
-		}
-		else{
-			transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
-			//			transfer=(sediment_storing*sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
-		}
+ 		transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage);
+		/* if(epl_head>sediment_head){ */
+		/* 	h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input); */
+ 		/* 	if(sediment_head>=h_max){ */
+		/* 		transfer=0.0; */
+		/* 	} */
+		/* 	else{ */
+		/* 		transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
+		/* 		//transfer=(epl_specificstoring*epl_thickness*sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
+		/* 	} */
+		/* } */
+		/* else{ */
+		/* 	transfer=(sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
+		/* 	//transfer=(sediment_storing*sediment_transmitivity*epl_head)/(sediment_thickness*leakage); */
+		/* } */
 		break;
 	default:
Index: /issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp	(revision 17963)
+++ /issm/trunk-jpl/src/c/classes/Loads/Pengrid.cpp	(revision 17964)
@@ -646,5 +646,5 @@
 		if(zigzag_counter>penalty_lock){
 			unstable=0;
-			active=0;
+			active=1;
 		}
 	}