-
Notifications
You must be signed in to change notification settings - Fork 8
/
Copy pathexample_bransdicke.txt
80 lines (70 loc) · 5.58 KB
/
example_bransdicke.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
% =============================== %
% This is the run parameters file %
% =============================== %
% =============================== %
% Modified gravity parameters
% =============================== %
modified_gravity_active 1 % DGP simulation [1] or LCDM [0]
wBD % The Brans-Dicke parameter
% The physical cosmological parameters. The codes ignores Omega etc. below and computes [h]
% from the requirement that Geff = 1 today
Omegah2 % The physical matter density parameter Omega * h^2
Omegarh2 % The physical radiation density parameter Omegar * h^2
Omegavh2 % The physical dark energy density parameter Omegav * h^2
include_screening 1 % Include screening [1] or not [0]
use_lcdm_growth_factors 0 % Use growth-factors computed in LCDM
input_pofk_is_for_lcdm 1 % If the input P(k) is for LCDM we get the correct MG P(k) by *= (D/DLCDM)^2
input_sigma8_is_for_lcdm 1 % If Sigma8 provided below is for LCDM then we rescale P(k)
inverted_initial_condition 0 % If 1 then delta(k) -> -delta(k) in IC
amplitude_fixed_initial_condition 0 % If 1 then enforce amplitude |delta(k)|^2 == P(k) in IC
% ==================
% Simulation outputs
% ==================
OutputDir output % Directory for output (must exist)
FileBase dgp % Base-filename of output files (appropriate additions are appended on at runtime)
OutputRedshiftFile files/output_redshifts.dat % The file containing the redshifts that we want snapshots for
NumFilesWrittenInParallel 1 % Limits the number of files that are written in parallel when outputting
% =========================
% Simulation Specifications
% =========================
UseCOLA 1 % Whether or not to use the COLA method (1 = true, 0 = false).
Buffer 1.5 % The amount of extra memory to reserve for particles moving between tasks during runtime.
Nmesh 64 % This is the size of the FFT grid used to compute the displacement field and gravitational forces.
Nsample 64 % This sets the total number of particles in the simulation, such that Ntot = Nsample^3.
Box 200.0 % The Periodic box size of simulation.
Init_Redshift 9.0 % The redshift to begin timestepping from (redshift = 9 works well for COLA)
Seed 5001 % Seed for IC-generator
SphereMode 0 % If "1" only modes with |k| < k_Nyquist are used to generate initial conditions (i.e. a sphere in k-space),
% otherwise modes with |k_x|,|k_y|,|k_z| < k_Nyquist are used (i.e. a cube in k-space).
WhichSpectrum 1 % "0" - Use transfer function, not power spectrum
% "1" - Use a tabulated power spectrum in the file 'FileWithInputSpectrum'
% otherwise, Eisenstein and Hu (1998) parametrization is used
% Non-Gaussian case requires "0" and that we use the transfer function
WhichTransfer 0 % "0" - Use power spectrum, not transfer function
% "1" - Use a tabulated transfer function in the file 'FileWithInputTransfer'
% otherwise, Eisenstein and Hu (1998) parameterization used
% For Non-Gaussian models this is required (rather than the power spectrum)
FileWithInputSpectrum files/input_power_spectrum.dat % filename of tabulated input spectrum (if used)
% expecting k Pk
FileWithInputTransfer files/input_transfer_function.dat % filename of tabulated transfer function (if used)
% expecting k and T (unnormalized)
% =======================
% Cosmological Parameters
% =======================
Omega 0.267 % Total matter density (CDM + Baryons at z=0).
OmegaBaryon 0.049 % Baryon density (at z=0).
HubbleParam 0.71 % Hubble parameter, 'little h' (may be used for power spectrum parameterization).
Sigma8 0.8 % Power spectrum normalization (power spectrum may already be normalized correctly).
PrimordialIndex 0.966 % Used to tilt the power spectrum for non-tabulated power spectra (if != 1.0 and nongaussian, generic flag required)
% =======================
% Units
% =======================
UnitLength_in_cm 3.085678e24 % defines length unit of output (in cm/h)
UnitMass_in_g 1.989e43 % defines mass unit of output (in g/cm)
UnitVelocity_in_cm_per_s 1e5 % defines velocity unit of output (in cm/sec)
InputSpectrum_UnitLength_in_cm 3.085678e24 % defines length unit of tabulated input spectrum in cm/h.
% Note: This can be chosen different from UnitLength_in_cm
% =============================================================== %
% Additional parameters are needed if we want non-gaussianity
% or lightcone sims. See original PICOLA for this
% =============================================================== %