# GALFIELD - Uniform galaxy field

file	image, dat

image = s1
dat = mktemp ("art")

gallist (dat, 1000, interactive=no, spatial="uniform", xmin=1.,
    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
    core_radius=50., base=0., sseed=i+1, luminosity="powlaw", minmag=-7.,
    maxmag=0., mzero=15., power=0.45, alpha=-1.24, mstar=-21.41, lseed=i+1,
    egalmix=0.4, ar=0.7, eradius=10., sradius=1., absorption=1.2, z=0.05,
    sfile="", nssample=100, sorder=10, lfile="", nlsample=100, lorder=10,
    rbinsize=10., mbinsize=0.5, dbinsize=0.5, ebinsize=0.1, pbinsize=20.,
    graphics="stdgraph", cursor="")

starlist (dat, 100, "", "", interactive=no, spatial="uniform", xmin=1.,
    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
    core_radius=30., base=0., sseed=i, luminosity="powlaw", minmag=-7.,
    maxmag=0., mzero=-4., power=0.6, alpha=0.74, beta=0.04, delta=0.294,
    mstar=1.28, lseed=i, nssample=100, sorder=10, nlsample=100, lorder=10,
    rbinsize=10., mbinsize=0.5, graphics="stdgraph", cursor="")

mkobjects (image, output="", ncols=512, nlines=512,
    title="Example artificial galaxy field", header="artdata$stdheader.dat",
    background=400., objects=dat, xoffset=0., yoffset=0., star="moffat",
    radius=1.0, beta=2.5, ar=1., pa=0., distance=1., exptime=1., magzero=5.5,
    gain=5., rdnoise=10., poisson=yes, seed=j, comments=b1)

delete (dat, verify=no)