Index: /issm/trunk-jpl/src/c/objects/Kriging/ExponentialVariogram.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Kriging/ExponentialVariogram.cpp	(revision 12366)
+++ /issm/trunk-jpl/src/c/objects/Kriging/ExponentialVariogram.cpp	(revision 12367)
@@ -71,5 +71,8 @@
 	h=sqrt(pow(deltax,2.)+pow(deltay,2.));
 
-	/*return covariance*/
+	/*If h is too small, return sill*/
+	if(h<0.0000001) return sill;
+
+	/*compute covariance*/
 	a     = 1./3.;
 	cova = (sill-nugget)*exp(-h/(a*range));
Index: /issm/trunk-jpl/src/c/objects/Kriging/GaussianVariogram.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Kriging/GaussianVariogram.cpp	(revision 12366)
+++ /issm/trunk-jpl/src/c/objects/Kriging/GaussianVariogram.cpp	(revision 12367)
@@ -71,9 +71,10 @@
 	h2=pow(deltax,2.)+pow(deltay,2.);
 
-	/*return covariance*/
+	/*If h is too small, return sill*/
+	if(h2<0.0000001) return sill;
+
+	/*compute covariance*/
 	a     = 1./3.;
 	cova = (sill-nugget)*exp(-h2/(a*range*range));
-
-	if(h2<0.0000001) cova = sill;
 
 	return cova;