FastSim/INPUT.cfg at master · vrastil/FastSim · GitHub

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
# **************************
# * SIMULATION BOX OPTIONS *
# **************************

mesh_num = 64		# number of mesh cells per dimension (potential)
mesh_num_pwr = 128 # number of mesh cells per dimension (power spectrum)
par_num = 64		# number of particles per dimension
box_size = 1000		# box size in units of Mpc/h

# ***********************
# * INTEGRATION OPTIONS *
# ***********************

redshift = 20		# redshift at the start of the simulation
redshift_0 = 0		# redshift at the end of the simulation
time_step = 0.1		# dimensionless time-step (scale factor)

# ******************
# * OUTPUT OPTIONS *
# ******************

print_every = 5     # save particle positions and power spectrum every n-th step, set 0 for no printing
pwr_bins = 10    # number of bins per decade in power spectrum
corr_pt = 10    # number of points per 10 Mpc in correlation function
out_dir = jobs/output/test_runs/	# output folder name

print_par_pos = 0   # print particles positions
print_dens = 0   # print density map and histogram
print_pwr = 1   # print power spectrum
print_extrap_pwr = 0   # print extrapolated power spectrum
print_corr = 0   # print correlation function
print_vel_pwr = 0   # print velocity power spectrum

# ******************
# * APPROXIMATIONS *
# ******************

comp_ZA = 1		# compute Zeldovich approximation
comp_TZA = 0		# compute Zeldovich approximation
comp_FF = 0		# compute Frozen-flow approximation
comp_FP = 0		# compute Frozen-potential approximation
comp_PM = 1		# compute Particle-mesh approximation
comp_AA = 0 	# compute Adhesion approximation
comp_FP_pp = 0	# compute Frozen-potential approximation (particle-particle interaction)

# ************************
# * CHAMELEON PARAMETERS *
# ************************

comp_chi = 0        # compute chameleon gravity (frozen-potential)
comp_chi_ff = 0        # compute chameleon gravity (frozen-flow)
comp_chi_lin = 0    # use only linear prediction in k-space
# chi_beta = 0,4    # coupling constant
chi_n = 0.5         # chameleon power-law potential exponent,0 < n < 1
chi_phi = 1E-5      # screening potential

# *******************
# * TEST PARAMETERS *
# *******************

R_sphere = 1
rho_sphere = 3E4
N_grid = 64
N_min = 16
fine_sweeps = 15
coarse_sweeps = 5
max_steps = 40
step_per_iter = 10
verbose = 0

# ***************************
# * COSMOLOGICAL PARAMETERS *
# ***************************

Omega_b = 0.0486  # density of baryons relative to the critical density
Omega_m = 0.3089   # density of all relative to the critical density
Hubble = 67.74     # Hubble constant in units of km/s/Mpc
n_s = 0.9667		    # power-law index of the scale-free power spectrum
sigma8 = 0.8159		# normalization of the power spectrum at R = 8 Mpc/h
smoothing_k = 0.5		# smoothing wavenumber of TZA in units of h/Mpc, set 0 for ZA
    # -------------------------
    # TRANSFER FUNCTION METHODS
    # -------------------------
    # ccl_emulator = 0
    # ccl_eisenstein_hu  = 1
    # ccl_bbks = 2
    # ccl_boltzmann_class = 3
    # ccl_boltzmann_camb = 4
    # -------------------------
    transfer_function = 0

    # ----------------------------
    # MATER POWER SPECTRUM METHODS
    # ----------------------------
    # ccl_linear = 0
    # ccl_halofit = 1
    # ccl_halo_model = 3
    # ccl_emu = 4
    # -----------------------------
    matter_power_spectrum = 4

    # -----------------------------
    # BARYON POWER SPECTRUM METHODS
    # -----------------------------
    # ccl_nobaryons = 0
    # ccl_bcm = 1
    # ------------------------------
    baryons_power_spectrum = 0

    # ---------------------
    # MASS FUNCTION METHODS
    # ---------------------
    # ccl_tinker = 1
    # ccl_tinker10 = 2
    # ccl_watson = 3
    # ccl_angulo = 4
    # ccl_shethtormen = 5
    # ---------------------
    mass_function = 2

# ***************
# * RUN OPTIONS *
# ***************

num_thread = 32		# number of threads the program will use
seed = 0 # seed to random number generator, use 0 for random
pair = 0 # if true run two simulations with opposite phases of random field
mlt_runs = 1 # how many runs should be simulated (only if seed = 0)