Changeset 38986 for trunk/Ohana/src/uniphot/src/update_catalog_setgalmodel.c

Timestamp:

Oct 27, 2015, 4:49:06 PM (11 years ago)

Author:

eugene

Message:

extensive work on relphot, relastro, uniphot, dvomerge aiming to the construction and calibration of PV3

Location:

trunk/Ohana

Files:

• . (modified) (1 prop)
• src/uniphot/src/update_catalog_setgalmodel.c (modified) (3 diffs)

trunk/Ohana/src/uniphot/src/update_catalog_setgalmodel.c

@@ -3 +3 @@
 int update_catalog_setgalmodel (Catalog *catalog) {
 
-  off_t i, j, m;
-
-  // choose this number more judiciously...
-  time_t timeRef = ohana_date_to_sec("2000/01/01,00:00:00");
+  off_t i, m;
 
   Average *average = catalog[0].average;
@@ -20 +17 @@
   for (i = 0; i < catalog[0].Naverage; i++) {
 
-    { 
-      m = average[i].measureOffset;
-      Measure *measure = &catalog[0].measure[m];
-      for (j = 0; j < average[i].Nmeasure; j++) {
-        measure[j].RoffGAL = 0.0; // XXX divide by dcos? project and calculate?
-        measure[j].DoffGAL = 0.0; // XXX divide by dcos? project and calculate?
-      }
-    }
+    average[i].uRgal = NAN;
+    average[i].uDgal = NAN;
 
     if (average[i].Nstarpar == 0) continue;
@@ -63 +54 @@
     starpar->uDEC = uD;
 
-    m = average[i].measureOffset;
-    Measure *measure = &catalog[0].measure[m];
-
-    // use the model uR,uD values to move the expected position based on the reference date;
-    for (j = 0; j < average[i].Nmeasure; j++) {
-
-      double dtime = (measure[j].t - timeRef) / 86400.0 / 365.25;
-
-      // do not modify the original Xccd,Yccd values
-      measure[j].RoffGAL = uR * dtime; // XXX divide by dcos? project and calculate?
-      measure[j].DoffGAL = uD * dtime; // XXX divide by dcos? project and calculate?
-    }
+    average[i].uRgal = uR;
+    average[i].uDgal = uD;
   }