#!/bin/tcsh
#

foreach colour (i)# r)

# Make an lc catalogue from the ascii shear catalogue

lc -C -n number -n flux -n flux_err -n mag -n mag_err -n isoarea -N '1 2 x' -n a -n b -n fr -n fwhm -n flags -n class -n rg -N '1 2 eobs' -N '1 2 ecor' -n Pgamma < {$CATDIR}/D1_{$colour}.shear.cat > D1_{$colour}.shear_lc.cat
   
# Make an "Emode" catalogue including S/N cuts and create a Emode mass map

lc +all 'e1 = %ecor[0] %Pgamma /' 'e2 = %ecor[1] %Pgamma /'< D1_{$colour}.shear_lc.cat | lc +all 'e = %e1 %e2 2 vector' | lc -i '%e %e dot 1 < %flux %flux_err / 10 > and ' > D1_{$colour}.shear_E.cat
makeshearimage -N 200 200 -X 25.0 < D1_{$colour}.shear_E.cat > D1_{$colour}.shear_E.fits
massmap_ft < D1_{$colour}.shear_E.fits | smooth -g 7 7 0 > D1_{$colour}.Emode.fits

# Make an "Bmode" catalogue including S/N cuts and create a Bmode mass map

lc +all 'e1 = %ecor[1] %Pgamma /' 'e2 = %ecor[0] %Pgamma / -1 mult'< D1_{$colour}.shear_lc.cat | lc +all 'e = %e1 %e2 2 vector' | lc -i '%e %e dot 1 < %flux %flux_err / 10 > and ' > D1_{$colour}.shear_B.cat
makeshearimage -N 200 200 -X 25.0 < D1_{$colour}.shear_B.cat > D1_{$colour}.shear_B.fits
massmap_ft < D1_{$colour}.shear_B.fits | smooth -g 7 7 0 > D1_{$colour}.Bmode.fits