diff --git a/splinterpolator.h b/splinterpolator.h
index 8781cb2427b73c9ea093401e7724dce1fe88ba7d..49b258a128d345faf67c9591bed6e9db5fbe21f3 100644
--- a/splinterpolator.h
+++ b/splinterpolator.h
@@ -608,17 +608,6 @@ double Splinterpolator<T>::value_at(const double *coord) const
   ni = get_start_indicies(coord,inds);
   get_wgts(coord,inds,wgts);
 
-  if (coord[0] < -1 && coord[0] > -2 && coord[1] < 89 && coord[1] > 88 && coord[2] < 65 && coord[2] > 64) {
-    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]);
-      }
-    }
-  }
-
-
   double val=0.0;
   for (unsigned int m=0, me=(_ndim>4)?ni:1; m<me; m++) { 
     for (unsigned int l=0, le=(_ndim>3)?ni:1; l<le; l++) { 
@@ -630,10 +619,6 @@ 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] < -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, coef[%d] = %f\n",i,j,k,ii,wgt2,i,*iiwgt,ii,cptr[ii]);
-	    }
 	    val += cptr[linear2+indx2indx(inds[0]+i,0)]*(*iiwgt)*wgt2;
           }
 	}