Index: /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/ConfigInit.c
===================================================================
--- /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/ConfigInit.c	(revision 37568)
+++ /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/ConfigInit.c	(revision 37569)
@@ -55,8 +55,8 @@
   }
   if (!ScanConfig (config, "FAKEASTRO_REF_EPOCH", "%s", 0, FAKEASTRO_REF_EPOCH)) {
-    strcpy (FAKEASTRO_REF_EPOCH, "2011/05/11,00:00:00"); // epoch of reference catalog
+    strcpy (FAKEASTRO_REF_EPOCH, "2000/01/01,00:00:00"); // epoch of truth positions
   }
   if (!ScanConfig (config, "FAKEASTRO_2MASS_EPOCH", "%s", 0, FAKEASTRO_2MASS_EPOCH)) {
-    strcpy (FAKEASTRO_2MASS_EPOCH, "2000/01/01,00:00:00"); // epoch of reference catalog
+    strcpy (FAKEASTRO_2MASS_EPOCH, "2000/01/01,00:00:00"); // epoch of 2MASS astrometry
   }
 
Index: /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_2mass_measures.c
===================================================================
--- /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_2mass_measures.c	(revision 37568)
+++ /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_2mass_measures.c	(revision 37569)
@@ -11,4 +11,6 @@
   Measure *measure = catalog->measure;
   
+  // we have starpar.R,D, which represent the true positions at the reference epoch, FAKEASTRO_REF_EPOCH
+  // the 2MASS stars are generated at locations for the 2MASS EPOCH
   time_t timeRef = ohana_date_to_sec(FAKEASTRO_REF_EPOCH);
   time_t tzero_2mass = ohana_date_to_sec(FAKEASTRO_2MASS_EPOCH);
Index: /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fake_stars_catalog.c
===================================================================
--- /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fake_stars_catalog.c	(revision 37568)
+++ /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fake_stars_catalog.c	(revision 37569)
@@ -20,4 +20,5 @@
   int Nstars = *nstars;
   
+  // we have starpar.R,D, which represent the true positions at the reference epoch, FAKEASTRO_REF_EPOCH
   time_t timeRef = ohana_date_to_sec(FAKEASTRO_REF_EPOCH);
   
Index: /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fakeqsos.c
===================================================================
--- /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fakeqsos.c	(revision 37568)
+++ /branches/eam_branches/ipp-20140904/Ohana/src/fakeastro/src/make_fakeqsos.c	(revision 37569)
@@ -14,4 +14,17 @@
 
   CoordTransform *transback = InitTransform (COORD_CELESTIAL, COORD_GALACTIC);
+
+  FILE *f = NULL;
+  if (ICRF) {
+    f = fopen ("fake.icrf.dat", "w");
+    if (!f) {
+      fprintf (stderr, "unable to open fake.icrf.dat\n");
+    }
+  } else {
+    f = fopen ("fake.zero.dat", "w");
+    if (!f) {
+      fprintf (stderr, "unable to open fake.zero.dat\n");
+    }
+  } 
 
   int i;
@@ -68,4 +81,41 @@
     stars[i].starpar.uRA  = 0.0;
     stars[i].starpar.uDEC = 0.0;
+
+/*
+#                       hr mn seconds   deg mn seconds      mas    mas                   Jy     Jy      Jy     Jy      Jy     Jy     Jy     Jy      Jy     Jy
+C  2357+080 J0000+0816  00 00 07.031141 +08 16 45.05175    0.46   0.85   0.758     41  -1.00  -1.00   -1.00  -1.00    0.020 <0.014  -1.00  -1.00   -1.00  -1.00   X    rfc_2014c
+01234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789
+0         1         2         3         4         5         6         7         8         9         0         1         2         3         4         5         6         
+*/
+
+    if (f) {
+      char Rline[32], Dline[32];
+      {
+	int flag = SIGN(R);
+	double hours = fabs(R / 15.0);  /* convert from degrees to hours */
+	int h = hours;
+	int m = 60.000001*(hours - h);
+	float s = 3600*(hours - h - m / 60.0);
+	if (flag > 0) {
+	  snprintf (Rline, 32, " %02d %02d %09.6f  ", h, m, s);
+	} else {
+	  snprintf (Rline, 32, "-%02d %02d %09.6f  ", h, m, s);
+	}	
+      }
+      {
+	int flag = SIGN(D);
+	double dec = fabs(D);
+	int d = dec;
+	int m = 60.000001*(dec - d);
+	float s = 3600*(dec - d - m / 60.0);
+	if (flag > 0) {
+	  snprintf (Dline, 32, " %02d %02d %09.6f  ", d, m, s);
+	} else {
+	  snprintf (Dline, 32, "-%02d %02d %09.6f  ", d, m, s);
+	}	
+      }
+
+      fprintf (f, "C  2357+080 J0000+0816 %16s %15s    0.46   0.85   0.758     41  -1.00  -1.00   -1.00  -1.00    0.020 <0.014  -1.00  -1.00   -1.00  -1.00   X    rfc_2014c\n", Rline, Dline);
+    }
   }
   fprintf (stderr, "\n");