Index: /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp	(revision 19930)
+++ /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp	(revision 19931)
@@ -484,4 +484,5 @@
 		IssmDouble pressure_water = rho_water*g*(head-bed);
 		if(pressure_water>pressure_ice) pressure_water = pressure_ice;
+      
 
 		meltrate = 1/latentheat*(G+frictionheat+rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1]));
@@ -501,4 +502,9 @@
 	if(newgap<mingap) newgap=mingap;
 
+	/*Limit gap height to grow to surface*/
+	if(newgap>thickness)
+	 newgap = thickness;
+	
+	 
 	/*Add new gap as an input*/
 	element->AddInput(HydrologyGapHeightEnum,&newgap,P0Enum);