Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 9594)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 9595)
@@ -1842,5 +1842,5 @@
 	double     Jdet,u,v,w,um,vm,wm;
 	double     gravity,rho_ice,rho_water;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     heatcapacity,thermalconductivity,dt;
 	double     pressure,enthalpy;
@@ -2118,5 +2118,5 @@
 	int        i,j,ig,found=0;
 	double     Jdet,u,v,w,um,vm,wm;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     gravity,rho_ice,rho_water;
 	double     heatcapacity,thermalconductivity,dt;
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 9594)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 9595)
@@ -1511,5 +1511,5 @@
 	double     obs_velocity_mag,velocity_mag;
 	double     dux,duy;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
 	double     scalex=0,scaley=0,scale=0,S=0;
@@ -1687,5 +1687,5 @@
 	double     obs_velocity_mag,velocity_mag;
 	double     dux,duy;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
 	double     scalex=0,scaley=0,scale=0,S=0;
@@ -4857,5 +4857,5 @@
 	double     Jelem=0;
 	double     misfit,Jdet;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
 	double     velocity_mag,obs_velocity_mag;
@@ -4922,5 +4922,5 @@
 	int        fit=-1;
 	double     Jelem=0, S=0;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
 	double     misfit, Jdet;
@@ -5012,5 +5012,5 @@
 	double     scalex=1,scaley=1;
 	double     misfit,Jdet;
-	double     epsvel=1.e-13;
+	double     epsvel=2.220446049250313e-16;
 	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
 	double     vx,vy,vxobs,vyobs,weight;
