Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 24278)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 24279)
@@ -1198,5 +1198,5 @@
 			dfall = min(max(1.29 - 0.17*V,0.20),1.0);
 			sfall = min(max(0.08*V + 0.38,0.5),0.9);
-			refall = (0.1 + (1.0-dfall)*0.25 + (0.5-sfall)*0.1)/2.0;
+			refall = max((0.1 + (1.0-dfall)*0.25 + (0.5-sfall)*0.1)/2.0,Gdntol);
 
 			// if snow depth is greater than specified min dz, new cell created
@@ -1228,5 +1228,5 @@
 				gdn[0] = (dfall * P + gdn[0] * mInit[0])/mass;
 				gsp[0] = (sfall * P + gsp[0] * mInit[0])/mass;
-				refall = (0.1 + (1.0-dfall)*0.25 + (0.5-sfall)*0.1)/2.0;
+				re[0] = max((0.1 + (1.0-gdn[0])*0.25 + (0.5-gsp[0])*0.1)/2.0,Gdntol);
 			}
 		}