schiff

commit d630e99171576f765361549d3ec09bcd193ace57
parent 419fa74ef25bdf88fd18d4de96c90e35f5f2620d
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Thu, 29 Oct 2015 08:15:10 +0100

Update the output of the schiff command

Add comments to the header and write the results as a list of double
precision floating point values, without any decoration. This should
make easier its treatment by other programs.

Diffstat:
M
src/schiff_sphere.c
 | 
15
++++++++++-----

1 file changed, 10 insertions(+), 5 deletions(-)
diff --git a/src/schiff_sphere.c b/src/schiff_sphere.c
@@ -224,7 +224,12 @@ schiff_run_sphere(const struct schiff_args* args)
   states = sa_add(states, nwlens);
   SSCHIFF(estimator_get_states(estimator, states));
 
-  fprintf(fp, "# Wavelength Extinction Absorption Scattering ProjArea\n");
+  fprintf(fp,
+"# Line format \"W E E' A A' S S' P P'\" with \"W\" the wavelength in meter,\n"
+"# \"E\", \"A\", and \"S\" the estimation of the extinction, absorption and\n"
+"# scattering cross section, respectively, and \"P\" the estimation of the\n"
+"# average projected area. The `prime' values, i.e. \"E'\", \"A'\", \"S'\" and \"P'\"\n"
+"# are the standard error of the aforementionned estimations.\n\n");
 
   /* Print the estimation results */
   FOR_EACH(iwlen, 0, nwlens) {
@@ -233,13 +238,13 @@ schiff_run_sphere(const struct schiff_args* args)
     fprintf(fp, "%g ", states[iwlen].wavelength);
 
     val = states[iwlen].values + SSCHIFF_EXTINCTION_CROSS_SECTION;
-    fprintf(fp, "%g +/- %g ", val->E, val->SE);
+    fprintf(fp, "%g %g ", val->E, val->SE);
     val = states[iwlen].values + SSCHIFF_ABSORPTION_CROSS_SECTION;
-    fprintf(fp, "%g +/- %g ", val->E, val->SE);
+    fprintf(fp, "%g %g ", val->E, val->SE);
     val = states[iwlen].values + SSCHIFF_SCATTERING_CROSS_SECTION;
-    fprintf(fp, "%g +/- %g ", val->E, val->SE);
+    fprintf(fp, "%g %g ", val->E, val->SE);
     val = states[iwlen].values + SSCHIFF_AVERAGE_PROJECTED_AREA;
-    fprintf(fp, "%g +/- %g\n", val->E, val->SE);
+    fprintf(fp, "%g %g\n", val->E, val->SE);
   }
 
 exit: