debug printouts for working on eddy

splinterpolator.h

@@ -606,7 +606,10 @@ 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]);
+  if (coord[0] < 41 && coord[0] > 40 && coord[1] < 71 && coord[1] > 70 && coord[2] < 91 && coord[2] > 90) {
+    for (unsigned int i=0; i<ni; i++) printf("coord[%d] = %f\n",i,coord[i]);
+    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;