Index: /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp	(revision 19925)
+++ /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp	(revision 19926)
@@ -240,5 +240,5 @@
 
 	/*Build friction element, needed later: */
-	Friction* friction=new Friction(element,3);
+	Friction* friction=new Friction(element,2);
 
 	/* Start  looping on the number of gaussian points: */
@@ -337,4 +337,9 @@
 		if(values[i]>rho_ice*thickness[i]/rho_water+bed[i]){
 			values[i] = rho_ice*thickness[i]/rho_water+bed[i];
+		}
+
+		/*Make sure that negative pressure is not allowed*/
+      if(values[i]<bed[i]){
+			values[i] = bed[i];
 		}
 
@@ -434,5 +439,5 @@
 
 	/*Build friction element, needed later: */
-	Friction* friction=new Friction(element,3);
+	Friction* friction=new Friction(element,2);
 
 	/*Keep track of weights*/