-
Notifications
You must be signed in to change notification settings - Fork 5
/
Model_H3.param
executable file
·170 lines (104 loc) · 6.65 KB
/
Model_H3.param
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
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
%------ File and path names, as well as output file format
OutputDir ./Model-H3-256K
OutputFile snap % Base filename of generated sequence of files
SnapFormat 1 % File format selection
%------ Basic structural parameters of model
CC 10.0 % halo concentration
V200 200.0 % circular velocity v_200 (in km/sec)
LAMBDA 0.0 % spin parameter
MD 0.0 % disk mass fraction
MB 0.0 % bulge mass fraction
MBH 0.0 % black hole mass fraction. If zero, no black
% hole is generated, otherwise one at the centre
% is added.
JD 0.00 % disk spin fraction, typically chosen equal to MD
DiskHeight 0.2 % thickness of stellar disk in units of radial scale length
BulgeSize 0.2 % bulge scale length in units of halo scale length
HaloStretch 1.0 % should be one for a spherical halo, smaller than one corresponds to prolate distortion, otherwise oblate
BulgeStretch 1.0 % should be one for a spherical bulge, smaller than one corresponds to prolate distortion, otherwise oblate
%------ Particle numbers in target model
N_HALO 100000 % desired number of particles in dark halo
N_DISK 0 % desired number of collisionless particles in disk
N_BULGE 0 % number of bulge particles
%------ Selection of symmetry constraints of velocity structure
TypeOfHaloVelocityStructure 1 % 0 = spherically symmetric, isotropic
% 1 = spherically symmetric, anisotropic (with beta parameter specified)
% 2 = axisymmetric, f(E, Lz), with specified net rotation
% 3 = axisymmetric, f(E, Lz, I_3), with <vz^2>/<vR^2> specified and net rotation specified
TypeOfDiskVelocityStructure 0 % 0 = spherically symmetric, isotropic
% 1 = spherically symmetric, anisotropic (with beta parameter specified)
% 2 = axisymmetric, f(E, Lz), with specified net rotation
% 3 = axisymmetric, f(E, Lz, I_3), with <vz^2>/<vR^2> specified and net rotation specified
TypeOfBulgeVelocityStructure 0 % 0 = spherically symmetric, isotropic
% 1 = spherically symmetric, anisotropic (with beta parameter specified)
% 2 = axisymmetric, f(E, Lz), with specified net rotation
% 3 = axisymmetric, f(E, Lz, I_3), with <vz^2>/<vR^2> specified and net rotation specified
HaloBetaParameter -1.0 % only relevant for TypeOfHaloVelocityStructure=1
BulgeBetaParameter 0 % only relevant for TypeOfBulgeVelocityStructure=1
HaloDispersionRoverZratio 1.0 % only relevant for TypeOfHaloVelocityStructure=3
DiskDispersionRoverZratio 1.0 % only relevant for TypeOfDiskVelocityStructure=3
BulgeDispersionRoverZratio 1.0 % only relevant for TypeOfBulgeVelocityStructure=3
HaloStreamingVelocityParameter 0.0 % gives the azimuthal streaming velocity in the TypeOf*VelocityStructure=2/3 cases ('k parameter')
DiskStreamingVelocityParameter 1.0 % gives the azimuthal streaming velocity in the TypeOf*VelocityStructure=2/3 cases ('k parameter')
BulgeStreamingVelocityParameter 0.0 % gives the azimuthal streaming velocity in the TypeOf*VelocityStructure=2/3 cases ('k parameter')
%------ Orbit integration accuracy
TorbitFac 10.0 % regulates the integration time of orbits
% (this is of the order of the typical number of orbits per particle)
TimeStepFactorOrbit 0.01
TimeStepFactorCellCross 0.25
%------ Iterative optimization parameters
FractionToOptimizeIndependendly 0.001
IndepenentOptimizationsPerStep 100
StepsBetweenDump 10
MaximumNumberOfSteps 200
MinParticlesPerBinForDispersionMeasurement 100
MinParticlesPerBinForDensityMeasurement 50
%------ Grid dimension and extenstion/resolution
DG_MaxLevel 7
EG_MaxLevel 7
FG_Nbin 256 % number of bins for the acceleration grid in the R- and z-directions
OutermostBinEnclosedMassFraction 0.999 % regulates the fraction of mass of the Hernquist
% halo that must be inside the grid (determines grid extension)
InnermostBinEnclosedMassFraction 0.0000001 % regulates the fraction of mass enclosed by the innermost
% bin (regulates size of innermost grid cells)
MaxVelInUnitsVesc 0.9999 % maximum allowed velocity in units of the local escape velocity
%------ Construction of target density field
SampleDensityFieldForTargetResponse 1 % if set to 1, the code will randomly sample points to construct the density field
SampleParticleCount 100000000 % number of points sampled for target density field
%------ Construction of force field
SampleForceNhalo 0 % number of points to use to for computing force field with a tree
SampleForceNdisk 0
SampleForceNbulge 0
Softening 0.05
%------ Accuracy settings of tree code used in construction of force field
TypeOfOpeningCriterion 1
ErrTolTheta 0.4
ErrTolForceAcc 0.0025
%------ Domain decomposition parameters used in parallel tree code
MultipleDomains 4
TopNodeFactor 4
%------ Parallel I/O paramaters, only affects writing of galaxy files
NumFilesPerSnapshot 1
NumFilesWrittenInParallel 1
%------ Memory allocation parameters
MaxMemSize 2300.0 % in MB
BufferSize 100.0
BufferSizeGravity 100.0
%------ Specification of internal system of units
UnitLength_in_cm 3.085678e21 % 1.0 kpc
UnitMass_in_g 1.989e43 % 1.0e10 solar masses
UnitVelocity_in_cm_per_s 1e5 % 1 km/sec
GravityConstantInternal 0
%------ Factors to tune merit function of time avaraged velocity structure (only ver 1.1)
HaloValueRsFac 1.0
HaloValueTsFac 1.0
HaloValueQsFac 1.0
HaloValuePsFac 1.0
DiskValueRsFac 1.0
DiskValueTsFac 1.0
DiskValueQsFac 1.0
DiskValuePsFac 1.0
BulgeValueRsFac 1.0
BulgeValueTsFac 1.0
BulgeValueQsFac 1.0
BulgeValuePsFac 1.0