Changeset 35355

Timestamp:

Apr 3, 2013, 8:39:10 PM (13 years ago)

Author:

eugene

Message:

add output variables

File:

• branches/eam_branches/ipp-20130307/Ohana/src/opihi/cmd.astro/fitplx.c (modified) (2 diffs)

branches/eam_branches/ipp-20130307/Ohana/src/opihi/cmd.astro/fitplx.c

@@ -95 +95 @@
   Tmean /= Npts;
 
-  fprintf (stderr, "R,D : %f,%f, T: %f, Trange: %f, Tmin: %f, Tmax: %f\n", Rmean, Dmean, Tmean, Tmax - Tmin, Tmin, Tmax);
+  float Trange = Tmax - Tmin;
+  // fprintf (stderr, "R,D : %f,%f, T: %f, Trange: %f, Tmin: %f, Tmax: %f\n", Rmean, Dmean, Tmean, Trange, Tmin, Tmax);
 
   /* project coordinates to a plane centered on the object with units of arcsec */
@@ -141 +142 @@
   float parRange = hypot (dXRange, dYRange);
         
-  fprintf (stderr, "par factor range: %f\n", parRange);
+  // fprintf (stderr, "par factor range: %f\n", parRange);
 
   PlxFit fit;
   FitPMandPar (&fit, X, dX, Y, dY, t, pX, pY, n, VERBOSE);
 
-  fprintf (stderr, "Roff, Doff: %f, %f; dRo, dDo: %f, %f\n", fit.Ro, fit.Do, fit.dRo, fit.dDo);
+  // fprintf (stderr, "Roff, Doff: %f, %f; dRo, dDo: %f, %f\n", fit.Ro, fit.Do, fit.dRo, fit.dDo);
   
   XY_to_RD (&Rmean, &Dmean, fit.Ro, fit.Do, &coords);
-  fprintf (stderr, "Ro, Do: %f, %f\n", Rmean, Dmean);
-
-  fprintf (stderr, "uR, uD: %f, %f; duR, duD: %f, %f\n", fit.uR, fit.uD, fit.duR, fit.duD);
-  fprintf (stderr, "par: %f +/- %f\n", fit.p, fit.dp);
-
-  fprintf (stderr, "chisq: %f Nfit %d\n", fit.chisq, fit.Nfit);
+  if (VERBOSE) {
+    fprintf (stderr, "Ro, Do: %f, %f +/- %f, %f\n", Rmean, Dmean, fit.dRo, fit.dDo);
+    fprintf (stderr, "uR, uD: %f, %f; duR, duD: %f, %f\n", fit.uR, fit.uD, fit.duR, fit.duD);
+    fprintf (stderr, "par: %f +/- %f\n", fit.p, fit.dp);
+    fprintf (stderr, "chisq: %f Nfit %d\n", fit.chisq, fit.Nfit);
+  }
+
+  set_variable ("RA",   Rmean);
+  set_variable ("DEC",  Dmean);
+  set_variable ("dR",   fit.dRo);
+  set_variable ("dD",   fit.dDo);
+  set_variable ("uR",   fit.uR);
+  set_variable ("uD",   fit.uD);
+  set_variable ("duR",   fit.duR);
+  set_variable ("duD",   fit.duD);
+  set_variable ("plx",  fit.p);
+  set_variable ("dplx", fit.dp);
+  
+  set_variable ("Tmean",  Tmean);
+  set_variable ("Trange", Trange);
+  set_variable ("Prange", parRange);
+
+  set_variable ("chisq", fit.chisq);
+  set_variable ("Nfit",  fit.Nfit);
 
   return (TRUE);