-
Notifications
You must be signed in to change notification settings - Fork 8
/
params.txt
306 lines (245 loc) · 16.7 KB
/
params.txt
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
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
%-------------------------------------------------------------------------
%---- This file contains the input parameters needed at run-time for
% simulations. It is based on and closely resembles the GADGET-3
% parameterfile (format of which and parsing routines written by
% Volker Springel [[email protected]]). It has been updated
% with new naming conventions and additional variables as needed by
% Phil Hopkins [[email protected]] for GIZMO.
%-------------------------------------------------------------------------
%---- Relevant files (filenames and directories)
InitCondFile NAME
OutputDir OUTFOLDER
%---- File formats (input and output)
ICFormat 3 % 1=unformatted (gadget) binary, 3=hdf5, 4=cluster
SnapFormat 3 % 1=unformatted (gadget) binary, 3=hdf5
%---- Output parameters
RestartFile restart
SnapshotFileBase snapshot
OutputListOn 0 % =1 to use list in "OutputListFilename"
OutputListFilename output_times.txt % list of times (in code units) for snaps
NumFilesPerSnapshot 1
NumFilesWrittenInParallel 1 % must be < N_processors & power of 2
%---- Output frequency
TimeOfFirstSnapshot 0. % time (code units) of first snapshot
TimeBetSnapshot DTSNAP % time between (if OutputListOn=0), code units
TimeBetStatistics DTSNAP % time between additional statistics (e.g. energy)
%---- CPU run-time and checkpointing time-limits
TimeLimitCPU 172800 % in seconds
CpuTimeBetRestartFile 7200 % in seconds
ResubmitOn 0
ResubmitCommand my-scriptfile
%---- Desired simulation beginning and end times (in code units) for run
TimeBegin 0.0 % Beginning of the simulation
TimeMax TMAX % End of the simulation
%---- Maximum and minimum timesteps allowed
MaxSizeTimestep MAXTIMESTEP % in code units, set for your problem
MinSizeTimestep 1.0e-15 % set this very low, or get the wrong answer
%---- System of units
UnitLength_in_cm UNIT_L % should be 1 pc
UnitMass_in_g UNIT_M % should be 1 solar mass
UnitVelocity_in_cm_per_s UNIT_V % should be 1 m/sec
UnitMagneticField_in_gauss UNIT_B % should be 1 Tesla
GravityConstantInternal 0 % calculated by code if =0
%---- Cosmological parameters
ComovingIntegrationOn 0 % is it cosmological? (yes=1, no=0)
BoxSize BOXSIZE % in code units
Omega0 0. % =0 for non-cosmological
OmegaLambda 0. % =0 for non-cosmological
OmegaBaryon 0. % =0 for non-cosmological
HubbleParam 1. % little 'h'; =1 for non-cosmological runs
%----- Memory allocation
%MaxMemSize 5000 % sets maximum MPI process memory use in MByte
PartAllocFactor 5.0 % memory load allowed for better cpu balance
BufferSize 500 % in MByte
%---- Rebuild domains when >this fraction of particles active
TreeDomainUpdateFrequency 0.005 % 0.0005-0.05, dept on core+particle number
%---- (Optional) Initial hydro temperature & temperature floor (in Kelvin)
InitGasTemp 20 % set by IC file if =0
MinGasTemp 2.73 % don't trust cooling function below ~5K, whether dust, molecular lines, or photoelectric heating
%---- Hydro reconstruction (kernel) parameters
DesNumNgb 32 % domain-reconstruction kernel number: 32 standard, 60-114 for quintic
MaxHsml 1.0e10 % minimum gas kernel length (some very large value to prevent errors)
MinGasHsmlFractional 0 % minimum kernel length relative to gas force softening (<= 1)
%---- Gravitational softening lengths
%----- Softening lengths per particle type. If ADAPTIVE_GRAVSOFT is set, these
%-------- are the minimum softening allowed for each type -------
%-------- (units are co-moving for cosmological integrations)
SofteningGas 1e-10 % gas (particle type=0) (in co-moving code units)
SofteningHalo 0.020 % dark matter/collisionless particles (type=1)
SofteningDisk 0.150 % collisionless particles (type=2)
SofteningBulge 0.500 % collisionless particles (type=3)
SofteningStars SOFTENING % stars spawned from gas (type=4)
SofteningBndry SOFTENING % black holes (if active), or collisionless (type=5)
%---- if these are set in cosmo runs, SofteningX switches from comoving to physical
%------- units when the comoving value exceeds the choice here
%------- (these are ignored, and *only* the above are used, for non-cosmo runs)
SofteningGasMaxPhys 0.0005 % e.g. switch to 0.5pc physical below z=1
SofteningHaloMaxPhys 0.010
SofteningDiskMaxPhys 0.075
SofteningBulgeMaxPhys 0.250
SofteningStarsMaxPhys 0.0005
SofteningBndryMaxPhys 0.0005
%----- parameters for adaptive gravitational softening
AGS_DesNumNgb 32 % neighbor number for calculating adaptive gravsoft
%-------------------------------------------------------------------------
%-------------------------------------------------------------------------
%---------- Physics Modules ----------------------------------------------
%-------------------------------------------------------------------------
%-------------------------------------------------------------------------
%------------------------------------------------------------
%------------------ Additional Fluid Physics ----------------
%------------------------------------------------------------
%---- Magneto-Hydrodynamics Parameters (MAGNETIC on)
%----- Initial B-Field Strengths (if MHD_B_SET_IN_PARAMS on, otherwise read from IC file)
BiniX 1.0e-8 % initial B_x, in code units
BiniY 1.0e-8 % initial B_y, in code units
BiniZ 1.0e-8 % initial B_z, in code units
%---- Thermal Conduction (CONDUCTION on)
%----- set coefficient kappa [code units] or, if CONDUCTION_SPITZER on, multiplies value
ConductionCoeff 1.0 % set/multiply conduction coefficient
%---- Navier-Stokes Viscosity (VISCOSITY on)
%--- set coefficients eta,zeta [code units] or, if VISCOSITY_BRAGINSKII on, multiplies value
ShearViscosityCoeff 1.0 % set/multiply shear viscosity coefficient
BulkViscosityCoeff 1.0 % set/multiply bulk viscosity coefficient
%---- Turbulent Diffusion Master Switch (TURB_DIFFUSION on)
TurbDiffusionCoefficient 1.0 % Normalizes diffusion rates relative to Smagorinsky-Lilly theory [best calibration] (~0.5-2)
%---- Cosmic Ray + Gas Fluids (COSMIC_RAYS on)
CosmicRayDiffusionCoeff 1.0 % multiplies anisotropic diffusion/streaming coefficients
%---- Dust-Gas Mixtures (GRAIN_FLUID on)
Grain_Internal_Density 1.0 % internal/material density of grains in g/cm^3
Grain_Size_Min 1.e-6 % minimum grain size in cm
Grain_Size_Max 1.e-4 % maximum grain size in cm
Grain_Size_Spectrum_Powerlaw 0.5 % power-law distribution of grain sizes (dm/dlnr~r^x)
%-------------------------------------------------------------------------
%------------------ Star, Black Hole, and Galaxy Formation ---------------
%-------------------------------------------------------------------------
%---- Star Formation parameters (GALSF on)
CritPhysDensity RHOMAX % critical physical density for star formation (cm^(-3))
SfEffPerFreeFall 1.0 % SFR/(Mgas/tfreefall) for gas which meets SF criteria
%---- sub-grid (Springel+Hernquist/GADGET/AREPO) "effective equation of state"
%------- star formation+feedback model (GALSF_EFFECTIVE_EQS on)
MaxSfrTimescale 4.0 % code units (SF timescale at 2-phase threshold)
TempSupernova 3.0e8 % in Kelvin (temp of hot gas in 2-phase model)
TempClouds 1000.0 % in Kelvin (temp of cold gas in 2-phase model)
FactorSN 0.1 % SNe coupling frac (frac of egy retained in hot)
FactorEVP 3000.0 % controls Kennicutt normalization
FactorForSofterEQS 1.0 % interpolate between 'stiff' and isothermal EOS
%------- the sub-grid "decoupled winds" model (GALSF_SUBGRID_WINDS on)
WindEfficiency 2.0 % mass-loading (Mdot_wind = SFR * WindEfficiency)
WindEnergyFraction 0.06 % fraction of SNe energy in winds (sets velocity)
WindFreeTravelMaxTime 0.1 % 'free-stream time' in units of t_Hubble(z)
WindFreeTravelDensFac 0.1 % 'free-stream' until density < this * CritPhysDensity
%------- alternative winds (set GALSF_SUBGRID_WIND_SCALING == 1 or 2)
%------- (scaling with local dark matter properties, as Dave/Oppenheimer/Mannucci/Illustris)
VariableWindVelFactor 1.0 % wind velocity relative to estimated halo v_escape
VariableWindSpecMomentum 5000. % wind momentum per unit stellar mass (code velocity units)
%-------------- FIRE (PFH) explicit star formation & feedback model (FIRE on)
%--- initial metallicity of gas & stars in simulation
InitMetallicity ZINIT % initial gas+stellar metallicity (in solar)
InitStellarAge 0.0 % initial mean age (in Gyr; for stars in sim ICs)
%--- local radiation-pressure driven winds (GALSF_FB_FIRE_RT_LOCALRP)
WindMomentumLoading 1.0 % fraction of photon momentum to couple
%--- SneII Heating Model (GALSF_FB_MECHANICAL)
SNeIIEnergyFrac 1.0 % fraction of mechanical energy to couple
%--- HII region photo-heating model (GALSF_FB_FIRE_RT_HIIHEATING)
HIIRegion_fLum_Coupled 1.0 % fraction of ionizing photons allowed to see gas
%--- long-range radiation pressure acceleration (GALSF_FB_FIRE_RT_LONGRANGE)
PhotonMomentum_Coupled_Fraction 1.0 % fraction of L to allow incident
PhotonMomentum_fUV 0.0 % incident SED f(L) in UV (minimum scattering)
PhotonMomentum_fOPT 0.0 % incident SED f(L) in optical/near-IR
%--- gas return/recycling
GasReturnFraction 1.0 % fraction of gas mass returned (relative to ssp)
GasReturnEnergy 1.0 % fraction of returned gas energy+momentum (relative to ssp)
%--- cosmic rays (COSMIC_RAYS)
CosmicRay_SNeFraction 0.1 % fraction of SNe ejecta kinetic energy into cosmic rays (~10%)
InterstellarRadiationFieldStrength ISRF % interstellar radiation field intensity relative to values measured in Solar neighborhood
%-------------- Black Hole accretion & formation (BLACK_HOLES on)
%--- formation/seeding
SeedBlackHoleMass BH_SEED_MASS % initial mass (on-the-fly or single galaxy)
SeedAlphaDiskMass 0.0 % initial mass in the alpha disk (BH_ALPHADISK_ACCRETION)
SeedBlackHoleMinRedshift 2.0 % minimum redshift where new BH particles are seeded (lower-z ceases seeding)
SeedBlackHoleMassSigma 0.5 % lognormal standard deviation (in dex) in initial BH seed masses
%----- (specific options for on-the-fly friends-of-friends based BH seeding: FOF on)
MinFoFMassForNewSeed 10. % minimum mass of FOF group (stars or DM) to get seed, in code units
TimeBetOnTheFlyFoF 1.01 % time (in code units, e.g. scale-factor) between on-the-fly FOF searches
%--- accretion
BlackHoleAccretionFactor 1.0 % multiplier for mdot (relative to model)
BlackHoleEddingtonFactor 1e100 % fraction of eddington to cap (can be >1)
BlackHoleNgbFactor 1.0 % multiplier for kernel neighbors for BH
BlackHoleMaxAccretionRadius BOXSIZE % max radius for BH neighbor search/accretion (code units)
BlackHoleRadiativeEfficiency 5e-7 % radiative efficiency (for accretion and feedback)
%--- feedback
BlackHoleFeedbackFactor 1.0 % generic feedback strength multiplier
BH_FluxMomentumFactor 0.0 % multiply radiation pressure (BH_PHOTONMOMENTUM), set it to zero to avoid launching gas from rad. pressure
BAL_f_accretion 0.7 % fraction of gas swallowed by BH (BH_WIND options)
BAL_v_outflow 100. % velocity (km/s) of BAL outflow (BH_WIND options)
BAL_internal_temperature 1.0e3 % internal temperature (K) of BAL outflow (BH_WIND_SPAWN)
BAL_wind_particle_mass JET_PART_MASS % mass of 'virtual wind particles' in code units (BH_WIND_SPAWN)
BAL_wind_particle_mass_MS MS_WIND_PART_MASS % mass of 'virtual wind particles' in MS stellar winds, if 0 then BAL_wind_particle_mass is used, in code units (SINGLE_STAR_FB_WINDS)
%-------------------------------------------------------------------------
%------------------ Grackle cooling module -----------------
%-------------------------------------------------------------------------
%-------------- Grackle UVB file (COOL_GRACKLE on)
GrackleDataFile CloudyData_UVB=HM2012.h5
%-------------------------------------------------------------------------
%------------------ Driven Turbulence (Large-Eddy boxes) -----------------
%-------------------------------------------------------------------------
%-------------- Legacy Turbulent stirring parameters (TURB_DRIVING on)
%ST_decay TURBDECAY % decay time for driving-mode phase correlations
%ST_energy TURBENERGY % energy of driving-scale modes: sets norm of turb
%ST_DtFreq TURBFREQ % time interval for driving updates (set by hand)
%ST_Kmin TURB_KMIN % minimum driving-k: should be >=2.*M_PI/All.BoxSize
%ST_Kmax TURB_KMAX % maximum driving-k: set to couple times Kmin or more if cascade desired
%ST_SolWeight 1.0 % fractional wt of solenoidal modes (wt*curl + (1-wt)*div)
%ST_AmplFac 1.0 % multiplies turb amplitudes
%ST_SpectForm 3 % driving pwr-spec: 0=Ek~const; 1=sharp-peak at kc; 2=Ek~k^(-5/3); 3=Ek~k^-2
%ST_Seed 42 % random number seed for modes (so you can reproduce it)
%IsoSoundSpeed 200. % initializes gas sound speed in box to this value
%TimeBetTurbSpectrum DTSNAP % time (code units) between evaluations of turb pwrspec
%-------------- Turbulent stirring parameters (TURB_DRIVING on)
TurbDrive_ApproxRMSVturb TURB_SIGMA
TurbDrive_MinWavelength TURB_MINLAMBDA
TurbDrive_MaxWavelength TURB_MAXLAMBDA
TurbDrive_SolenoidalFraction 1
TurbDrive_CoherenceTime TURB_COHERENCE_TIME
TurbDrive_DrivingSpectrum 1
TurbDrive_RandomNumberSeed 42
TurbDrive_TimeBetweenTurbUpdates DTSNAP
TurbDrive_TimeBetTurbSpectrum DTSNAP
%-------------------------------------------------------------------------------------------------
%------------------ Non-Standard Dark Matter, Dark Energy, Gravity, or Expansion -----------------
%-------------------------------------------------------------------------------------------------
%-------------- Parameters for non-standard or time-dependent Gravity/Dark Energy/Expansion (GR_TABULATED_COSMOLOGY on)
DarkEnergyConstantW -1 % time-independent DE parameter w, used only if no table
TabulatedCosmologyFile CosmoTbl % table with cosmological parameters
%-------------------------------------------------------------
%------------------ Solid bodies and Impacts -----------------
%-------------------------------------------------------------
%-------------- Parameters for custom Tillotson equation-of-state (EOS_TILLOTSON on)
%--- In ICs, set "CompositionType": 0=custom,1=granite,2=basalt,3=iron,4=ice,5=olivine/dunite,6=water;
%--- their EOS parameters will be set accordingly. If CompositionType=0, the custom parameters below
%--- are used, matched to the definitions in Table A1 of Reinhardt+Stadel 2017,MNRAS,467,4252 (below is iron)
Tillotson_EOS_params_a 0.5 % a parameter [dimensionless]
Tillotson_EOS_params_b 1.5 % b parameter [dimensionless]
Tillotson_EOS_params_u_0 9.5e10 % u_0 parameter in [erg/g]
Tillotson_EOS_params_rho_0 7.86 % rho_0 parameter in [g/cm^3]
Tillotson_EOS_params_A 1.28e12 % A parameter in [erg/cm^3]
Tillotson_EOS_params_B 1.05e12 % B parameter in [erg/cm^3]
Tillotson_EOS_params_u_s 1.42e10 % u_s parameter in [erg/g]
Tillotson_EOS_params_u_s_prime 8.45e10 % u_s^prime parameter in [erg/g]
Tillotson_EOS_params_alpha 5.0 % alpha parameter [dimensionless]
Tillotson_EOS_params_beta 5.0 % beta parameter [dimensionless]
Tillotson_EOS_params_mu 7.75e11 % elastic shear modulus in [erg/cm^3] (used if EOS_ELASTIC is on)
Tillotson_EOS_params_Y0 8.5e10 % hugoniot elastic limit in [erg/cm^3] (used if EOS_ELASTIC is on)
%--- Developer-Mode Parameters (usually hard-coded, but set manually if DEVELOPER_MODE is on) --------
ErrTolTheta 0.5 % 0.7=standard
ErrTolForceAcc 0.0025 % 0.0025=standard
ErrTolIntAccuracy 0.01 % <0.02
CourantFac 0.2 % <0.20
MaxRMSDisplacementFac 0.125 % <0.25
MaxNumNgbDeviation 0.05 % <<DesNumNgb (values<1 are fine)
AGS_MaxNumNgbDeviation 0.10 % tolerance in Nngb (make larger than gas)
%---- Magneto-Hydrodynamics Developer-Mode Parameters (MAGNETIC on)
%--- Dedner 2002 div-cleaning parameters
DivBcleaningParabolicSigma 1.0 % (3D~0.7-1.0,2D~0.2-0.3)
DivBcleaningHyperbolicSigma 1.0 % (~1)