Index: /branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro
===================================================================
--- /branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro	(revision 34054)
+++ /branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro	(revision 34055)
@@ -22,4 +22,36 @@
 # PSF_MODEL apply (GAUSS PS1_V1)
 # variance mode: CONSTANT, IMAGE_VAR, MODEL_VAR
+
+macro go.one
+  mkdir test
+
+  local sky fwhm psf_in psf_out
+
+  # mkref creates refimage.* and catdir.ref
+  file catdir.ref found
+  if (not($found))
+    mkref
+  end
+
+  $KAPA = kapa -noX
+  if (not($?RefConfig)) init
+
+  foreach sky 20.0
+    foreach fwhm 1.0
+      foreach psf_in PS1_V1
+
+        sprintf name "test/test.%02d.%02d.%s" $sky {10*$fwhm} $psf_in
+        mkexp $name $sky $fwhm $psf_in
+        
+        foreach psf_out GAUSS PS1_V1
+          foreach mode CONSTANT IMAGE_VAR MODEL_VAR
+            runphot $name $name.$psf_out.$mode "-D LINEAR_FIT_VARIANCE_MODE $mode -D PSF_MODEL PS_MODEL_$psf_out"
+            ckchip.mags $name.in.cmf $name.$psf_out.$mode.cmf $name.$psf_out.$mode 0.0
+          end
+        end
+      end
+    end
+  end
+end
 
 macro go
@@ -203,4 +235,8 @@
   style -pt 0 -sz 0.4
   show.pair 0
+  reindex ap = AP_MAG_out using index1
+  set dap = ap - v2
+  plot -c red -pt 2 -sz 0.5 rv dap
+
   delete -q imag_V dmag_V smag_V Dmag_V dmag_Vraw
   for imag -11 -6 0.2
@@ -807,5 +843,5 @@
   $RealConfig = $RealConfig -Db PSF.CONVOLVE F
   $RealConfig = $RealConfig -Df STARS.DENSITY 10.0
-  $RealConfig = $RealConfig -Df STARS.SIGMA.LIM 2.5
+  $RealConfig = $RealConfig -Df STARS.SIGMA.LIM 1.0
   $RealConfig = $RealConfig -Db GALAXY.FAKE F 
   $RealConfig = $RealConfig -Db GALAXY.GRID F