diff --git a/examples/EAGLE_low_z/EAGLE_100/eagle_100.yml b/examples/EAGLE_low_z/EAGLE_100/eagle_100.yml index 14a42441f33462ab3ca698b35d69074010b11309..002c617461af025b725fe6576d192542ea9a63dd 100644 --- a/examples/EAGLE_low_z/EAGLE_100/eagle_100.yml +++ b/examples/EAGLE_low_z/EAGLE_100/eagle_100.yml @@ -47,7 +47,7 @@ Statistics: Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. MAC: adaptive - epsilon_fmm: 0.01 + epsilon_fmm: 0.001 theta_cr: 0.7 # Opening angle (Multipole acceptance criterion) mesh_side_length: 1024 comoving_DM_softening: 0.0026994 # Comoving DM softening length (in internal units). diff --git a/examples/EAGLE_low_z/EAGLE_12/README b/examples/EAGLE_low_z/EAGLE_12/README index 0530021202191bac74690e47106ca3488011d1b5..2e9a053c20ea6859ea3ce0869797f6aabfe25213 100644 --- a/examples/EAGLE_low_z/EAGLE_12/README +++ b/examples/EAGLE_low_z/EAGLE_12/README @@ -14,3 +14,21 @@ running these ICs on 8 cores. MD5 checksum of the ICs: 6f96624175df28f2c4a02265b871b3d5 EAGLE_ICs_12.hdf5 + + +A second example is set up to run with `./run_rt_test.sh`. It makes use +of the `eagle_12_rt_test.yml` file, and is intended to stress-test the +RT implementation rather than produce actual physically meaningful results. + + +o To use GEAR-RT, configure SWIFT with + + --with-stars=basic --with-hydro=gizmo-mfv --with-riemann-solver=hllc --with-rt=GEAR_1 --with-rt-riemann-solver=GLF --with-feedback=none + +[technically, any other feedback scheme should work as well.] + + +o To use the DEBUG RT scheme, configure SWIFT with + --with-stars=basic --with-rt=debug --with-feedback=none + +[technically, any other feedback scheme should work as well.] diff --git a/examples/EAGLE_low_z/EAGLE_12/eagle_12.yml b/examples/EAGLE_low_z/EAGLE_12/eagle_12.yml index 24e6c7a53aae26663049cf838a930b626d6ab4f6..4c0a6496de31a30f5b532c98c5af526a159d9776 100644 --- a/examples/EAGLE_low_z/EAGLE_12/eagle_12.yml +++ b/examples/EAGLE_low_z/EAGLE_12/eagle_12.yml @@ -48,7 +48,7 @@ Statistics: Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. MAC: adaptive - epsilon_fmm: 0.01 + epsilon_fmm: 0.001 theta_cr: 0.7 # Opening angle (Multipole acceptance criterion) mesh_side_length: 128 comoving_DM_softening: 0.0026994 # Comoving DM softening length (in internal units). diff --git a/examples/EAGLE_low_z/EAGLE_12/eagle_12_rt_test.yml b/examples/EAGLE_low_z/EAGLE_12/eagle_12_rt_test.yml new file mode 100644 index 0000000000000000000000000000000000000000..bd5c4139fc2a52cd57d41c6c61f7906029262ba6 --- /dev/null +++ b/examples/EAGLE_low_z/EAGLE_12/eagle_12_rt_test.yml @@ -0,0 +1,102 @@ +MetaData: + run_name: EAGLE-L0012N0188-Ref-RT-TEST + +# Define the system of units to use internally. +InternalUnitSystem: + UnitMass_in_cgs: 1.98841e43 # 10^10 M_sun in grams + UnitLength_in_cgs: 3.08567758e24 # Mpc in centimeters + UnitVelocity_in_cgs: 1e5 # km/s in centimeters per second + UnitCurrent_in_cgs: 1 # Amperes + UnitTemp_in_cgs: 1 # Kelvin + +# Cosmological parameters +Cosmology: + h: 0.6777 # Reduced Hubble constant + a_begin: 0.9090909 # Initial scale-factor of the simulation + a_end: 1.0 # Final scale factor of the simulation + Omega_cdm: 0.2587481 # Cold Dark Matter density parameter + Omega_lambda: 0.693 # Dark-energy density parameter + Omega_b: 0.0482519 # Baryon density parameter + +# Parameters governing the time integration +TimeIntegration: + max_nr_rt_subcycles: 32 + time_begin: 0. # The starting time of the simulation (in internal units). + time_end: 1e-2 # The end time of the simulation (in internal units). + dt_min: 1e-12 # The minimal time-step size of the simulation (in internal units). + dt_max: 1e-2 # The maximal time-step size of the simulation (in internal units). + +# Parameters governing the snapshots +Snapshots: + basename: eagle # Common part of the name of output files + scale_factor_first: 0.91 # Scale-factor of the first snaphot (cosmological run) + time_first: 0.01 # Time of the first output (non-cosmological run) (in internal units) + delta_time: 1.01 # Time difference between consecutive outputs (in internal units) + compression: 1 + +Scheduler: + links_per_tasks: 500 + dependency_graph_frequency: 0 + +# Parameters governing the conserved quantities statistics +Statistics: + scale_factor_first: 0.92 # Scale-factor of the first stat dump (cosmological run) + time_first: 0.01 # Time of the first stat dump (non-cosmological run) (in internal units) + delta_time: 1.05 # Time between statistics output + +# Parameters for the self-gravity scheme +Gravity: + eta: 0.025 # Constant dimensionless multiplier for time integration. + MAC: adaptive + epsilon_fmm: 0.001 + theta_cr: 0.7 # Opening angle (Multipole acceptance criterion) + mesh_side_length: 128 + comoving_DM_softening: 0.0026994 # Comoving DM softening length (in internal units). + max_physical_DM_softening: 0.0007 # Max physical DM softening length (in internal units). + comoving_baryon_softening: 0.0026994 # Comoving DM softening length (in internal units). + max_physical_baryon_softening: 0.0007 # Max physical DM softening length (in internal units). + +# Parameters for the hydrodynamics scheme +SPH: + resolution_eta: 1.2348 # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel). + h_min_ratio: 0.1 # Minimal smoothing length in units of softening. + h_max: 0.5 # Maximal smoothing length in co-moving internal units. + CFL_condition: 0.1 # Courant-Friedrich-Levy condition for time integration. + minimal_temperature: 100 # (internal units) + particle_splitting: 0 + particle_splitting_mass_threshold: 7e-4 # Internal units (i.e. 7e6 Msun ~ 4 times the initial gas mass) + +# Parameters of the stars neighbour search +Stars: + resolution_eta: 1.1642 # Target smoothing length in units of the mean inter-particle separation + h_tolerance: 7e-3 + overwrite_birth_time: 1 + birth_time: 0.33333 # Pretend all the stars were born at z = 2 + +# Parameters related to the initial conditions +InitialConditions: + file_name: ./EAGLE_ICs_12.hdf5 # The file to read + periodic: 1 + cleanup_h_factors: 1 # Remove the h-factors inherited from Gadget + cleanup_velocity_factors: 1 # Remove the sqrt(a) factor in the velocities inherited from Gadget + +GEARRT: + f_reduce_c: 0.01 # reduce the speed of light for the RT solver by multiplying c with this factor + CFL_condition: 0.9 + f_limit_cooling_time: 0.1 # (Optional) multiply the cooling time by this factor when estimating maximal next time step. Set to 0.0 to turn computation of cooling time off. + photon_groups_Hz: [3.288e15, 5.944e15, 13.157e15] # Photon frequency group bin edges in Hz. Needs to be 1 less than the number of groups (N) requested during the configuration (--with-RT=SPHM1RT_N). Outer edges of zero and infinity are assumed. + stellar_luminosity_model: const # Which model to use to determine the stellar luminosities. + const_stellar_luminosities_LSol: [1.764e+04, 3.631e+04, 8.037e+03 ] # (Conditional) constant star emission rates for each photon frequency group to use if use_constant_emission_rates is set, in units of Solar Luminosity. + hydrogen_mass_fraction: 0.75 # total hydrogen (H + H+) mass fraction in the metal-free portion of the gas + set_equilibrium_initial_ionization_mass_fractions: 0 # (Optional) set the initial ionization fractions depending on gas temperature assuming ionization equilibrium. + set_initial_ionization_mass_fractions: 1 # (Optional) manually overwrite initial mass fraction of each species (using the values you set below) + mass_fraction_HI: 0.75 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HI mass fractions to this value + mass_fraction_HII: 1.e-6 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HII mass fractions to this value + mass_fraction_HeI: 0.25 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HeI mass fractions to this value + mass_fraction_HeII: 1.e-6 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HeII mass fractions to this value + mass_fraction_HeIII: 1.e-6 # (Conditional) If set_initial_ionization_fractions=1, needed to set initial HeIII mass fractions to this value + stellar_spectrum_type: 1 # Which radiation spectrum to use. 0: constant from 0 until some max frequency set by stellar_spectrum_const_max_frequency_Hz. 1: blackbody spectrum. + stellar_spectrum_blackbody_temperature_K: 1.e5 # (Conditional) if stellar_spectrum_type=1, use this temperature (in K) for the blackbody spectrum. + stars_max_timestep: 5.468750e-05 # (Optional) restrict the maximal timestep of stars to this value (in internal units) + + diff --git a/examples/EAGLE_low_z/EAGLE_12/run_rt_test.sh b/examples/EAGLE_low_z/EAGLE_12/run_rt_test.sh new file mode 100755 index 0000000000000000000000000000000000000000..eed49dc5227a7615fa30556d1e26caa93c3c6524 --- /dev/null +++ b/examples/EAGLE_low_z/EAGLE_12/run_rt_test.sh @@ -0,0 +1,27 @@ +#!/bin/bash + +# Run a RT test with EAGLE ICs. +# +# To use GEAR-RT, configure SWIFT with +# +# --with-stars=basic --with-hydro=gizmo-mfv --with-riemann-solver=hllc --with-rt=GEAR_1 --with-rt-riemann-solver=GLF --with-feedback=none +# [technically, any other feedback scheme should work as well.] +# +# +# To use the DEBUG RT scheme, configure SWIFT with +# --with-stars=basic --with-rt=debug --with-feedback=none +# [technically, any other feedback scheme should work as well.] + + +# Generate the initial conditions if they are not present. +if [ ! -e EAGLE_ICs_12.hdf5 ] +then + echo "Fetching initial conditions for the EAGLE 12Mpc example..." + ./getIC.sh +fi + +../../../swift \ + --hydro --threads=8 --stars --self-gravity \ + --feedback --radiation \ + eagle_12_rt_test.yml + diff --git a/examples/EAGLE_low_z/EAGLE_25/eagle_25.yml b/examples/EAGLE_low_z/EAGLE_25/eagle_25.yml index d79359fc684ce4e042c7827ccceb28c7e3d8db15..03ab1520bebcd1a6c91bd520ad33728551f63462 100644 --- a/examples/EAGLE_low_z/EAGLE_25/eagle_25.yml +++ b/examples/EAGLE_low_z/EAGLE_25/eagle_25.yml @@ -55,7 +55,7 @@ Statistics: Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. MAC: adaptive - epsilon_fmm: 0.01 + epsilon_fmm: 0.001 theta_cr: 0.7 # Opening angle (Multipole acceptance criterion) mesh_side_length: 256 comoving_DM_softening: 0.0026994 # Comoving DM softening length (in internal units). diff --git a/examples/EAGLE_low_z/EAGLE_50/eagle_50.yml b/examples/EAGLE_low_z/EAGLE_50/eagle_50.yml index c84621a62512f0429defe6bbbff18abd4296e281..d93d2d8cefdf1e5f92ecbffe81d91daff021bb25 100644 --- a/examples/EAGLE_low_z/EAGLE_50/eagle_50.yml +++ b/examples/EAGLE_low_z/EAGLE_50/eagle_50.yml @@ -47,7 +47,7 @@ Statistics: Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. MAC: adaptive - epsilon_fmm: 0.01 + epsilon_fmm: 0.001 theta_cr: 0.7 # Opening angle (Multipole acceptance criterion) mesh_side_length: 512 comoving_DM_softening: 0.0026994 # Comoving DM softening length (in internal units). diff --git a/examples/EAGLE_low_z/EAGLE_6/eagle_6.yml b/examples/EAGLE_low_z/EAGLE_6/eagle_6.yml index bc1b78866b0785edbb12871285f192989480c159..963821e177f8baec2ef6c78f58552b213be16fe3 100644 --- a/examples/EAGLE_low_z/EAGLE_6/eagle_6.yml +++ b/examples/EAGLE_low_z/EAGLE_6/eagle_6.yml @@ -58,7 +58,7 @@ Statistics: Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. MAC: adaptive - epsilon_fmm: 0.01 + epsilon_fmm: 0.001 theta_cr: 0.7 # Opening angle (Multipole acceptance criterion) mesh_side_length: 64 comoving_DM_softening: 0.0026994 # Comoving DM softening length (in internal units).