debug printouts for working on eddy

splinterpolator.h

@@ -632,7 +632,7 @@ double Splinterpolator<T>::value_at(const double *coord) const
           for (unsigned int i=0; i<ni; i++, iiwgt++) {
 	    if (coord[0] < -1 && coord[0] > -2 && coord[1] < 89 && coord[1] > 88 && coord[2] < 65 && coord[2] > 64) {
 	      int ii = linear2+indx2indx(inds[0]+i,0);
-              // printf("i = %d, j = %d, k = %d, ii = %d, wgt2 = %f, wgts[%d] = %f\n",i,j,k,ii,wgt2,i,*iiwgt);
+              printf("i = %d, j = %d, k = %d, ii = %d, wgt2 = %f, wgts[%d] = %f\n",i,j,k,ii,wgt2,i,*iiwgt);
 	    }
 	    val += cptr[linear2+indx2indx(inds[0]+i,0)]*(*iiwgt)*wgt2;
           }