MetaData: run_name: RT Cooling Test # Define the system of units to use internally. InternalUnitSystem: UnitMass_in_cgs: 1.98848e43 # 10^10 M_sun in grams UnitLength_in_cgs: 3.08567758e21 # 1 kpc in cm UnitVelocity_in_cgs: 1e5 # 1 km/s in cm/s UnitCurrent_in_cgs: 1 # Amperes UnitTemp_in_cgs: 1 # Kelvin # Parameters governing the time integration TimeIntegration: max_nr_rt_subcycles: 1 time_begin: 0. # The starting time of the simulation (in internal units). time_end: 0.100 # The end time of the simulation (in internal units). dt_min: 1.e-8 # The minimal time-step size of the simulation (in internal units). dt_max: 4.882814e-03 # The maximal time-step size of the simulation (in internal units). # Parameters governing the snapshots Snapshots: basename: output # Common part of the name of output files output_list_on: 0 output_list: snaplist_cooling scale_factor_first: 0.047601 # Time of the first output (in internal units) delta_time: 1.006 # Parameters governing the conserved quantities statistics Statistics: scale_factor_first: 0.047601 delta_time: 1.006 # Time between statistics output # 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). CFL_condition: 0.6 # Courant-Friedrich-Levy condition for time integration. minimal_temperature: 0. # Kelvin # Parameters related to the initial conditions InitialConditions: file_name: ./cooling_test.hdf5 # The file to read periodic: 1 # periodic ICs GEARRT: f_reduce_c: 1.e-9 # This test is without actual radiation, so we don't care about this CFL_condition: 0.9 # CFL condition for RT, independent of hydro photon_groups_Hz: [3.288e15, 5.945e15, 13.157e15] # Lower photon frequency group bin edges in Hz. Needs to have exactly N elements, where N is the configured number of bins --with-RT=GEAR_N stellar_luminosity_model: const # Which model to use to determine the stellar luminosities. const_stellar_luminosities_LSol: [0., 0., 0.] # (Conditional) constant star luminosities for each photon frequency group to use if stellar_luminosity_model:const is set, in units of Solar Luminosity. hydrogen_mass_fraction: 0.76 # total hydrogen mass fraction stellar_spectrum_type: 0 # 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_const_max_frequency_Hz: 1.e17 # (Conditional) if stellar_spectrum_type=0, use this maximal frequency for the constant spectrum. case_B_recombination: 0 # reference solution assumes case A recombination GrackleCooling: cloudy_table: CloudyData_UVB=HM2012.h5 # Name of the Cloudy Table (available on the grackle bitbucket repository) with_UV_background: 0 # Enable or not the UV background redshift: 0 # Redshift to use (-1 means time based redshift) with_metal_cooling: 0 # Enable or not the metal cooling provide_volumetric_heating_rates: 0 # (optional) User provide volumetric heating rates provide_specific_heating_rates: 0 # (optional) User provide specific heating rates max_steps: 10000 # (optional) Max number of step when computing the initial composition convergence_limit: 1 # (optional) Convergence threshold (relative) for initial composition self_shielding_method: 0 # (optional) Grackle (1->3 for Grackle's ones, 0 for none and -1 for GEAR) primordial_chemistry: 1 thermal_time_myr: 5 maximal_density_Hpcm3: -1 # Maximal density (in hydrogen atoms/cm^3) for cooling. Higher densities are floored to this value to ensure grackle works properly when interpolating beyond the cloudy_table maximal density. A value < 0 deactivates this parameter. Scheduler: tasks_per_cell: 128 Cosmology: # Planck13 (EAGLE flavour) a_begin: 0.0476 # z~20 a_end: 0.2 # z~4 h: 0.6777 Omega_cdm: 0.2587481 Omega_lambda: 0.693 Omega_b: 0.0482519