Index: /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 24068)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 24069)
@@ -4,4 +4,6 @@
 #include "../shared/shared.h"
 #include "../modules/modules.h"
+
+#define AEPS 2.22e-15
 
 /*Model processing*/
@@ -256,5 +258,5 @@
 		/*Get norm of gradient of hydraulic potential and make sure it is >0*/
 		IssmDouble normgradphi = sqrt(dphi[0]*dphi[0] + dphi[1]*dphi[1]);
-		if(normgradphi < 1.e-12) normgradphi = 1.e-12;
+		if(normgradphi < AEPS) normgradphi = AEPS;
 
 		IssmDouble coeff = k*pow(h,alpha)*pow(normgradphi,beta-2.);
@@ -501,5 +503,5 @@
 
 		/*Make sure it is positive*/
-		if(h_new[iv]<1.e-12) h_new[iv] = 1.e-12;
+		if(h_new[iv]<AEPS) h_new[iv] = AEPS;
 	}
 
Index: /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp	(revision 24068)
+++ /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp	(revision 24069)
@@ -25,4 +25,5 @@
 #define ALPHA_S     5./4.
 #define BETA_S      3./2.
+#define AEPS 2.22e-15
 
 /*Channel constructors and destructor*/
@@ -423,5 +424,5 @@
 		dphimds = rho_water*g*(db[0]*tx + db[1]*ty);
 		Ngrad   = fabs(dphids);
-		if(Ngrad<1.e-12) Ngrad = 1.e-12;
+		if(Ngrad<AEPS) Ngrad = AEPS;
 
 		/*Compute the effective conductivity Kc = k h^alpha |grad Phi|^{beta-2} (same for sheet)*/
@@ -550,5 +551,5 @@
 		dphimds = rho_water*g*(db[0]*tx + db[1]*ty);
 		Ngrad   = fabs(dphids);
-		if(Ngrad<1.e-12) Ngrad = 1.e-12;
+		if(Ngrad<AEPS) Ngrad = AEPS;
 
 		/*Compute the effective conductivity Ks = k h^alpha |grad Phi|^{beta-2} (same for sheet)*/
@@ -652,5 +653,5 @@
 	dphimds = rho_water*g*(db[0]*tx + db[1]*ty);
 	Ngrad   = fabs(dphids);
-	if(Ngrad<1.e-12) Ngrad = 1.e-12;
+	if(Ngrad<AEPS) Ngrad = AEPS;
 
 	/*d(phi - phi_m)/ds*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp	(revision 24068)
+++ /issm/trunk-jpl/src/c/classes/Loads/Moulin.cpp	(revision 24069)
@@ -230,5 +230,5 @@
 /*}}}*/
 bool  Moulin::IsPenalty(void){/*{{{*/
-	return true;
+	return false;
 }
 /*}}}*/