debug printouts for working on eddy

splinterpolator.h

@@ -623,6 +623,10 @@ double Splinterpolator<T>::value_at(const double *coord) const
           int linear2 = add2linear(linear1,inds[1]+j);
           double *iiwgt=iwgt;
           for (unsigned int i=0; i<ni; i++, iiwgt++) {
+	    if (coord[0] < -4 && coord[0] > -5 && coord[1] < 143 && coord[1] > 142 && coord[2] < 92 && coord[2] > 91) {
+	      int ii = linear2+indx2indx(inds[0]+i,0);
+              printf("i = %d, j = %d, k = %d, ii = %d\n",i,j,k,ii);
+	    }
 	    val += cptr[linear2+indx2indx(inds[0]+i,0)]*(*iiwgt)*wgt2;
           }
 	}