Index: /branches/eam_branches/ipp-20130711/psphot/test/tap_psphot_galaxygrid.pro
===================================================================
--- /branches/eam_branches/ipp-20130711/psphot/test/tap_psphot_galaxygrid.pro	(revision 36052)
+++ /branches/eam_branches/ipp-20130711/psphot/test/tap_psphot_galaxygrid.pro	(revision 36053)
@@ -392,4 +392,8 @@
 
 macro grid.plots.sersic
+  if ($0 != 2)
+    echo "USAGE: grid.plot.devexp (version)"
+    break
+  end
   # things to examine: theta, Rmajor, AR
 
@@ -403,5 +407,5 @@
   subset dTx = dT if (ARin < 0.95)
   histogram dTx NdT -8 8 0.1 -range dTi
-  lim -n 0 dTi NdT; clear; box; plot -x 1 dTi NdT
+  lim -n 0$1 dTi NdT; clear; box; plot -x 1 dTi NdT
   label -x "angle offset (degrees, only non-circular)" -y "number count"
 
@@ -413,5 +417,5 @@
   create n 0 dR[]
   set dRf = dR / Rin_s
-  lim -n 1 n dRf; clear; box; plot n dRf
+  lim -n 1$1 n dRf; clear; box; plot n dRf
   label -x sequence -y "1 - R_out| / R_in|"
 
@@ -420,6 +424,16 @@
   # lim ARin ARot; clear; box; plot ARin ARot
   set fAR = ARot / ARin
-  lim -n 2 n fAR; clear; box; plot n fAR
+  lim -n 2$1 n fAR; clear; box; plot n fAR
   label -x sequence -y "AR_out| / AR_in|"
+
+  # check on magnitude
+  set dM = Mot_s - Min_s
+  lim -n 3$1 n dM; clear; box; plot n dM    
+  label -x sequence -y "M_out| - M_in|"
+
+  # check on index
+  set dI = Iot_s - Iin_s
+  lim -n 4$1 n dI; clear; box; plot n dI    
+  label -x sequence -y "I_out| - I_in|"
 end
 
@@ -486,6 +500,9 @@
 
   # ppImage / psphot on the output
+  break -auto off
   echo ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname
-  exec ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname
+  exec echo ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname >& $outname.log
+  exec ppImage $ppImageConfig $psphotConfig -file $basename.fits $outname >>& $outname.log
+  break -auto on
 end