Index: /branches/eam_branches/ipp-20110213/psphot/test/tap_psphot_stackphot.pro
===================================================================
--- /branches/eam_branches/ipp-20110213/psphot/test/tap_psphot_stackphot.pro	(revision 30720)
+++ /branches/eam_branches/ipp-20110213/psphot/test/tap_psphot_stackphot.pro	(revision 30721)
@@ -763,4 +763,27 @@
 end
 
+macro ckradialflux
+ data $1
+ read -fits Chip.xrad X_APER Y_APER PSF_FWHM APER_FLUX APER_FLUX_ERR APER_FLUX_STDEV APER_FILL
+ read -fits Chip.psf X_PSF Y_PSF PSF_INST_MAG
+ set mag = PSF_INST_MAG
+ # XXX include nradii in header
+ for i 0 11
+   set mA$i = -2.5*log(APER_FLUX:$i)
+   set dM$i = mA$i - PSF_INST_MAG
+   vstat dM$i
+   $DM$i = $MEDIAN
+ end
+ for i 1 11
+   $j = $i - 1
+   echo {$DM$i - $DM$j}
+ end
+
+ lim mag -2 2; clear; box
+ for i 0 11
+   plot mag dM$i -c black
+ end
+end
+
 if ($SCRIPT)
   fulltest 4