-
Notifications
You must be signed in to change notification settings - Fork 26
/
Copy pathsample_galaxycatalog_run.cfg
167 lines (143 loc) · 9.56 KB
/
sample_galaxycatalog_run.cfg
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
#suggested configuration file for galaxy catalog hydro run where search star particles
#Configuration file for analysing all particles
#runs 3DFOF algorithm to find initial stellar list, keeps 3dfof info, calculates many properties
#Units currently set to take in as input, Mpc, 1e10 solar masses, km/s, output in same units
#To set temporally unique halo ids, alter Snapshot_value=SNAP to appropriate value. Ie: for snapshot 12, change SNAP to 12
################################
#input options
#set up to use EAGLE HDF input, load gas, star, bh and dark matter
################################
HDF_name_convention=2 #HDF EAGLE naming convention
Input_includes_dm_particle=0 #include dark matter particles in hydro input
Input_includes_gas_particle=0 #include gas particles in hydro input
Input_includes_star_particle=1 #include star particles in hydro input
Input_includes_bh_particle=0 #include bh particles in hydro input
Input_includes_wind_particle=0 #include wind particles in hydro input (used by Illustris and moves particle type 0 to particle type 3 when decoupled from hydro forces). Here shown as example
Input_includes_tracer_particle=0 #include tracer particles in hydro input (used by Illustris). Here shown as example
Input_includes_extradm_particle=0 #include extra dm particles stored in particle type 2 and type 3, useful for zooms
#cosmological run
Cosmological_input=1
################################
#unit options, should always be provided
################################
#units conversion from input input to desired internal unit
Length_input_unit_conversion_to_output_unit=1.0 #default code unit,
Velocity_input_unit_conversion_to_output_unit=1.0 #default velocity unit,
Mass_input_unit_conversion_to_output_unit=1.0 #default mass unit,
#assumes input is in 1e10 msun, Mpc and km/s and output units are the same
#converting hydro quantities
Stellar_age_input_is_cosmological_scalefactor=1
Metallicity_input_unit_conversion_to_output_unit=1.0
Stellar_age_input_unit_conversion_to_output_unit=1.0
Star_formation_rate_input_unit_conversion_to_output_unit=1.0
#set the units of the output by providing conversion to a defined unit
#conversion of output length units to kpc
Length_unit_to_kpc=1000.0
#conversion of output velocity units to km/s
Velocity_to_kms=1.0
#conversion of output mass units to solar masses
Mass_to_solarmass=1.0e10
Metallicity_to_solarmetallicity=1.0
Star_formation_rate_to_solarmassperyear=1.0
Stellar_age_to_yr=1.0
#ensures that output is physical and not comoving distances per little h
Comoving_units=0
#sets the total buffer size in bytes used to store temporary particle information
#of mpi read threads before they are broadcast to the appropriate waiting non-read threads
#if not set, default value is equivalent to 1e6 particles per mpi process, quite large
#but significantly minimises the number of send/receives
#in this example the buffer size is roughly that for a send/receive of 10000 particles
#for 100 mpi processes
MPI_particle_total_buf_size=100000000
################################
#search related options
################################
#how to search a simulation
Particle_search_type=3 #search star particles
#for baryon search
Baryon_searchflag=0 #if 1 search for baryons separately using phase-space search when identifying substructures, 2 allows special treatment in field FOF linking and phase-space substructure search, 0 treat the same as dark matter particles
#for search for substruture
Search_for_substructure=1 #if 0, end search once field objects are found
#also useful for zoom simulations or simulations of individual objects, setting this flag means no field structure search is run
Singlehalo_search=0 #if file is single halo in which one wishes to search for substructure
#additional option for field haloes
Keep_FOF=1 #if field 6DFOF search is done, allows to keep structures found in 3DFOF (can be interpreted as the inter halo stellar mass when only stellar search is used).\n
#minimum size for structures
Minimum_size=100 #min 100 particles
Minimum_halo_size=-1 #if field halos have different minimum sizes, otherwise set to -1.
#for field fof halo search
FoF_Field_search_type=3 #5 3DFOF search for field halos, 4 for 6DFOF clean up of field halos, 3 for 6DFOF with velocity scale distinct for each halo
Halo_3D_linking_length=0.20 #3DFOF linking length in interparticle spacing
Halo_6D_linking_length_factor=0.20 #for 6d fof halo search, decreases the physical linking length to separate galaxies unless strongly interacting
Halo_6D_vel_linking_length_factor=1.0 #for 6d fof halo search, use the dispersion of the initial 3DFOF linked galaxy when running the 6DFOF algorithm
#To search for galaxies defined as well separated phase-space density maxima
FoF_search_type=6 #for galaxies, just run phase-space core search.
Halo_core_search=2 #for galaxies, grow phase-space density cores found using iterative Gaussian mixture model like routine
Use_adaptive_core_search=1 #calculate dispersions in configuration & vel space to determine linking lengths
Use_phase_tensor_core_growth=2 #use full stepped phase-space tensor assignment
Halo_core_ellx_fac=0.7 #how linking lengths are changed when searching for local 6DFOF cores,
Halo_core_ellv_fac=2.0 #how velocity lengths based on dispersions are changed when searching for local 6DFOF cores
Halo_core_ncellfac=0.0001 #fraction of total halo particle number setting min size of a local 6DFOF core
Halo_core_adaptive_sigma_fac=3.0 #used when running fully adaptive core search with phase-space tensors, specifies the width of the physical linking length in configuration space dispersion (think of this as how many sigma to include). Typically values are 2
Halo_core_num_loops=8 #allows the core search to iterate, shrinking the velocity linking length to used till the number of cores identified decreases or this limit is reached. Allows apative search with larger linking length to be robust. Typically values are 5-10
Halo_core_loop_ellx_fac=1.0 #Factor by which physical linking length is changed when running loops for core search. Typically values are 0.5-1.0
Halo_core_loop_ellv_fac=0.8 #Factor by which velocity linking length is changed when running loops for core search. Typically values are 0.5-1.0
Halo_core_loop_elln_fac=1.5 #how much to change the min number of particles per iteration
Halo_core_phase_significance=2.0 #how significant a core must be in terms of dispersions (sigma) significance
################################
#Unbinding options (VELOCIraptor is able to accurately identify tidal debris so particles need not be bound to a structure)
################################
#unbinding related items
Unbind_flag=0 #run unbinding
################################
#Cosmological parameters
#this is typically overwritten by information in the gadget/hdf header if those input file types are read
################################
h_val=1.0
Omega_m=0.3
Omega_Lambda=0.7
Critical_density=1.0
Virial_density=200 #so-called virial overdensity value
Omega_b=0. #no baryons
################################
#Calculation of properties related options
################################
Inclusive_halo_masses=0 #calculate inclusive masses for halos using full Spherical overdensity apertures once all substructures have been found (if substructures are searched for).
#when calculating properties, for field objects calculate inclusive masses
Iterate_cm_flag=1 #do not interatively find the centre-of-mass, giving bulk centre of mass and centre of mass velocity.
Sort_by_binding_energy=1 #sort by binding energy
Reference_frame_for_properties=0 #use the position of the particle with the minimum potential as the point about which properties should be calculated.
#calculate more (sub)halo properties (like angular momentum in spherical overdensity apertures, both inclusive and exclusive)
#aperture related (list must be in increasing order and terminates with , ie: 1,2,3, )
#calculate aperture masses
Calculate_aperture_quantities=1
Number_of_apertures=6
Aperture_values_in_kpc=3,5,10,30,50,100,
Number_of_projected_apertures=3
Projected_aperture_values_in_kpc=10,50,100,
#spherical overdensity related quantities
#number of spherical overdensity thresholds
Number_of_overdensities=0
#calculate radial profiles
Calculate_radial_profiles=1
Number_of_radial_profile_bin_edges=20
#default radial normalisation log rad bins, in proper kpc
Radial_profile_norm=0
Radial_profile_bin_edges=-2.,-1.87379263,-1.74758526,-1.62137789,-1.49517052,-1.36896316,-1.24275579,-1.11654842,-0.99034105,-0.86413368,-0.73792631,-0.61171894,-0.48551157,-0.3593042,-0.23309684,-0.10688947,0.0193179,0.14552527,0.27173264,0.39794001,
################################
#output related
################################
Write_group_array_file=0 #write a group array file
Separate_output_files=0 #separate output into field and substructure files similar to subfind
Binary_output=2 #binary output 1, ascii 0, and HDF 2
#do not output particles residing in the spherical overdensity apertures of halos, only the particles exclusively belonging to halos
Spherical_overdensity_halo_particle_list_output=0
#halo ids are adjusted by this value * 1000000000000 (or 1000000 if code compiled with the LONGINTS option turned off)
#to ensure that halo ids are temporally unique. So if you had 100 snapshots, for snap 100 set this to 100 and 100*1000000000000 will
#be added to the halo id as set for this snapshot, so halo 1 becomes halo 100*1000000000000+1 and halo 1 of snap 0 would just have ID=1
#ALTER THIS as part of a script to get temporally unique ids
Snapshot_value=SNAP
################################
#other options
################################
Verbose=0 #how talkative do you want the code to be, 0 not much, 1 a lot, 2 chatterbox