Index: /trunk/ippToPsps/test/checktest.sh
===================================================================
--- /trunk/ippToPsps/test/checktest.sh	(revision 38959)
+++ /trunk/ippToPsps/test/checktest.sh	(revision 38959)
@@ -0,0 +1,112 @@
+#!/usr/bin/env dvo
+
+macro go
+ echo "--- g batches ---"
+ check.mags.p2 testdata/P2/catdir.cam/B00000001/00000001.FITS g
+ check.mags.p2 testdata/P2/catdir.cam/B00000002/00000002.FITS g
+ check.mags.p2 testdata/P2/catdir.cam/B00000003/00000003.FITS g
+
+ echo "--- r batches ---"
+ check.mags.p2 testdata/P2/catdir.cam/B00000005/00000013.FITS r
+ check.mags.p2 testdata/P2/catdir.cam/B00000006/00000014.FITS r
+ check.mags.p2 testdata/P2/catdir.cam/B00000007/00000015.FITS r
+ check.mags.p2 testdata/P2/catdir.cam/B00000008/00000016.FITS r
+end
+
+macro check.mags.p2
+  if ($0 != 3)
+    echo "USAGE: check.mags (batchfile) (filter)"
+    break
+  end
+
+  local filter
+  $filter = $2
+
+  catdir testdata/catdir.cam
+  skyregion 9 11 19 21
+  avextract RA DEC g:psf:mag g:psf:flux r:psf:flux g:kron:mag g:kron:flux r:kron:flux g:ap:mag g:ap:flux r:ap:flux  
+
+  data $1
+  read -fits Detection ra dec zp expTime airMass psfFlux psfFluxErr apFlux apFluxErr kronFlux kronFluxErr kronRad 
+
+  match2d ra dec RA DEC {0.5/3600}
+
+  reindex psf_bat = psfFlux using index1
+  reindex psf_dvo = $filter\:psf:flux using index2
+
+  reindex kron_bat = kronFlux using index1
+  reindex kron_dvo = $filter\:kron:flux using index2
+
+  reindex ap_bat = apFlux using index1
+  reindex ap_dvo = $filter\:ap:flux using index2
+
+  reindex ra_bat = ra using index1
+  reindex dec_bat = dec using index1
+
+  reindex ra_dvo = RA using index1
+  reindex dec_dvo = DEC using index1
+
+  set dR = 3600.0*(ra_bat - ra_dvo)
+  set dD = 3600.0*(dec_bat - dec_dvo)
+
+  set dF_kron = kron_bat - kron_dvo
+  set dF_psf  = psf_bat  - psf_dvo
+  set dF_ap   = ap_bat   - ap_dvo
+
+  vstat dR
+  vstat dD
+  vstat dF_kron
+  vstat dF_psf
+  vstat dF_ap
+end
+
+
+macro check.mags.st
+  if ($0 != 2)
+    echo "USAGE: check.mags (batchfile)"
+    break
+  end
+
+  local filter
+
+  catdir testdata/catdir.stk
+  skyregion 9 11 19 21
+  avextract RA DEC g:psf:mag r:psf:mag g:kron:mag g:kron:mag r:kron:mag g:ap:mag g:ap:mag r:ap:mag  
+
+  data $1
+  read -fits StackObjectThin gra gdec gPSFMag gApMag gKronMag rPSFMag rApMag rKronMag
+
+  match2d gra gdec RA DEC {0.5/3600}
+
+  reindex ra_bat = gra using index1
+  reindex dec_bat = gdec using index1
+
+  reindex ra_dvo = RA using index1
+  reindex dec_dvo = DEC using index1
+
+  set dR = 3600.0*(ra_bat - ra_dvo)
+  set dD = 3600.0*(dec_bat - dec_dvo)
+
+  vstat dR
+  vstat dD
+
+  foreach filter g r 
+    reindex $filter\_psf_bat = $filter\PSFMag using index1
+    reindex $filter\_psf_dvo = $filter\:psf:mag using index2
+    
+    reindex $filter\_kron_bat = $filter\KronMag using index1
+    reindex $filter\_kron_dvo = $filter\:kron:mag using index2
+    
+    reindex $filter\_ap_bat = $filter\ApMag using index1
+    reindex $filter\_ap_dvo = $filter\:ap:mag using index2
+
+    set dM_$filter\_kron = $filter\_kron_bat - $filter\_kron_dvo
+    set dM_$filter\_psf  = $filter\_psf_bat  - $filter\_psf_dvo
+    set dM_$filter\_ap   = $filter\_ap_bat   - $filter\_ap_dvo
+
+    vstat dM_$filter\_kron
+    vstat dM_$filter\_psf
+    vstat dM_$filter\_ap
+  end
+end
+