Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 17978)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 17979)
@@ -449,10 +449,9 @@
 			if(values[i]>h_max) {
 				residual[i] = kappa*(values[i]-h_max);
-				pressure[i]=(rho_ice*g*thickness[i])-(rho_freshwater*g*(values[i]-base[i]));
 			}
 			else{
 				residual[i] = 0.;
-				pressure[i]=(rho_ice*g*thickness[i])-(rho_freshwater*g*(h_max-base[i]));
 			}
+			pressure[i]=(rho_ice*g*thickness[i])-(rho_freshwater*g*(min(hmax,values[i])-base[i]));
 		}
 		xDelete<IssmDouble>(thickness);