Index: trunk/psastro/src/psastroRemoveClumps.c
===================================================================
--- trunk/psastro/src/psastroRemoveClumps.c	(revision 31333)
+++ trunk/psastro/src/psastroRemoveClumps.c	(revision 41285)
@@ -35,4 +35,10 @@
     }
 
+    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);
     pmFPA *fpa = input->fpa;
@@ -58,5 +64,5 @@
 		// XXX do we need both RAWSTARS and SUBSET? 
 		// XXX put these parameters in the recipe, please
-		psArray *subset = psastroRemoveClumpsIterate(rawstars, 150, 3);
+		psArray *subset = psastroRemoveClumpsIterate(rawstars, clumpScale, nIter);
 		psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.RAWSTARS.SUBSET", PS_DATA_ARRAY, "astrometry objects", subset);
 		psFree (subset);
@@ -65,5 +71,5 @@
                 if ((gridstars == rawstars) || (gridstars == NULL)) { continue; }
 
-		psArray *gridstars_subset = psastroRemoveClumpsIterate(gridstars, 150, 3);
+		psArray *gridstars_subset = psastroRemoveClumpsIterate(gridstars, clumpScale, nIter);
 		psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PSASTRO.GRID.RAWSTARS.SUBSET", PS_DATA_ARRAY, "astrometry objects", gridstars_subset);
 		psFree (gridstars_subset);
@@ -79,5 +85,5 @@
     psArray *newset = psMemIncrRefCounter (input);
     for (int i = 0; i < nIter; i++) {
-	psArray *subset = psastroRemoveClumps (newset, 150);
+	psArray *subset = psastroRemoveClumps (newset, scale);
 	psFree (newset);
 	newset = subset;