# FINDSTARS -- Find stars based on the first radial moments.

procedure t_findstars ()

pointer	input			# Input catalog
pointer	output			# Output catalog
int	nstars			# Minimum number of stars to find
real	sigmamin		# Minimum <ir1>/sigma
int	nstart			# Maximum initial number of objects
real	lsigma			# Low sigma threshold
real	hsigma			# High sigma threshold
pointer	filter			# Filter options

pointer	sp
int	clgeti(), stflt()
real	clgetr()


begin
	call smark (sp)
	call salloc (input, SZ_FNAME, TY_CHAR)
	call salloc (output, SZ_FNAME, TY_CHAR)
	call salloc (filter, SZ_LINE, TY_CHAR)

	call clgstr ("input", Memc[input], SZ_FNAME)
	call clgstr ("output", Memc[output], SZ_FNAME)
	nstars = clgeti ("nstars")
	sigmamin = clgetr ("sigmamin")
	nstart = clgeti ("nstart")
	lsigma = clgetr ("lsigma")
	hsigma = clgetr ("hsigma")
	call clgstr ("filter", Memc[filter], SZ_LINE)

	call strpak (Memc[input], Memc[input], SZ_FNAME)
	call strpak (Memc[output], Memc[output], SZ_FNAME)
	call strpak (Memc[filter], Memc[filter], SZ_LINE)

	if (stflt (Memc[filter]) == 0)
	    call findstars (Memc[input], Memc[output], nstars, sigmamin,
		nstart, lsigma, hsigma)

	call sfree (sp)
end