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
|
# MULTIFIBER - Multifiber example
file out, obj, dat
out = s1
obj = mktemp ("art")
dat = mktemp ("art")
mk1dspec (obj, output="", ap=1, rv=0., z=no, title="", header="",
ncols=512, naps=1, wstart=4000., wend=8000., continuum=1000., slope=0.,
temperature=5700., lines="", nlines=50, peak=-0.5, profile="gaussian",
gfwhm=24, seed=i, comments=b1)
print (obj, " .1 gauss 3 0 10.1 .002", > dat)
print (obj, " .2 gauss 3 0 20.2 .002", >> dat)
print (obj, " .3 gauss 3 0 30.3 .002", >> dat)
print (obj, " .4 gauss 3 0 40.4 .002", >> dat)
print (obj, " .5 gauss 3 0 50.5 .002", >> dat)
print (obj, " .6 gauss 3 0 60.6 .002", >> dat)
print (obj, " .7 gauss 3 0 70.7 .002", >> dat)
print (obj, " .8 gauss 3 0 80.8 .002", >> dat)
print (obj, " .9 gauss 3 0 90.9 .002", >> dat)
mk2dspec (out, output="", model=dat, comments=b1, ncols=100, nlines=512,
title="Example artificial multifiber image", header="artdata$stdheader.dat")
mknoise (out, output="", title="", header="", ncols=512, nlines=512,
background=0., gain=1., rdnoise=3., poisson=yes, seed=j, cosrays="",
ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
imdelete (obj, verify=no)
delete (dat, verify=no)
|
