Index: /branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro
===================================================================
--- /branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro	(revision 34055)
+++ /branches/eam_branches/ipp-20120601/psphot/test/tap_psphot_varmodel.pro	(revision 34056)
@@ -23,4 +23,14 @@
 # variance mode: CONSTANT, IMAGE_VAR, MODEL_VAR
 
+macro ppsimtest
+  data $1
+  # sf1 : measured flux inserted in image
+  # sm2 : theoretical magnitude for star
+  read x 1 y 2 sf1 3 sm2 5 Io 12
+  set sm1 = -2.5*log(sf1)
+  set dsm = sm1 - sm2
+  lim -n 2 sm1 -0.2 0.2; clear; box; plot sm1 dsm
+end
+
 macro go.one
   mkdir test
@@ -43,4 +53,8 @@
         sprintf name "test/test.%02d.%02d.%s" $sky {10*$fwhm} $psf_in
         mkexp $name $sky $fwhm $psf_in
+	data $name.dat
+	read fluxIn 3 MagPredIn 5
+	set MagRealIn = -2.5*log(fluxIn)
+	set dMagIn = MagRealIn - MagPredIn
         
         foreach psf_out GAUSS PS1_V1
@@ -256,4 +270,5 @@
   section a1 0.0 0.50 1.0 0.25
   lim rv -0.025 0.075; label -fn courier 14; box; 
+  plot -c black -pt 0 -sz 0.3 MagPredIn dMagIn
   plot -c red   -pt 7 -sz 2.0 imag_V dmag_V
   plot -c darkgreen -pt 3 -sz 2.0 imag_V Dmag_V
@@ -264,4 +279,5 @@
   sprintf line "OFFSET: %6.3f" $offset
   textline -frac 0.1 0.8 -fn courier 24 "$line"
+
 
   set lChiNorm = log(PSF_CHISQ_out / PSF_NDOF_out)