debug printouts for working on eddy

splinterpolator.h

@@ -102,7 +102,7 @@ public:
   }
   T ValAndDerivs(double x, double y, double z, std::vector<T>& rderiv) const;
 
-  // Return continous derivative at voxel centres (only works for order<1)
+  // Return continous derivative at voxel centres (only works for order>1)
   T Deriv(const std::vector<unsigned int>& indx, unsigned int ddir) const;
   T Deriv1(const std::vector<unsigned int>& indx) const {return(Deriv(indx,0));}
   T Deriv2(const std::vector<unsigned int>& indx) const {return(Deriv(indx,1));}
@@ -606,6 +606,7 @@ double Splinterpolator<T>::value_at(const double *coord) const
   const T      *cptr = coef_ptr();
 
   ni = get_start_indicies(coord,inds);
+  for (unsigned int i=0; i<ni; i++) printf("inds[%d] = %d\n",i,inds[i]);
   get_wgts(coord,inds,wgts);
 
   double val=0.0;