Index: trunk/Ohana/src/opihi/cmd.data/fit2d.c
===================================================================
--- trunk/Ohana/src/opihi/cmd.data/fit2d.c	(revision 5850)
+++ trunk/Ohana/src/opihi/cmd.data/fit2d.c	(revision 7917)
@@ -48,5 +48,5 @@
 
   if (argc != 5) {
-    fprintf (stderr, "USAGE: fit x y z order [-dz wt]\n");
+    gprint (GP_ERR, "USAGE: fit x y z order [-dz wt]\n");
     return (FALSE);
   }
@@ -57,13 +57,13 @@
 
   if (xvec[0].Nelements != yvec[0].Nelements) {
-    fprintf (stderr, "vectors must have same length\n");
+    gprint (GP_ERR, "vectors must have same length\n");
     return (FALSE);
   }
   if (xvec[0].Nelements != zvec[0].Nelements) {
-    fprintf (stderr, "vectors must have same length\n");
+    gprint (GP_ERR, "vectors must have same length\n");
     return (FALSE);
   }
   if (Weight && xvec[0].Nelements != dzvec[0].Nelements) {
-    fprintf (stderr, "vectors must have same length\n");
+    gprint (GP_ERR, "vectors must have same length\n");
     return (FALSE);
   }
@@ -160,20 +160,20 @@
     }
 
-    fprintf (stderr, "matrix:\n");
+    gprint (GP_ERR, "matrix:\n");
     for (k = 0; k < wterm; k++) {
 	for (K = 0; K < wterm; K++) {
-	    fprintf (stderr, "%f ", c[k][K]);
-	}
-	fprintf (stderr, "\n");
+	    gprint (GP_ERR, "%f ", c[k][K]);
+	}
+	gprint (GP_ERR, "\n");
     }
 
     gaussj (c, wterm, b, 1);
 
-    fprintf (stderr, "inverse:\n");
+    gprint (GP_ERR, "inverse:\n");
     for (k = 0; k < wterm; k++) {
 	for (K = 0; K < wterm; K++) {
-	    fprintf (stderr, "%f ", c[k][K]);
-	}
-	fprintf (stderr, "\n");
+	    gprint (GP_ERR, "%f ", c[k][K]);
+	}
+	gprint (GP_ERR, "\n");
     }
 
@@ -182,5 +182,5 @@
       for (i = ny = 0; ny < nterm; ny++) {
 	for (nx = 0; nx < nterm - ny; nx++, i++) {
-	  fprintf (stderr, "x^%d y^%d: %g\n", nx, ny, b[i][0]);
+	  gprint (GP_ERR, "x^%d y^%d: %g\n", nx, ny, b[i][0]);
 	}
       }
@@ -228,5 +228,5 @@
     maxsigma = ClipNSigma * sigma;
 
-    if (VERBOSE) fprintf (stderr, "mean: %g, sigma: %g, maxsigma: %g\n", mean, sigma, maxsigma);
+    if (VERBOSE) gprint (GP_ERR, "mean: %g, sigma: %g, maxsigma: %g\n", mean, sigma, maxsigma);
 
     /* mask outlier points */
@@ -245,13 +245,13 @@
       }	
     }
-    if (VERBOSE) fprintf (stderr, "pass: %d, Nmask: %d\n", N, Nmask);
-  }
-
-  if (!Quiet) fprintf (stderr, "z = ");
+    if (VERBOSE) gprint (GP_ERR, "pass: %d, Nmask: %d\n", N, Nmask);
+  }
+
+  if (!Quiet) gprint (GP_ERR, "z = ");
   for (N = ny = 0; ny < nterm; ny++) {
     for (nx = 0; nx < nterm - ny; nx++, N++) {
       sprintf (name, "CX%dY%d", nx, ny);
       set_variable (name, b[N][0]);
-      if (!Quiet) fprintf (stderr, "%f x^%d y^%d  ", b[N][0], nx, ny);
+      if (!Quiet) gprint (GP_ERR, "%f x^%d y^%d  ", b[N][0], nx, ny);
     }
   }
@@ -259,5 +259,5 @@
   set_variable (name, (double) order);
   
-  if (!Quiet) fprintf (stderr, "\n");
+  if (!Quiet) gprint (GP_ERR, "\n");
 
   /* free internal data */