1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
# 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)
|
