KSBf90: A Fortran 90 Implementation of Kaiser, Squires and Broadhurst 1995.

Practical 1:


In this Practical you will use KSBf90 to analyse your SkyMaker simulated image.


  1. 1)Download and unpack KSBf90.tar


This contains a Makefile and pre-compiled 32-bit Linux executables.  At the IAP only some of the computers have the cfitsio routines installed so the Makefile will not compile on all machines.  However you should be able to use the pre-compiled executables.



  1. 2)Set your directories as follows (edit)


setenv dataDIR /Users/heymans/Paris_DUEL/images        # Where your images are

setenv CATDIR /Users/heymans/Paris_DUEL/catalogues     # Where your catalogues are

setenv KSBDIR /Users/heymans/Paris_DUEL/KSBf90         # Where your KSBf90 software is




  1. 3)Run SExtractor on your SkyMaker image. 


Put the resulting Sextractor ascii catalogue in your $CATDIR



4) Download and unpack scripts.tar.


This directory has two example cshell scripts that you will need to edit to include the names of the data your wish to analyse.  It also contains your KSB parameter file KSBf90.param.  You also need to edit this file to update it with the SExtractor catalogue format (see  KSBf90.param website).




5) Run runstars.scr

This script runs findstars.f90 and psffit.f90.




6) Run rungals.scr

This script runs gal_correct.f90 and gal_select.f90.




  1. 7)Did you measure the shear that you put into the simulated image? If you have time, share your simulated data with other students.  Don’t tell the other students what shear you put into your image so you can analyse the shear blindly.  With many simulated images you can make a plot of (g - g_true) vs g_true


Page maintained by Catherine Heymans, IfA, University of Edinburgh