diff --git a/examples/Makefile.am b/examples/Makefile.am index a36be246e74496757b0a02db01753e5c7140dafd..ff32a0bdb6d450c6e5e1418342ceaec68c5d4f4d 100644 --- a/examples/Makefile.am +++ b/examples/Makefile.am @@ -84,7 +84,8 @@ EXTRA_DIST = CoolingBox/coolingBox.yml CoolingBox/energy_plot.py CoolingBox/make PerturbedBox_3D/makeIC.py PerturbedBox_3D/perturbedBox.yml PerturbedBox_3D/run.sh \ PMillennium-384/p-mill-384.yml \ PMillennium-768/p-mill-768.yml \ - SantaBarbara/README SantaBarbara/getIC.sh SantaBarbara/santa_barbara.yml \ + SantaBarbara/README SantaBarbara/getIC.sh SantaBarbara/santa_barbara.yml SantaBarbara/run.sh \ + SantaBarbara_low/README SantaBarbara_low/getIC.sh SantaBarbara_low/santa_barbara.yml SantaBarbara_low/run.sh \ SedovBlast_1D/makeIC.py SedovBlast_1D/plotSolution.py SedovBlast_1D/run.sh SedovBlast_1D/sedov.yml \ SedovBlast_2D/getGlass.sh SedovBlast_2D/makeIC.py SedovBlast_2D/plotSolution.py SedovBlast_2D/run.sh SedovBlast_2D/sedov.yml \ SedovBlast_3D/getGlass.sh SedovBlast_3D/makeIC.py SedovBlast_3D/plotSolution.py SedovBlast_3D/run.sh SedovBlast_3D/sedov.yml \ diff --git a/examples/SantaBarbara_low/README b/examples/SantaBarbara_low/README new file mode 100644 index 0000000000000000000000000000000000000000..f86f1a4a4e1d16c3f4011c9e3ed8f35f643bd47e --- /dev/null +++ b/examples/SantaBarbara_low/README @@ -0,0 +1,21 @@ +Initital conditions for the Santa-Barbara cluster comparison project. +These have been regenerated from the orinigal Frenk et al. 1999 paper. + +The cosmology is Omega_m = 1, Omega_b = 0.1, h = 0.5 and sigma_8 = 0.9. + +The ICs are 128^3 particles in a 64^3 Mpc^3 volume. This is about 10x +higher resolution than in the original paper. The ICs have been +created for Gadget and the positions and box size are hence expressed +in h-full units (e.g. box size of 32 / h Mpc). Similarly, the peculiar +velocitites contain an extra sqrt(a) factor. + +We will use SWIFT to cancel the h- and a-factors from the ICs. Gas +particles will be generated at startup. + +MD5 check-sum of the ICS: +1a1600b41002789b6057b1fa6333f3f0 SantaBarbara_128.hdf5 + +You can use the script run_velociraptor.sh to also run a basic 3D FoF +with VELOCIraptor on your output data. You will need to set the +VELOCIRAPTOR_PATH environment variable to tell us where the stf-gas +binary lives. diff --git a/examples/SantaBarbara_low/getIC.sh b/examples/SantaBarbara_low/getIC.sh new file mode 100755 index 0000000000000000000000000000000000000000..759cef50dcfc346b389b1400054fe38358793fdd --- /dev/null +++ b/examples/SantaBarbara_low/getIC.sh @@ -0,0 +1,2 @@ +#!/bin/bash +wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/SantaBarbara_128.hdf5 diff --git a/examples/SantaBarbara_low/run.sh b/examples/SantaBarbara_low/run.sh new file mode 100755 index 0000000000000000000000000000000000000000..044f9f4ef340ab95ffe643104294552469a676ce --- /dev/null +++ b/examples/SantaBarbara_low/run.sh @@ -0,0 +1,4 @@ +#!/bin/bash +# Run SWIFT +../swift -c -s -G -t 28 santa_barbara.yml + diff --git a/examples/SantaBarbara_low/santa_barbara.yml b/examples/SantaBarbara_low/santa_barbara.yml new file mode 100644 index 0000000000000000000000000000000000000000..0e3c66b1a1c1e04bf6fad30e806327b83f03737e --- /dev/null +++ b/examples/SantaBarbara_low/santa_barbara.yml @@ -0,0 +1,60 @@ +# Define the system of units to use internally. +InternalUnitSystem: + UnitMass_in_cgs: 1.98848e43 # 10^10 M_sun in grams + UnitLength_in_cgs: 3.08567758e24 # Mpc in centimeters + UnitVelocity_in_cgs: 1e5 # 1 km/s + UnitCurrent_in_cgs: 1 # Amperes + UnitTemp_in_cgs: 1 # Kelvin + +# Cosmological parameters +Cosmology: + h: 0.5 + a_begin: 0.047619048 # z_ini = 20 + a_end: 1.0 # z_end = 0 + Omega_m: 1.0 + Omega_lambda: 0.0 + Omega_b: 0.1 + +# Parameters governing the time integration +TimeIntegration: + dt_max: 0.01 + dt_min: 1e-10 + +Scheduler: + max_top_level_cells: 16 + cell_split_size: 100 + +# Parameters governing the snapshots +Snapshots: + basename: santabarbara_low + scale_factor_first: 0.05 + delta_time: 1.02 + +# Parameters governing the conserved quantities statistics +Statistics: + delta_time: 1.02 + scale_factor_first: 0.05 + +# Parameters for the self-gravity scheme +Gravity: + eta: 0.025 + theta: 0.5 + comoving_softening: 0.02 # 20 kpc = 1/25 mean inter-particle separation + max_physical_softening: 0.00526 # 20 ckpc = 5.26 pkpc at z=2.8 (EAGLE-like evolution of softening). + mesh_side_length: 64 + +# Parameters of the hydro scheme +SPH: + resolution_eta: 1.2348 # "48 Ngb" with the cubic spline kernel + CFL_condition: 0.1 + initial_temperature: 1200. # (1 + z_ini)^2 * 2.72K + minimal_temperature: 100. + +# Parameters related to the initial conditions +InitialConditions: + file_name: ./SantaBarbara_128.hdf5 + periodic: 1 + cleanup_h_factors: 1 # ICs were generated for Gadget, we need to get rid of h-factors + cleanup_velocity_factors: 1 # ICs were generated for Gadget, we need to get rid of sqrt(a) factors in the velocity + generate_gas_in_ics: 1 # Generate gas particles from the DM-only ICs + cleanup_smoothing_lengths: 1 # Since we generate gas, make use of the (expensive) cleaning-up procedure. \ No newline at end of file