Index: trunk/psastro/src/psastroChooseRefstars.c
===================================================================
--- trunk/psastro/src/psastroChooseRefstars.c	(revision 37573)
+++ trunk/psastro/src/psastroChooseRefstars.c	(revision 41285)
@@ -46,4 +46,10 @@
 
     bool matchLumFunc = psMetadataLookupBool (&status, recipe, "PSASTRO.MATCH.LUMFUNC");
+
+    int nIter = psMetadataLookupS32 (&status, recipe, "PSASTRO.REFSTAR.CLUMP.NITER");
+    if (!status) nIter = 3;
+
+    psF32 clumpScale = psMetadataLookupS32 (&status, recipe, "PSASTRO.REFSTAR.CLUMP.SCALE");
+    if (!status) clumpScale = 150;
 
     pmFPAview *view = pmFPAviewAlloc (0);
@@ -178,5 +184,5 @@
 		// generate a reduced subset excluding the clumps
 		// XXX do we need both REFSTARS and SUBSET? 
-		psArray *subset = psastroRemoveClumpsIterate(refstars, 150, 3);
+		psArray *subset = psastroRemoveClumpsIterate(refstars, clumpScale, nIter);
 		psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.REFSTARS.SUBSET", PS_DATA_ARRAY, "astrometry objects", subset);