Index: /issm/trunk-jpl/src/c/classes/kriging/Covertree.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/kriging/Covertree.cpp	(revision 18918)
+++ /issm/trunk-jpl/src/c/classes/kriging/Covertree.cpp	(revision 18919)
@@ -6,9 +6,9 @@
 
 	/*Constructors/Destructors*/
-Covertree::Covertree(const double& maxDist,const std::vector<Observation>& points){
+Covertree::Covertree(int maxLevel,const std::vector<Observation>& points){
 	this->base = 2.;
 	_root=NULL;
 	_numNodes=0;
-	_maxLevel=ceilf(log(maxDist)/log(base));
+	_maxLevel=maxLevel;//ceilf(log(maxDist)/log(base));
 	_minLevel=_maxLevel-1;
 	std::vector<Observation>::const_iterator it;
Index: /issm/trunk-jpl/src/c/classes/kriging/Covertree.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/kriging/Covertree.h	(revision 18918)
+++ /issm/trunk-jpl/src/c/classes/kriging/Covertree.h	(revision 18919)
@@ -89,5 +89,5 @@
 	 */
 
-	Covertree(const double& maxDist,
+	Covertree(int maxDist,
 				const std::vector<Observation>& points=std::vector<Observation>()); 
 	~Covertree();
Index: /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp	(revision 18918)
+++ /issm/trunk-jpl/src/c/classes/kriging/Observations.cpp	(revision 18919)
@@ -157,13 +157,11 @@
 void Observations::InitCovertree(IssmPDouble* observations_list,IssmPDouble* x,IssmPDouble* y,int n,Options* options){/*{{{*/
 
-    int maxdepth = 10;
+    int maxdepth = 20;
 	 _printf0_("Generating covertree with a maximum depth " <<  maxdepth <<"... ");
     this->covertree=new Covertree(maxdepth);
 
     for(int i=0;i<n;i++){
-		 if(i%1000) printf("progress = %g \n",double(i)/double(n)*100.);
 		 this->covertree->insert(Observation(x[i],y[i],observations_list[i]));
     }
-
 	 _printf0_("done\n");
 }
@@ -619,5 +617,4 @@
 }/*}}}*/
 
-
 void Observations::ObservationList2(double **px,double **py,double **pobs,int* pnobs,double x_interp,double y_interp,double radius,int maxdata){/*{{{*/
     
@@ -630,5 +627,5 @@
 
 	 kNN=(this->covertree->kNearestNeighbors(observation, maxdata));
-	 //cout << "kNN's size: " << kNN.size() << endl;
+	 //cout << "kNN's size: " << kNN.size() << " (maxdata = " <<maxdata<<")"<<endl;
 	
 	//kNN is sort from closest to farthest neighbor