Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 23242)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 23243)
@@ -427,5 +427,5 @@
 			IssmDouble a0 = min(0.98, 1 - 1.58 *pow(gsz,0.5));
 			// 0.8 - 1.5um
-			IssmDouble a1 = max(0, 0.95 - 15.4 *pow(gsz,0.5));
+			IssmDouble a1 = max(0., 0.95 - 15.4 *pow(gsz,0.5));
 			// 1.5 - 2.8um
 			IssmDouble a2 = max(0.127, 0.88 + 346.3*gsz - 32.31*pow(gsz,0.5));