Changeset 13779

Timestamp:

Jun 12, 2007, 6:00:35 PM (19 years ago)

Author:

eugene

Message:

adding FRNG values to output header; fixing error scaling

File:

• trunk/ppImage/src/ppImageDetrendFringe.c (modified) (2 diffs)

trunk/ppImage/src/ppImageDetrendFringe.c

@@ -45 +45 @@
         expTime = 1.0;
     }
+    // normalize by exposure time : NOTE df is 1/sigma
     psBinaryOp(measurements->f, measurements->f, "*", psScalarAlloc(1.0 / expTime, PS_TYPE_F32));
-    psBinaryOp(measurements->df, measurements->df, "*", psScalarAlloc(1.0 / expTime, PS_TYPE_F32));
+    psBinaryOp(measurements->df, measurements->df, "*", psScalarAlloc(expTime, PS_TYPE_F32));
 
     // Science fringe measurements
@@ -189 +190 @@
                          md5string, "");
         psFree(md5string);
+
+        // write metadata header value
+        psString keyword = NULL;
+        psStringAppend (&keyword, "FRNG_%02d", i);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, keyword, PS_META_REPLACE, "Fringe Amplitude", solution->coeff->data.F32[i + 1]);
+        psFree (keyword);
+
+        keyword = NULL;
+        psStringAppend (&keyword, "FRNG_%02dD", i);
+        psMetadataAddF32(hdu->header, PS_LIST_TAIL, keyword, PS_META_REPLACE, "Fringe Amplitude error", solution->coeffErr->data.F32[i + 1]);
+        psFree (keyword);
     }
     if (expTime != 1.0) {