Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 18182)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 18183)
@@ -559,5 +559,5 @@
 				if(EPL_N<0.0)EPL_N=0.0;
 				/*Get then the square of the gradient of EPL heads*/
-				EPLgrad2 = (epl_slopeX[i]*epl_slopeX[i]+epl_slopeY[i]*epl_slopeY[i]);
+				EPLgrad2 = (epl_slopeX[i]+epl_slopeY[i])*(epl_slopeX[i]+epl_slopeY[i]);
 				
 				/*And proceed to the real thing*/
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 18182)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 18183)
@@ -604,5 +604,7 @@
 		sed_trans_input->GetInputValue(&sediment_transmitivity,gauss);
 		element->FindParam(&leakage,HydrologydcLeakageFactorEnum);
+		
 		transfer=(sediment_transmitivity)/(sediment_thickness*leakage);
+		
 		/* if(epl_head>sed_head){ */
 		/* 	h_max=GetHydrologyDCInefficientHmax(element,gauss,thickness_input,base_input); */
@@ -619,4 +621,5 @@
 		/* 	//transfer=(sediment_storing*sediment_transmitivity)/(sediment_thickness*leakage); */
 		/* } */
+		
 		break;
 	default: