Commit 2987a26b authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Added the small cosmological box example provided by Chris Power.

parent 566ab02f
Small LCDM cosmological simulation generated by C. Power. Cosmology
is WMAP9 and the box is 100Mpc/h in size with 64^3 particles.
We use a softening length of 1/25th of the mean inter-particle separation.
The ICs have been generated to run with Gadget-2 so we need to switch
on the options to cancel the h-factors and a-factors at reading time.
MD5 checksum of the ICs:
2a9c603ffb1f6d29f3d98a3ecb9d3238 small_cosmo_volume.hdf5
#!/bin/bash
wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/small_cosmo_volume.hdf5
#!/bin/bash
# Generate the initial conditions if they are not present.
if [ ! -e small_cosmo_volume.hdf5 ]
then
echo "Fetching initial conditions for the small cosmological volume example..."
./getIC.sh
fi
../swift -a -c -G -t 16 small_cosmo_volume.yml 2>&1 | tee output.log
# Define the system of units to use internally.
InternalUnitSystem:
UnitMass_in_cgs: 1.989e43 # 10^10 M_sun in grams
UnitLength_in_cgs: 3.085678e24 # Mpc in centimeters
UnitVelocity_in_cgs: 1e5 # km/s in centimeters per second
UnitCurrent_in_cgs: 1 # Amperes
UnitTemp_in_cgs: 1 # Kelvin
# WMAP9 cosmology
Cosmology:
Omega_m: 0.276
Omega_lambda: 0.724
Omega_b: 0.0455
h: 0.703
a_begin: 0.0196078
a_end: 1.0
# Parameters governing the time integration
TimeIntegration:
dt_min: 1e-6
dt_max: 1e-2
# Parameters for the self-gravity scheme
Gravity:
eta: 0.025
theta: 0.3
comoving_softening: 0.08
max_physical_softening: 0.08
mesh_side_length: 32
# Parameters governing the snapshots
Snapshots:
basename: snap
delta_time: 1.02
scale_factor_first: 0.02
# Parameters governing the conserved quantities statistics
Statistics:
delta_time: 1.02
scale_factor_first: 0.02
Scheduler:
max_top_level_cells: 8
cell_split_size: 50
# Parameters related to the initial conditions
InitialConditions:
file_name: small_cosmo_volume.hdf5
cleanup_h_factors: 1
cleanup_velocity_factors: 1
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment