Index: ../trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 19003)
+++ ../trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 19004)
@@ -466,7 +466,7 @@
 				residual[i] = 0.;
 			}
 			//adding base in min to take into account heads under bed wich don't change N
-			pressure[i]=(rho_ice*g*thickness[i])-(rho_freshwater*g*(min(h_max,values[i],base[i])-base[i]));
+			pressure[i]=(rho_ice*g*thickness[i])-(rho_freshwater*g*(min((min(h_max,values[i])-base[i]),0.0)));
 		}
 		xDelete<IssmDouble>(thickness);
 		xDelete<IssmDouble>(base);
