mirror of
https://github.com/cosmo-sims/monofonIC.git
synced 2024-09-19 17:03:45 +02:00
updated example conf file
This commit is contained in:
parent
1dce202a1a
commit
143aaaf280
3 changed files with 72 additions and 86 deletions
90
example.conf
90
example.conf
|
@ -1,54 +1,106 @@
|
|||
#########################################################################################
|
||||
# Example conf file for MUSIC2 - monofonIC single resolution simulation ICs
|
||||
#########################################################################################
|
||||
|
||||
#########################################################################################
|
||||
[setup]
|
||||
# number of grid cells per linear dimension for calculations = particles for sc initial load
|
||||
GridRes = 65536
|
||||
GridRes = 128
|
||||
|
||||
# length of the box in Mpc/h
|
||||
BoxLength = 4000
|
||||
BoxLength = 300
|
||||
|
||||
# starting redshift
|
||||
zstart = 49.0
|
||||
zstart = 24.0
|
||||
|
||||
# order of the LPT to be used (1,2 or 3)
|
||||
LPTorder = 1
|
||||
LPTorder = 3
|
||||
|
||||
# also do baryon ICs?
|
||||
DoBaryons = yes
|
||||
DoBaryons = no
|
||||
|
||||
# do mode fixing à la Angulo&Pontzen
|
||||
DoFixing = yes
|
||||
# particle load, can be 'sc' (1x), 'bcc' (2x) or 'fcc' (4x) (increases number of particles by factor!)
|
||||
|
||||
# invert phases (for paired simulations)
|
||||
DoInversion = no
|
||||
|
||||
# particle load, can be 'sc' (1x), 'bcc' (2x) or 'fcc' (4x) (increases number of particles by factor!) or 'glass'
|
||||
ParticleLoad = sc
|
||||
|
||||
# Add a possible constraint field here:
|
||||
#ConstraintFieldFile = initial_conditions.h5
|
||||
#ConstraintFieldName = ic_white_noise
|
||||
# if `ParticleLoad = glass' then specify here where to load the glass distribution from
|
||||
#GlassFileName = glass128.hdf5
|
||||
#GlassTiles = 1
|
||||
|
||||
|
||||
#########################################################################################
|
||||
[cosmology]
|
||||
transfer = CLASS
|
||||
ztarget = 2.5
|
||||
# transfer = eisenstein
|
||||
# transfer = file_CAMB
|
||||
## transfer = ... specifies the Einstein-Boltzmann plugin module
|
||||
# transfer = eisenstein # Eisenstein&Hu fitting formula
|
||||
|
||||
# transfer = file_CAMB # CAMB file to be specified as 'transfer_file = ...'
|
||||
# transfer_file = wmap5_transfer_out_z0.dat
|
||||
|
||||
transfer = CLASS # CLASS module (if enabled in CMake file)
|
||||
|
||||
ztarget = 2.5 # target redshift for CLASS module, output at ztarget will be back-scaled to zstart
|
||||
|
||||
Omega_m = 0.302
|
||||
Omega_b = 0.045
|
||||
Omega_L = 0.698
|
||||
H0 = 70.3
|
||||
sigma_8 = 0.811
|
||||
nspec = 0.961
|
||||
sigma_8 = 0.811
|
||||
# A_s = 2.148752e-09 # can use instead of sigma_8
|
||||
|
||||
ZeroRadiation = true # For Back-scaling: set to false if your simulation code can deal with Omega_r!=0
|
||||
|
||||
# Additional cosmological parameters (set by default to the given values)
|
||||
# w0 = -1.0
|
||||
# wa = 0.0
|
||||
# Tcmb = 2.7255
|
||||
# Neff = 3.046
|
||||
|
||||
# anisotropic large scale tidal field
|
||||
# see Stuecker+2020
|
||||
# LSS_aniso_lx = +0.1
|
||||
# LSS_aniso_ly = +0.1
|
||||
# LSS_aniso_lz = -0.2
|
||||
|
||||
#########################################################################################
|
||||
[random]
|
||||
## generator = ... specifies the random field generator plugin module
|
||||
generator = NGENIC
|
||||
seed = 9001
|
||||
|
||||
# generator = PANPHASIA
|
||||
# descriptor = [Panph1,L10,(800,224,576),S9,CH1564365824,MXXL]
|
||||
|
||||
# Add a possible constraint field here:
|
||||
# ConstraintFieldFile = initial_conditions.hdf5
|
||||
# ConstraintFieldName = ic_white_noise
|
||||
|
||||
|
||||
#########################################################################################
|
||||
[testing]
|
||||
# enables diagnostic output
|
||||
# can be 'none' (default), 'potentials_and_densities', 'velocity_displacement_symmetries', or 'convergence'
|
||||
test = none
|
||||
|
||||
|
||||
#########################################################################################
|
||||
[execution]
|
||||
NumThreads = 8
|
||||
|
||||
|
||||
#########################################################################################
|
||||
[output]
|
||||
## format = .... specifies the output plugin module
|
||||
|
||||
# format = grafic2
|
||||
# filename = ics_ramses
|
||||
# grafic_use_SPT = no # if no then uses PPT, otherwise linear SPT
|
||||
|
||||
# format = gadget2
|
||||
# filename = ics_gadget.dat
|
||||
# UseLongids = false
|
||||
|
@ -56,14 +108,10 @@ NumThreads = 8
|
|||
#format = gadget_hdf5
|
||||
#filename = ics_gadget.hdf5
|
||||
|
||||
format = AREPO
|
||||
filename = ics_arepo.hdf5
|
||||
# format = AREPO
|
||||
# filename = ics_arepo.hdf5
|
||||
|
||||
# format = generic
|
||||
# filename = debug.hdf5
|
||||
# generic_out_eulerian = yes
|
||||
|
||||
# format = grafic2
|
||||
# filename = ics_ramses
|
||||
# grafic_use_SPT = yes
|
||||
# generic_out_eulerian = yes # if yes then uses PPT for output
|
||||
|
||||
|
|
62
ics.conf
62
ics.conf
|
@ -1,62 +0,0 @@
|
|||
[setup]
|
||||
# number of grid cells per linear dimension for calculations = particles for sc initial load
|
||||
GridRes = 128
|
||||
# length of the box in Mpc/h
|
||||
BoxLength = 200
|
||||
# starting redshift
|
||||
zstart = 24.0
|
||||
# order of the LPT to be used (1,2 or 3)
|
||||
LPTorder = 1
|
||||
# also do baryon ICs?
|
||||
DoBaryons = no
|
||||
# do mode fixing à la Angulo&Pontzen
|
||||
DoFixing = yes
|
||||
# particle load, can be 'sc' (1x), 'bcc' (2x), 'fcc' (4x), or 'rsc' (8x)
|
||||
ParticleLoad = sc
|
||||
|
||||
[testing]
|
||||
# enables diagnostic output
|
||||
# can be 'none' (default), 'potentials_and_densities', 'velocity_displacement_symmetries', or 'convergence'
|
||||
#test = potentials_and_densities
|
||||
#test = convergence
|
||||
test = none
|
||||
|
||||
[execution]
|
||||
NumThreads = 1
|
||||
|
||||
[output]
|
||||
fname_hdf5 = output.hdf5
|
||||
fbase_analysis = output
|
||||
|
||||
#format = gadget2
|
||||
#filename = ics_gadget.dat
|
||||
|
||||
format = generic
|
||||
filename = debug.hdf5
|
||||
#generic_out_eulerian = yes
|
||||
|
||||
#format = grafic2
|
||||
#filename = ics_ramses
|
||||
#grafic_use_SPT = yes
|
||||
|
||||
[random]
|
||||
generator = NGENIC
|
||||
seed = 9001
|
||||
|
||||
[cosmology]
|
||||
transfer = eisenstein
|
||||
#transfer = CLASS
|
||||
#transfer = eisenstein_wdm
|
||||
#WDMmass = 0.1
|
||||
Omega_m = 0.302
|
||||
Omega_b = 0.045
|
||||
Omega_L = 0.698
|
||||
H0 = 70.3
|
||||
sigma_8 = 0.811
|
||||
nspec = 0.961
|
||||
|
||||
# anisotropic large scale tidal field
|
||||
#LSS_aniso_lx = 0.1
|
||||
#LSS_aniso_ly = 0.1
|
||||
#LSS_aniso_lz = -0.2
|
||||
|
|
@ -88,7 +88,7 @@ int Run( config_file& the_config )
|
|||
|
||||
//--------------------------------------------------------------------------------------------------------
|
||||
//! do constrained ICs?
|
||||
const bool bAddConstrainedModes = the_config.contains_key("setup", "ConstraintFieldFile" );
|
||||
const bool bAddConstrainedModes = the_config.contains_key("random", "ConstraintFieldFile" );
|
||||
|
||||
//--------------------------------------------------------------------------------------------------------
|
||||
//! add beyond box tidal field modes following Schmidt et al. (2018) [https://arxiv.org/abs/1803.03274]
|
||||
|
@ -195,8 +195,8 @@ int Run( config_file& the_config )
|
|||
//--------------------------------------------------------------------
|
||||
if( bAddConstrainedModes ){
|
||||
Grid_FFT<real_t,false> cwnoise({8,8,8}, {boxlen,boxlen,boxlen});
|
||||
cwnoise.Read_from_HDF5( the_config.get_value<std::string>("setup", "ConstraintFieldFile"),
|
||||
the_config.get_value<std::string>("setup", "ConstraintFieldName") );
|
||||
cwnoise.Read_from_HDF5( the_config.get_value<std::string>("random", "ConstraintFieldFile"),
|
||||
the_config.get_value<std::string>("random", "ConstraintFieldName") );
|
||||
cwnoise.FourierTransformForward();
|
||||
|
||||
size_t ngrid_c = cwnoise.size(0), ngrid_c_2 = ngrid_c/2;
|
||||
|
|
Loading…
Reference in a new issue