debug printouts for working on eddy

splinterpolator.h

@@ -606,11 +606,18 @@ double Splinterpolator<T>::value_at(const double *coord) const
   const T      *cptr = coef_ptr();
 
   ni = get_start_indicies(coord,inds);
+  get_wgts(coord,inds,wgts);
+
   if (coord[0] < -4 && coord[0] > -5 && coord[1] < 143 && coord[1] > 142 && coord[2] < 92 && coord[2] > 91) {
     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]);
+    for (int j=0; j<3; j++) {
+      for (int i=0; i<4; i++) {
+	printf("wgts[%d] = %f\n",j*4+i,wgts[j][i]);
+      }
+    }
   }
-  get_wgts(coord,inds,wgts);
+
 
   double val=0.0;
   for (unsigned int m=0, me=(_ndim>4)?ni:1; m<me; m++) { 
@@ -625,7 +632,7 @@ double Splinterpolator<T>::value_at(const double *coord) const
           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, 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;
           }