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
31
32
33
34
35
36
37
38
39
|
# GALCLUSTER - Galaxy cluster
file image, dat
image = s1
dat = mktemp ("art")
gallist (dat, 100, interactive=no, spatial="hubble", xmin=1., xmax=512.,
ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF, core_radius=50.,
base=0., sseed=i+1, luminosity="schecter", minmag=-7., maxmag=0.,
mzero=15., power=0.6, alpha=-1.24, mstar=-21.41, lseed=i+1, egalmix=0.8,
ar=0.7, eradius=20., 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="")
gallist (dat, 500, 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.6, alpha=-1.24, mstar=-21.41, lseed=i+1,
egalmix=0.4, ar=0.7, eradius=20., 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 cluster", header="artdata$stdheader.dat",
background=1000., objects=dat, xoffset=0., yoffset=0., star="moffat",
radius=1.0, beta=2.5, ar=1., pa=0., distance=1., exptime=1., magzero=7.,
gain=3., rdnoise=10., poisson=yes, seed=j, comments=b1)
delete (dat, verify=no)
|
