Index: /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp	(revision 17326)
+++ /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp	(revision 17327)
@@ -202,5 +202,5 @@
 
 	/*If radius is not provided or is 0, return all observations*/
-	if(radius==0.) radius=this->quadtree->root->length;
+	if(radius==0.) radius=this->quadtree->root->length*2.;
 
 	/*Compute radius square*/
@@ -466,7 +466,4 @@
 	_assert_(pprediction);
 
-	/*If radius is not provided or is 0, return all observations*/
-	if(radius==0) radius=this->quadtree->root->length;
-
 	/*Get list of observations for current point*/
 	this->ObservationList(&x,&y,&obs,&n_obs,x_interp,y_interp,radius,maxdata);
@@ -505,4 +502,6 @@
 				Green[i*n_obs+j] = Green[j*n_obs+i];
 			}
+			/*Zero diagonal (should be done already, but just in case)*/
+			Green[i*n_obs+i] = 0.;
 		}