-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathEXAMPLES
127 lines (90 loc) · 4.26 KB
/
EXAMPLES
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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
How to perform specific common tasks using Starlab tools.
--------------------------------------------------------
Note: Most tools are just simple interfaces onto the corresponding
library functions, so the tasks listed below could in principle
also be carried out by compiled programs rather than by
pipes. However, the use of pipes is in many ways clearer and
much more flexible.
For more information on Starlab tools, see the file TOOLS in
this directory. For details on a specific tool, type
tool-name --help
* Create a linked list of 100 equal-mass nodes of unit total mass
makenode -n 100 -m 1
* Create a system of 100 nodes with a Salpeter mass spectrum with masses
in the range 0.5 to 10
makenode -n 100 | makemass -f 1 -x -2.35 -l 0.5 -u 10
* Create a system of 100 nodes with a mass spectrum and evolve the
stars without dynamics
makenode -n 100 | makemass -f 1 -x -2.35 -l 0.5 -u 10 | ???Simon???
* Create a 500-particle Plummer model, with numbered stars, scaled to
standard dynamical units
makeplummer -n 500 -i
* Create a 500-particle W0 = 5 King model, with numbered stars,
unscaled
makeking -n 500 -w 5 -i -u
* Create a 500-particle W0 = 5 King model with a Miller-Scalo mass
spectrum between 0.1 and 20 solar masses, then rescale to unit total
mass, total energy -0.25, and virial ratio 0.5 and display the
results graphically
makeking -n 500 -w 5 -i -u \
| makemass -f 2 -l 0.1 -u 20 \
| scale -m 1 -e -0.25 -q 0.5 \
| xstarplot -l 5 -P .5
* Create a 500-particle W0 = 5 King model with a Miller-Scalo mass
spectrum between 0.1 and 20 solar masses, add in a 10 percent 1-10 kT
binary population, then rescale to unit total mass, total energy
(top-level nodes) -0.25, and virial ratio (top-level nodes) 0.5, and
finally verify the results by analyzing the final snapshot
makeking -n 500 -w 5 -i -u \
| makemass -f 2 -l 0.1 -u 20 \
| makesecondary -f 0.1 -l 0.25 \
| makebinary -l 1 -u 10 \
| scale -m 1 -e -0.25 -q 0.5 \
| sys_stats -n -s
* Evolve this model without stellar evolution for 100 dynamical times,
with log output ever dynamical time and snapshot output every 10
dynamical times, with a self-consistent tidal field, removing
escapers when they are more than two Jacobi radii from the cluster
center
makeking -n 500 -w 5 -i -u \
| makemass -f 2 -l 0.1 -u 20 \
| makesecondary -f 0.1 -l 0.25 \
| makebinary -l 1 -u 10 \
| scale -m 1 -e -0.25 -q 0.5 \
| dstar_kira -t 100 -d 1 -D 10 -Q -G 2
* Create a King model with a power-law mass spectrum and a binary
population, then evolve it with stellar and binary evolution
makeking -n 500 -w 5 -i -u \
| makemass -f 1 -x -2.0 -l 0.1 -u 20 \
| makesecondary -f 0.1 -l 0.1 \
| add_star -Q 0.5 -R 5 \
| scale -M 1 -E -0.25 -Q 0.5 \
| makebinary -f 1 -l 1 -u 1000 -o 2 \
| dstar_kira -t 100 -d 1 -D 10 -f 0.3 -n 10 -q 0.5 -Q -G 2 -S -B
* Perform a series of 100 3-body scattering experiments involving an
equal-mass circular binary and a double-mass incomer, with impact
parameter equal to the binary semimajor axis, relative velocity at
infinity half that needed for zero total energy, and all other
parameters chosen randomly, and display the results as a movie
scatter3 -m 0.5 -e 0 -M 1 -r 1 -v 0.5 -n 100 -C 5 -D 0.1 \
| xstarplot -l 4
* Compute cross-sections for interactions between a circular binary
with component masses 0.75 and 0.25 and an incoming star of mass 1
and velocity at infinity 0.1, all stars having radius 0.05 binary
semimajor axes
sigma3 -d 100 -m 0.25 -e 0 -M 1 -v 0.1 -x 0.05 -y 0.05 -z 0.05
* Create a scattering configuration involving a head-on collision
between a circular binary and a stable hierarchical triple, and
verify the result
mkscat -M 1.5 -r 0 -v 1 -t -a 1 -e 0 -p -a 1 -e 0 -p1 -a 0.1 -e 0 \
| flatten | make_tree -D 1 | pretty_print_tree
* Create a scattering configuration involving a head-on collision
between a circular binary and a stable hierarchical triple, and
integrate it forward in time
scatter -i "-M 1.5 -r 0 -v 1 -t -a 1 -e 0 -p -a 1 -e 0 -p1 -a 0.1 -e 0" \
-t 100 -d 1 -v
* Run kira by calling
kira < [infile] -t 1
* Run kira using the parallel MPI-enabeled version
mpirun -np 16 kira -R [infile] -t 1
(etc.)