Index: trunk/ippScripts/Build.PL
===================================================================
--- trunk/ippScripts/Build.PL	(revision 37052)
+++ trunk/ippScripts/Build.PL	(revision 37053)
@@ -135,4 +135,5 @@
         scripts/psphot_fullforce_warp.pl
         scripts/psphot_fullforce_summary.pl
+        scripts/stack_bkg_mk_mdc.pl
         scripts/sc_mk_stack_mdc.pl
         scripts/sc_prepare_camera.pl
Index: trunk/ippScripts/scripts/psphot_fullforce_warp.pl
===================================================================
--- trunk/ippScripts/scripts/psphot_fullforce_warp.pl	(revision 37052)
+++ trunk/ippScripts/scripts/psphot_fullforce_warp.pl	(revision 37053)
@@ -186,4 +186,5 @@
     $command .= " -F PSPHOT.OUTPUT PSPHOT.OUT.CMF.MEF";
     $command .= " -F PSPHOT.BACKMDL PSPHOT.BACKMDL.MEF";
+    $command .= " -photcode-rule '{DETECTOR}.{FILTER.ID}.ForcedWarp'";
     $command .= " -tracedest $traceDest -log $logDest";
     $command .= " -dbname $dbname" if defined $dbname;
Index: trunk/psastro/doc/outline.txt
===================================================================
--- trunk/psastro/doc/outline.txt	(revision 37052)
+++ trunk/psastro/doc/outline.txt	(revision 37053)
@@ -3,9 +3,12 @@
 
   psastroParseCamera
+    psastroDefineFiles
   psastroDataLoad
-    psastroConvertReadout (save PSASTRO.RAWSTARS & PSASTRO.RAWSTARS.SUBSET on readout->analysis)
+    psastroConvertReadout (saves PSASTRO.RAWSTARS on readout->analysis)
+      psastroCorrectKH
   psastroAnalysis
     psastroUseModel
     psastroAstromGuess (plots PSASTRO.RAWSTARS)
+    psastroRemoveClumpsRawstars (saves PSASTRO.RAWSTARS.SUBSET on readout->analysis)
     psastroLoadRefstars (loads full reference star set)
     psastroChooseRefstars (save PSASTRO.REFSTARS & PSASTRO.REFSTARS.SUBSET on readout->analysis)
Index: trunk/psphot/test/tap_psphot_galaxies.pro
===================================================================
--- trunk/psphot/test/tap_psphot_galaxies.pro	(revision 37052)
+++ trunk/psphot/test/tap_psphot_galaxies.pro	(revision 37053)
@@ -32,4 +32,5 @@
 $FakeOptions = $BaseOptions
 $FakeOptions = $FakeOptions -exptime 30.0
+# $FakeOptions = $FakeOptions -D PSF.MODEL PS_MODEL_PS1_V1
 $FakeOptions = $FakeOptions -D PSF.MODEL PS_MODEL_GAUSS
   
@@ -44,27 +45,56 @@
 
 # basic config for ppSim with randomly distributed stars and gridded galaxies
-$FakeConfig = -camera SIMTEST
-$FakeConfig = $FakeConfig -recipe PPSIM STACKTEST.RUN
-$FakeConfig = $FakeConfig -D PSASTRO:PSASTRO.CATDIR catdir.ref
-$FakeConfig = $FakeConfig -Db STARS.FAKE F                         ; # only use stars from catdir.ref
-$FakeConfig = $FakeConfig -Db MATCH.DENSITY F
-$FakeConfig = $FakeConfig -Db PSF.CONVOLVE T
-$FakeConfig = $FakeConfig -Db GALAXY.FAKE T                        ; # generate a "realistic" distribution of galaxies
-$FakeConfig = $FakeConfig -Df GALAXY.MAG 17.0
-$FakeConfig = $FakeConfig -Db GALAXY.GRID T                        ; # generate a grid of galaxies (constant mag)
-#$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_GAUSS
-#$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_EXP
-#$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_SERSIC
-#$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_DEV
-$FakeConfig = $FakeConfig -Df GALAXY.RMAJOR.MIN 10.0
-$FakeConfig = $FakeConfig -Df GALAXY.RMAJOR.MAX 10.0
-$FakeConfig = $FakeConfig -Df GALAXY.ARATIO.MIN 0.25
-$FakeConfig = $FakeConfig -Df GALAXY.ARATIO.MAX 0.25
-$FakeConfig = $FakeConfig -Df GALAXY.THETA.MIN 0 
-$FakeConfig = $FakeConfig -Df GALAXY.THETA.MAX 180
-$FakeConfig = $FakeConfig -Df GALAXY.INDEX.MIN 1.66
-$FakeConfig = $FakeConfig -Df GALAXY.INDEX.MAX 1.66
-$FakeConfig = $FakeConfig -Di GALAXY.GRID.DX 120
-$FakeConfig = $FakeConfig -Di GALAXY.GRID.DY 120
+if (1)
+  $FakeConfig = -camera SIMTEST
+  $FakeConfig = $FakeConfig -recipe PPSIM STACKTEST.RUN
+  $FakeConfig = $FakeConfig -D PSASTRO:PSASTRO.CATDIR catdir.ref
+  $FakeConfig = $FakeConfig -Db STARS.FAKE F                         ; # only use stars from catdir.ref
+  $FakeConfig = $FakeConfig -Db MATCH.DENSITY F
+  $FakeConfig = $FakeConfig -Db PSF.CONVOLVE T
+  $FakeConfig = $FakeConfig -Db GALAXY.FAKE T                        ; # generate a "realistic" distribution of galaxies
+  $FakeConfig = $FakeConfig -Df GALAXY.MAG 17.0
+  $FakeConfig = $FakeConfig -Db GALAXY.GRID F                        ; # generate a grid of galaxies (constant mag)
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_GAUSS
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_EXP
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_SERSIC
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_DEV
+  $FakeConfig = $FakeConfig -Df GALAXY.RMAJOR.MIN 10.0
+  $FakeConfig = $FakeConfig -Df GALAXY.RMAJOR.MAX  1.0
+  $FakeConfig = $FakeConfig -Df GALAXY.ARATIO.MIN 0.25
+  $FakeConfig = $FakeConfig -Df GALAXY.ARATIO.MAX 1.00
+  $FakeConfig = $FakeConfig -Df GALAXY.THETA.MIN 0 
+  $FakeConfig = $FakeConfig -Df GALAXY.THETA.MAX 180
+  $FakeConfig = $FakeConfig -Df GALAXY.INDEX.MIN 1.66
+  $FakeConfig = $FakeConfig -Df GALAXY.INDEX.MAX 1.66
+  $FakeConfig = $FakeConfig -Di GALAXY.GRID.DX 120
+  $FakeConfig = $FakeConfig -Di GALAXY.GRID.DY 120
+end
+
+# basic config for ppSim with randomly distributed stars and gridded galaxies
+if (0) 
+  $FakeConfig = -camera SIMTEST
+  $FakeConfig = $FakeConfig -recipe PPSIM STACKTEST.RUN
+  $FakeConfig = $FakeConfig -D PSASTRO:PSASTRO.CATDIR catdir.ref
+  $FakeConfig = $FakeConfig -Db STARS.FAKE F                         ; # only use stars from catdir.ref
+  $FakeConfig = $FakeConfig -Db MATCH.DENSITY F
+  $FakeConfig = $FakeConfig -Db PSF.CONVOLVE T
+  $FakeConfig = $FakeConfig -Db GALAXY.FAKE T                        ; # generate a "realistic" distribution of galaxies
+  $FakeConfig = $FakeConfig -Df GALAXY.MAG 17.0
+  $FakeConfig = $FakeConfig -Db GALAXY.GRID T                        ; # generate a grid of galaxies (constant mag)
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_GAUSS
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_EXP
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_SERSIC
+  #$FakeConfig = $FakeConfig -D GALAXY.MODEL PS_MODEL_DEV
+  $FakeConfig = $FakeConfig -Df GALAXY.RMAJOR.MIN 10.0
+  $FakeConfig = $FakeConfig -Df GALAXY.RMAJOR.MAX 10.0
+  $FakeConfig = $FakeConfig -Df GALAXY.ARATIO.MIN 0.25
+  $FakeConfig = $FakeConfig -Df GALAXY.ARATIO.MAX 0.25
+  $FakeConfig = $FakeConfig -Df GALAXY.THETA.MIN 0 
+  $FakeConfig = $FakeConfig -Df GALAXY.THETA.MAX 180
+  $FakeConfig = $FakeConfig -Df GALAXY.INDEX.MIN 1.66
+  $FakeConfig = $FakeConfig -Df GALAXY.INDEX.MAX 1.66
+  $FakeConfig = $FakeConfig -Di GALAXY.GRID.DX 120
+  $FakeConfig = $FakeConfig -Di GALAXY.GRID.DY 120
+end
 
 list fwhm 
@@ -164,5 +194,10 @@
   $psphotConfig = $psphotConfig -Di PSPHOT:LMM_FIT_GAIN_FACTOR_MODE 2
   $psphotConfig = $psphotConfig -Db PSPHOT:SAVE.RESID T
+  $psphotConfig = $psphotConfig -D  PSPHOT:PSF_MODEL PS_MODEL_PS1_V1
+  # $psphotConfig = $psphotConfig -D  PSPHOT:PSF_MODEL PS_MODEL_GAUSS
   $psphotConfig = $psphotConfig -D  PSPHOT:EXTENDED_SOURCE_MODELS_SELECTION $fitModel
+  # $psphotConfig = $psphotConfig -D  PSPHOT:OUTPUT.FORMAT PS1_V5
+  $psphotConfig = $psphotConfig -D  PSPHOT:OUTPUT.FORMAT PS1_SV3
+  $psphotConfig = $psphotConfig -Db PSPHOT:LENSING_PARAMETERS T
 
   # ppImage / psphot on the output
@@ -217,4 +252,16 @@
 
   # plot (input - output) vs mag
+end
+
+macro ck.poserrors
+  reindex X_raw_m = X_raw using index2
+  reindex X_out_m = X_PSF_out using index1
+  reindex dX_out_m = X_PSF_SIG_out using index1
+  set Xoff = X_raw_m - X_out_m
+  reindex M_raw_m = PSF_INST_MAG_raw using index2
+
+  # lim M_raw_m Xoff; clear; box; plot M_raw_m Xoff
+  set dC = Xoff^2 / dX_out_m^2
+  lim M_raw_m dX_out_m; clear; box; plot M_raw_m dX_out_m
 end
 
Index: trunk/pswarp/src/pswarpCleanup.c
===================================================================
--- trunk/pswarp/src/pswarpCleanup.c	(revision 37052)
+++ trunk/pswarp/src/pswarpCleanup.c	(revision 37053)
@@ -55,4 +55,5 @@
 	goto DONE;
     }
+
     pmChip *chip;
     while ((chip = pmFPAviewNextChip (view, output->fpa, 1)) != NULL) {