#!/bin/csh -f
#
#
# Use Figaro to correct ESO interquadrant row cross talk
# from SOFI
#

if ($#argv < 3) then
   echo "Usage : esoxtalk input factor output"
   echo "        input, output are expected to be Figaro files"
   echo "        factor is a correction factor"
   echo "            1.4e-5 is recommended by ESO"
   echo "            2.4e-5 has been found by me to be useful"
   exit
endif


set file = $argv[1]:r
set factor = $argv[2]
set out = $argv[3]:r

ystract $file 1 1024 /tmp/cut
icmult /tmp/cut $factor /tmp/cut
isubset /tmp/cut xs=1 xe=512 out=/tmp/cut1
isubset /tmp/cut xs=513 xe=1024 out=/tmp/cut2
iadd /tmp/cut1 /tmp/cut2 /tmp/cut
icdiv /tmp/cut 2.0 /tmp/cut
isuper /tmp/cut xs=1024 ys=1 xp=1 out=/tmp/correct1
isuper /tmp/cut xs=1024 ys=1 xp=513 out=/tmp/correct2
iadd /tmp/correct1 /tmp/correct2 /tmp/correct
isysub $file /tmp/correct $out
exit