Make sky images of sheared galaxy catalog for all pass bands.
To do this we use a publicly available skymaker code.
Note: If you compile skymaker with multi-threads mode, you may have
a bug on the random seed of star position: with exactly same parameters set
including non zero random seed, the star positions differ run by run.
This problem is fatal since we are going to stack images.
If you have this problem, you can quickly avoid this by compile skymaker
with single thread mode, i.e. configure --disable-threads. Also be sure
to have compiled single precision FFTW library with single thread mode.
-- config files:
sky_g.conf,
sky_r.conf,
sky_i.conf,
sky_z.conf,
sky_y.conf,
Below we show the example of cshell script to make many images,
#! /bin/csh
set ns = (3 3 15 15 15)
set filt = ("g" "r" "i" "z" "y")
@ i = 1
while ( $i <= 5 )
@ g = 1
while ( $g <= $ns[$i] )
set n = `printf "%03d" ${g}`
set f = skyimages/sky_${filt}[$i]${n}
sky stuff_${filt}[$i].list4 -c ${f}.conf -IMAGE_NAME ${f}.fits
@ g ++
end
@ i ++
end
The exposure time for each shot is (g,r,i,z,y)=(300,300,60,60,60) secnds, and we try to
make stacking images equivalent to 900 seconds exposure: stack by average, instead of sum.
$ stackn sky_i001.fits sky_i002.fits ...... sky_i015.fits sky_i.fits
-- program:
stackn.f90 /
binary(Linux i386)
-- images :
sky_g.fits,
sky_r.fits,
sky_i.fits,
sky_z.fits,
sky_y.fits,
riz image
-- object list :
sky_g.list,
sky_r.list,
sky_i.list,
sky_z.list,
sky_y.list