Index: /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/include/fakeastro.h
===================================================================
--- /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/include/fakeastro.h	(revision 38706)
+++ /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/include/fakeastro.h	(revision 38707)
@@ -103,4 +103,5 @@
 
 float  RADIUS;
+float  MAX_MAG_2MASS;
 
 SkyRegion UserPatch;
Index: /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/args.c
===================================================================
--- /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/args.c	(revision 38706)
+++ /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/args.c	(revision 38707)
@@ -104,4 +104,11 @@
   if ((N = get_argument (*argc, argv, "-uniform-radec"))) {
     UNIFORM_RADEC = TRUE;
+    remove_argument (N, argc, argv);
+  }
+
+  MAX_MAG_2MASS = 16.0;
+  if ((N = get_argument (*argc, argv, "-2mass-limit"))) {
+    remove_argument (N, argc, argv);
+    MAX_MAG_2MASS = atof(argv[N]);
     remove_argument (N, argc, argv);
   }
Index: /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/make_2mass_measures.c
===================================================================
--- /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/make_2mass_measures.c	(revision 38706)
+++ /branches/eam_branches/ipp-20150625/Ohana/src/fakeastro/src/make_2mass_measures.c	(revision 38707)
@@ -24,6 +24,6 @@
   // PhotCode *codeK = GetPhotcodebyName ("2MASS_K");
 
-  float ZP = 21.0; // this is chosen to give 16.0 mag stars a flux of 100.0 counts (crude, yes)
-  float SkyCts = 300.0; // this is chosen to make SN @ 16.0 = 5.0 (
+  float ZP = MAX_MAG_2MASS + 6.5; // this is chosen to give stars at the limit a flux of 400.0 counts (crude, yes)
+  float SkyCts = 1000.0; // this is chosen to make SN @ limit ~ 10.0 (
 
   int Nmeasure = catalog->Nmeasure;
@@ -45,5 +45,5 @@
     // make a crude JHK color for now:
     double J_PS1 = secfilt[i*Nsecfilt + Nsec].M - 1.0; // make all stars have y-J = 1.0
-    if (J_PS1 > 16.0) continue; // only generate 2MASS detections for objects with J_PS1 < 16.0
+    if (J_PS1 > MAX_MAG_2MASS) continue; // only generate 2MASS detections for objects with J_PS1 < 16.0
 
     if (isnan(secfilt[i*Nsecfilt + Nsec].M)) {
