Import fixes to parameter files for low-z EAGLE examples from the no_marktasks...

Import fixes to parameter files for low-z EAGLE examples from the no_marktasks branch. Also adds an EAGLE-12 RT example

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).

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.]

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).

examples/EAGLE_low_z/EAGLE_12/eagle_12_rt_test.yml

+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)
+
+

examples/EAGLE_low_z/EAGLE_12/run_rt_test.sh

+#!/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
+

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).

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).

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).