Index: /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp	(revision 18914)
+++ /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp	(revision 18915)
@@ -628,16 +628,7 @@
     Observation observation=Observation(x_interp,y_interp,0.);
     std::vector<Observation> kNN;
-    
-    /*If radius is not provided or is 0, return all observations*/
-    if(radius==0.)
-    {
-        kNN=(this->covertree->getRoot())->getObservations();
-    }
-    else
-    {
-        kNN=(this->covertree->kNearestNeighbors(observation, maxdata));
-		//cout << "kNN's size: " << kNN.size() << endl;
-		
-    }
+
+	 kNN=(this->covertree->kNearestNeighbors(observation, maxdata));
+	 //cout << "kNN's size: " << kNN.size() << endl;
 	
 	//kNN is sort from closest to farthest neighbor
@@ -645,12 +636,14 @@
 	//deletes and resizes the kNN vector
 	vector<Observation>::iterator it;
-	for (it = kNN.begin(); it != kNN.end(); ++it) {
-		//(*it).print();
-		//cout << "\n" << (*it).distance(observation) << endl;
-		if ((*it).distance(observation) > radius) {
-			break;
-		}
-	}
-	kNN.erase(it, kNN.end());
+	if(radius>0.){
+		for (it = kNN.begin(); it != kNN.end(); ++it) {
+			//(*it).print();
+			//cout << "\n" << (*it).distance(observation) << endl;
+			if ((*it).distance(observation) > radius) {
+				break;
+			}
+		}
+		kNN.erase(it, kNN.end());
+	}
     
 	/*Allocate vectors*/
@@ -661,5 +654,5 @@
 	/*Loop over all observations and fill in x, y and obs*/
 	int i = 0;
-	for (it = kNN.begin(); it != kNN.end(); ++it) {
+	for(it = kNN.begin(); it != kNN.end(); ++it) {
 		(*it).WriteXYObs((*it), &x[i], &y[i], &obs[i]);
 		i++;
@@ -669,4 +662,4 @@
     *py=y;
     *pobs=obs;
-	*pnobs = kNN.size();
-}/*}}}*/
+	 *pnobs = kNN.size();
+}/*}}}*/