# Written by Phil Massey, Lowell Obs.
# Modified to use keywords about keywords March 21, 2002
# Brought sky in on 7/20/05 and re-ran on everything
        procedure prelim3 (img)
        string img {prompt="Images for daofind and phot"}

        begin

        string s1, s2, s3

	real fwhmv,sigpv,datanv,dataxv,srd,srdlo,srdhi

# 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,"sround")
	srd=real(imgets.value)
	srdlo=srd-1.
	srdhi=srd+1.
	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),"sround=",(srdlo)," to ",(srdhi))
	delete(s1//".coo.*")
daofind(s1,"default",thresh=10.0,datamax=dataxv,datamin=datanv,fwhmpsf=fwhmv,sigma=sigpv,interact-,verify-,verbose-,filter=s2,obstime=s3,roundlo=srdlo,roundhi=srdhi)
flpr
	delete(s1//".mag.*")
phot(s1,"default","default",datamax=dataxv,datamin=datanv,apertur="3,5,10",calgori="none",interact-,verify-,verbose-,filter=s2,obstime=s3,annulu=10,dannu=8)

flpr
end