III. 3. B. Doing the preliminaries: prelim.cl

The source code for prelim.cl can be obtained here.

We execute the script prelim.cl, which runs daofind and phot on an image. We adopt an aperture radius of 3.5, and a sky annulus of 12 pixels to 20 pixels.

# Written by Phil Massey, Lowell Obs.
# Modified to use keywords about keywords March 21, 2002
        procedure prelim (img)
        string img {prompt="Images for daofind and phot"}

        begin

        string s1, s2, s3

	real fwhmv,sigpv,datanv,dataxv

# Get parameters
	s1=img
	imgets(s1,"datami")
	datanv=real(imgets.value)
	imgets(s1,"datama")
	dataxv=real(imgets.value)
	imgets(s1,"fwhm")
	fwhmv=real(imgets.value)
	imgets(s1,"sigma")
	sigpv=real(imgets.value)
	imgets(s1,"fname")
	s2=imgets.value
	imgets(s1,"timestamp")
	s3=imgets.value
	print("Running daofind and phot on:",(s1)," with datamin:",(datanv)," datamax:",(dataxv)," fwhm:",(fwhmv),"sigma=",(sigpv))

	delete(s1//".coo.*")
daofind(s1,"default",thresh=4.0,datamax=dataxv,datamin=datanv,fwhmpsf=fwhmv,sigma=sigpv,interact-,verify-,verbose-,filter=s2,obstime=s3)
	delete(s1//".mag.*")
phot(s1,"default","default",datamax=dataxv,datamin=datanv,apertur=3.5,calgori="none",interact-,verify-,verbose-,filter=s2,obstime=s3,annulu=12,dannu=8)

end