Index: /branches/czw_branch/20160809/Ohana/src/fakeastro/src/args.c
===================================================================
--- /branches/czw_branch/20160809/Ohana/src/fakeastro/src/args.c	(revision 39747)
+++ /branches/czw_branch/20160809/Ohana/src/fakeastro/src/args.c	(revision 39748)
@@ -128,5 +128,5 @@
   }
 
-  MAX_MAG_GAIA = 16.0;
+  MAX_MAG_GAIA = 21.0;
   if ((N = get_argument (*argc, argv, "-gaia-limit"))) {
     remove_argument (N, argc, argv);
Index: /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_fake_stars_catalog.c
===================================================================
--- /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_fake_stars_catalog.c	(revision 39747)
+++ /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_fake_stars_catalog.c	(revision 39748)
@@ -117,5 +117,5 @@
 
     // XXX TEST
-    dDsee = dRsee = 0.0;
+    // dDsee = dRsee = 0.0;
 
     double dRoff = (dRpm + dRsee) / 3600.0;
Index: /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_fakestars.c
===================================================================
--- /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_fakestars.c	(revision 39747)
+++ /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_fakestars.c	(revision 39748)
@@ -35,4 +35,6 @@
       // we can generate a distribution which is uniform on the sky, in which case
       // we can limit to the selected patch analyically
+
+      // note: r & z are generated in parsec
       r = pow(drand48(), 0.33333) * FAKEASTRO_RGAL;
       z = 0.0;
@@ -48,4 +50,5 @@
       int inPatch = FALSE;
       while (!inPatch) {
+	// note: r & z are generated in parsec
 	z = ohana_gaussdev_rnd (0.0, FAKEASTRO_ZGAL);
 	r = sqrt(drand48()) * FAKEASTRO_RGAL;
@@ -68,4 +71,5 @@
     // double y = r*sin(L);
 
+    // distance here is in parsec
     double distance = sqrt (SQ(r) + SQ(z));
 
@@ -74,5 +78,6 @@
 
     double uL_sol, uB_sol;
-    SolarMotionModel_radians(&uL_sol, &uB_sol, Lrad, Brad, distance);
+    SolarMotionModel_radians(&uL_sol, &uB_sol, Lrad, Brad, distance / 1000.0);
+    // note: SolarMotionModel wants distance in kpc
 
     double uL = uL_gal + uL_sol;
Index: /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_gaia_measures.c
===================================================================
--- /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_gaia_measures.c	(revision 39747)
+++ /branches/czw_branch/20160809/Ohana/src/fakeastro/src/make_gaia_measures.c	(revision 39748)
@@ -42,5 +42,5 @@
 
     // make a crude G-i color for now:
-    double G_PS1 = secfilt[i*Nsecfilt + Nsec].M + 1.0; // make all stars have G-i = 1.0
+    double G_PS1 = secfilt[i*Nsecfilt + Nsec].M; // make all stars have G-i = 1.0
     if (G_PS1 > MAX_MAG_GAIA) continue; // only generate GAIA detections for objects with G_PS1 < 21.0
 
@@ -78,9 +78,9 @@
 
     // uR,uD in linear arcsec
-    double dRsee = ohana_gaussdev_rnd (0.0, 1.0 / SN);
-    double dDsee = ohana_gaussdev_rnd (0.0, 1.0 / SN);
+    double dRsee = ohana_gaussdev_rnd (0.0, 0.01 / SN);
+    double dDsee = ohana_gaussdev_rnd (0.0, 0.01 / SN);
 
     // XXX TEST
-    dDsee = dRsee = 0.0;
+    // dDsee = dRsee = 0.0;
 
     double dRoff = (dRpm + dRsee) / 3600.0;