Changeset 41608 for branches/eam_branches/relphot.20210521/src/GridOps.c

Timestamp:

May 22, 2021, 1:44:24 PM (5 years ago)

Author:

eugene

Message:

fix inconsistency between Mgrid and Mflat

File:

• branches/eam_branches/relphot.20210521/src/GridOps.c (modified) (3 diffs)

branches/eam_branches/relphot.20210521/src/GridOps.c

@@ -188 +188 @@
         if (isnan(Mmos)) continue;
 
+        float Mflat = getMflat (nm, nc, catalog);
+        if (isnan(Mflat)) continue;
+
         // Mrel* is the average magnitude for this star.  For PS1 stacks, we have too much
         // PSF variability.  We need to calibrate the PSF magnitudes separately from the
@@ -200 +203 @@
         if (isnan(MsysPSF)) continue;
 
-        // what about Mflat?
-        float Moff =  Mcal + Mgrp + Mmos;
+        float Moff =  Mcal + Mgrp + Mmos + Mflat;
 
         // Msys = Mrel + Moff + Mgrid
@@ -307 +309 @@
   for (off_t j = 0; j < catalog->Nmeasure; j++) {
     float Mgrid = getMgrid (&catalog->measure[j]);
-    catalog->measure[j].Mflat -= Mgrid;
-  }
-}
-
+    catalog->measure[j].Mflat += Mgrid;
+  }
+}
+