Main Code

ASLOTH main program

program  asloth
Use :

config, defined_types, eps_wrapper, tree_memory_arrays_passable, numerical_parameters, resolution_parameters, to_file, input_files, check_flags, utility, volume_fractions, feedback_routines (node_feedback()), bubble_tree (start_btree()), tree_reader, tree_initialization, star_formation, read_yields, metal_functions, populations, sf_routines, crit_masses (set_mcrit()), cosmic_time (set_cosmic_time())

Call to:

set_default_config(), read_config(), check_compiler_flags(), assert(), init_outputs(), set_cosmic_time(), set_mcrit(), read_tree(), make_eps_tree(), init_stellar_populations(), assign_yields(), output_pop(), start_btree(), init_tree(), print_bar(), real_to_string3(), mem_report(), clean_tree(), get_fracs(), node_feedback(), sf_step(), add_parent_to_tree(), add_to_sfr(), add_to_feedback_volumes(), add_to_base_mdf(), add_to_base_mstariisurv(), write_files_reali(), outputs_satellites(), write_mw_properties(), write_files(), close_outputs()

Configuration

config

Quick access

Variables:

config_file, cubic_box, nthreads

Routines:

read_config(), set_default_config()

Needed modules

  • numerical_parameters

  • input_files: List of input files

  • random_object: This is an object oriented version of the ran3 algorithm from numerical recipes Before using a rng type object it needs to be initialized by calling rng%init(seed) where seed is an integer seed Afterwards rng%get_ran() will return a random real in the range [0,1)

Variables

  • config/config_file [character,protected]

    e.g. param_template.nml

  • config/cubic_box [logical,protected]
  • config/nthreads [integer,protected]

Subroutines and functions

subroutine  config/set_default_config()
Called from:

asloth

subroutine  config/read_config(fname)
Parameters:

fname [character,in]

Use :

omp_lib

Called from:

asloth

Type definitions

defined_types

Quick access

Types:

node_arrays, stars_highmass, treenode

Routines:

add_baryons(), del_high_mass(), make_stars_high_mass(), print_all_props(), set_baryons()

Needed modules

Types

  • type  defined_types/treenode

    class for storing information of each halo at each time

    Type fields:

    • % base [treenode,pointer] :: which halo will this one be at at the the lowest redshift

    • % child [treenode,pointer/optional/default=>null()]

    • % desc_id [integer] :: descendent id

    • % firstpopiisf [logical,optional/default=.true.]

    • % highstars [stars_highmass,pointer/optional/default=>]

    • % highstars_first [stars_highmass,pointer/optional/default=>]

    • % highstars_last [stars_highmass,pointer/optional/default=>]

    • % i_checked [logical] :: has this halo already been checked for stochastic feedback?

    • % i_enr [logical,optional/default=.false.] :: is this halo enriched?

    • % i_ion [logical,optional/default=.false.] :: is the halo affected by external ionising radiation

    • % i_mcrit [logical,optional/default=.false.] :: is the halo above the critical mass?

    • % i_sf [logical,optional/default=.false.] :: this flag controls whether the SF_step subroutine is called

    • % i_sub [logical,optional/default=.false.] :: is this a subhalo?

    • % id [integer] :: halo id

    • % igm_z (*) [real,allocatable] :: ambient metallicity

    • % jlevel [integer] :: which time-step is the halo at

    • % l_ion [real] :: ionising radiation emission rate [Photons/second]

    • % llstars [stars_highmass,pointer/optional/default=>]

    • % llstars_first [stars_highmass,pointer/optional/default=>]

    • % llstars_last [stars_highmass,pointer/optional/default=>]

    • % lp_id [integer] :: last progenitor depthfirst id

    • % m_cold [real] :: cold gas mass [M_sun]

    • % m_filter [real,optional/default=0]

    • % m_hot [real] :: hot gas mass [M_sun]

    • % m_metals (,) [real,allocatable]

    • % m_out [real] :: outflow mass [M_sun]

    • % m_peak [real] :: peak mass [Msun]

    • % m_peak_filter [real,optional/default=0]

    • % m_starii [real] :: stellar mass Pop II [M_sun]

    • % m_starii_surv [real,optional/default=0.0] :: Msun

    • % m_stariii [real] :: stellar mass Pop III [M_sun]

    • % mhalo [real] :: mass of halo

    • % nchild [integer] :: number of child nodes

    • % null [node_arrays,pointer]

    • % num_popiiprogenitor [integer,optional/default=0] :: how many Pop II forming progenitors

    • % parent [treenode,pointer/optional/default=>null()]

    • % pop [integer]

    • % r_en [real,optional/default=0] :: enrichment radius [physical Mpc]

    • % r_ion [real,optional/default=0] :: ionised radius [physical Mpc]

    • % sibling [treenode,pointer/optional/default=>null()]

    • % this_array [node_arrays,pointer/optional/default=>]

    • % x (3) [real] :: position

  • type  defined_types/node_arrays
    Type fields:

    • % mdf_array (,) [real,allocatable]

  • type  defined_types/stars_highmass
    Type fields:

    • % bin_id [integer,optional/default=0] :: which bin it belongs, should be consistent with IMF arrays

    • % next [stars_highmass,pointer/optional/default=>]

    • % null [stars_highmass,pointer]

    • % num [integer,optional/default=0]

    • % pop [integer,optional/default=0] :: PopII (2) or PopIII (3) stars

    • % prev [stars_highmass,pointer/optional/default=>]

    • % timeborn [real,optional/default=0.0] :: in order to know when the stars die

    • % z (*) [real,allocatable] :: Zgas [Z/Zsun] when the star is born. Related with SN yields.

Subroutines and functions

function  defined_types/make_stars_high_mass(n, i, timeborn, zgas, pop)
Parameters:

  • n [integer,in]

  • i [integer,in]

  • timeborn [real,in]

  • zgas (*) [real,in]

  • pop [integer,in]

Return:

res [stars_highmass]

subroutine  defined_types/del_high_mass(this)
Parameters:

this [stars_highmass,inout]

subroutine  defined_types/set_baryons(this, other)

copies baryonic properties from one halo to another (overwrites)

Parameters:

  • this [real]

  • other [real]

subroutine  defined_types/add_baryons(this, other)

adds baryonic properties of one halo to another

Parameters:

  • this [real]

  • other [real]

subroutine  defined_types/print_all_props(a)

Printing function that prints various porperties of a halo in order to allow easier tracing of errors and unexpected changes

Parameters:

a [treenode,in]

Tree Reader

tree_reader

Quick access

Routines:

read_tree()

Needed modules

Subroutines and functions

subroutine  tree_reader/read_tree(count, read_mass)
Parameters:

  • count [integer,out] :: reading tree

  • read_mass [real,out]

Called from:

asloth

Call to:

print_bar(), mem_report()

Tree initialization

tree_initialization

Description

prepares tree for running asloth

Quick access

Routines:

filter_virial_masses(), init_tree(), set_tree_params()

Subroutines and functions

subroutine  tree_initialization/init_tree(count, n_jcurrent)
Parameters:

  • count [integer,in]

  • n_jcurrent (nlev_max) [integer,out]

Use :

defined_types, numerical_parameters, tree_memory_arrays_passable, utility, crit_masses, cosmic_time, resolution_parameters, metalmix_dz, filter_module, trace_mdf, metals

Called from:

asloth

Call to:

print_bar(), mem_report(), read_lookups(), filter_virial_masses(), init_mdf()

subroutine  tree_initialization/filter_virial_masses()
Use :

defined_types, numerical_parameters, tree_memory_arrays_passable, filter_module

Called from:

init_tree()

Call to:

compute_savgol_coeff(), print_bar(), savgol_filter(), savgol_free()

subroutine  tree_initialization/set_tree_params()
Use :

power_spectrum_parameters, modified_merger_tree, numerical_parameters, input_files, time_parameters

Called from:

make_eps_tree()

numerical and resolution parameters

numerical parameters

Quick access

Variables:

a_radiation, ab_mag_norm, ab_norm, alpha_b, alpha_ion, alpha_outflow, ang2m, atomic_mass_helium, atomic_mass_hydrogen, bt_delta, bt_max, bt_min, btnbin, btnpad, c_cgs, c_hii, c_light, c_light_angstroms, c_light_cm, clump, cm32m3, crit_freq, deca, delta200, deltaeds, e_51, em_rdisk_half_scale_o2, eps3, eps4, eps6, epsm, epsr, erg2j, escape_fraction, etaii, etaiii, ev2erg, ev2j, f_escii, f_esciii, fx_peak_nfw, g, g_cgs, g_gyrkms3, g_mpcgyr, g_mpcgyr2, g_si, gpi, gyr2s, gyr2yr, h0100, h0100pgyr, h_planck, h_planck_erg, hecto, hubble, hubble_cgs, imf_max, imf_min, iso_fac, j2erg, k_boltzmann, k_boltzmann_erg, khorizon, kilo, kind99, km2m, kms2mpcgyr, kom, kpc2cm, kstrc_1, kstrc_2, l_ab0, l_box, ln10, ln2, log4, log_lsun_erg, logang2m, logc_light, logc_light_angstroms, logdeca, logev2erg, loggyr2s, loghecto, logj2erg, logk_boltzmann, logkilo, logl_ab0, loglsun, loglsun_bc, logm2cm, logmega, logmpc2asat10pc, logpc2m, logpi, lograd2as, long_bn, lsun, lsun40_bc, lsun_bc, lsun_erg, lw_mode, m2cm, m32cm3, m8crit, m_atomic, m_atomic_g, m_be, m_cha_out, m_electron, mean_mol, mega, mend, milli, mp_cgs, mpc2asat10pc, mpc2cm, mpc2m, mpc_cgs, mpckm2gyr, mres, msolar, msolar_1030kg, msolar_g, mstart, msun_cgs, msun_mpc3_cgs, mu_primordial, myr2s, n_b_cgs, n_cloud_default, n_cold_default, n_h_cgs, n_ionii, n_ioniii, n_spec, nextoutputid, nlev, nlev_max, nstarmassbinii, nstarmassbiniii, omega_b, omega_l, omega_m, pc2cm, pc2m, pi, pi4, pio2, pio4, pisq, popii_slope, rad2as, rad2degrees, rdisk_half_scale, rho_b_ast, rho_b_cgs, rho_crit_cgs, rho_m_cgs, rhocrit, right_angle_degrees, semi_circle_degrees, sigma_8, sigma_t, sigma_thomson, slope, sqrt2, sqrt2opi, sqrt2pi, sqrtatomic_mass_hydrogen, sqrtdelta200, sqrtg, sqrtg_si, sqrtk_boltzmann, sqrtm_atomic, sqrtmpc2m, sqrtmsolar, sqrtpi, sqrtrhocrit, stdout_bn, surf_dens_norm_half, t_crit, t_ion, t_univ, tree_mass, v_com, v_enrich_popiii, v_out, vbc, x_hydrogen, x_hydrogen_solar, xh, y_helium, y_helium_solar, year_cgs, yhe, z_crit, z_metals, z_metals_solar

Variables

  • numerical_parameters/a_radiation [real,parameter=7.565743419463647e-16]
  • numerical_parameters/ab_mag_norm [real,parameter=48.6]
  • numerical_parameters/ab_norm [real,parameter=3.6307805477e-20]
  • numerical_parameters/alpha_b [real,parameter=2.6e-13]
  • numerical_parameters/alpha_ion [real,parameter=2.6e-13]
  • numerical_parameters/alpha_outflow [real]
  • numerical_parameters/ang2m [real,parameter=1e-10]
  • numerical_parameters/atomic_mass_helium [real,parameter=4.002602]
  • numerical_parameters/atomic_mass_hydrogen [real,parameter=1.00794]
  • numerical_parameters/bt_delta [real,optional/default=0.1]
  • numerical_parameters/bt_max [real,optional/default=log10(4.e17)]
  • numerical_parameters/bt_min [real,optional/default=log10(4.e15)]
  • numerical_parameters/btnbin [integer,optional/default=10]
  • numerical_parameters/btnpad [integer,optional/default=2]
  • numerical_parameters/c_cgs [real,parameter=29980000000.0]
  • numerical_parameters/c_hii [real,parameter=3]
  • numerical_parameters/c_light [real,parameter=299792458.0]
  • numerical_parameters/c_light_angstroms [real,parameter=2.99792458e+18]
  • numerical_parameters/c_light_cm [real,parameter=29979245800.0]
  • numerical_parameters/clump [real,parameter=4.0]
  • numerical_parameters/cm32m3 [real,parameter=1e-06]
  • numerical_parameters/crit_freq [real,parameter=3.969127815104263e-18]
  • numerical_parameters/deca [real,parameter=10.0]
  • numerical_parameters/delta200 [real,parameter=200.0]
  • numerical_parameters/deltaeds [real,parameter=177.65287922076283]
  • numerical_parameters/e_51 [real,parameter=1.2]
  • numerical_parameters/em_rdisk_half_scale_o2 [real,parameter=0.432067481]
  • numerical_parameters/eps3 [real,parameter=0.001]
  • numerical_parameters/eps4 [real,parameter=0.0001]
  • numerical_parameters/eps6 [real,parameter=1e-06]
  • numerical_parameters/epsm [real,parameter=0.0001]
  • numerical_parameters/epsr [real,parameter=1e-06]
  • numerical_parameters/erg2j [real,parameter=1e-07]
  • numerical_parameters/escape_fraction [real]
  • numerical_parameters/etaii [real]
  • numerical_parameters/etaiii [real]
  • numerical_parameters/ev2erg [real,parameter=1.60217733e-12]
  • numerical_parameters/ev2j [real,parameter=1.60217733e-19]
  • numerical_parameters/f_escii [real]
  • numerical_parameters/f_esciii [real]
  • numerical_parameters/fx_peak_nfw [real,parameter=0.2162165956]
  • numerical_parameters/g [real,parameter=4.301307710181788e-09]
  • numerical_parameters/g_cgs [real,parameter=6.674e-08]
  • numerical_parameters/g_gyrkms3 [real,parameter=4.205785220992724e-06]
  • numerical_parameters/g_mpcgyr [real,parameter=4.398999721935946e-12]
  • numerical_parameters/g_mpcgyr2 [real,parameter=4.498910530810333e-15]
  • numerical_parameters/g_si [real,parameter=6.67259e-11]
  • numerical_parameters/gpi [real,parameter=1.3212863952890194e-05]
  • numerical_parameters/gyr2s [real,parameter=3.15576e+16]
  • numerical_parameters/gyr2yr [real,parameter=1000000000.0]
  • numerical_parameters/h0100 [real,parameter=100.0]
  • numerical_parameters/h0100pgyr [real,parameter=0.10227121653079844]
  • numerical_parameters/h_planck [real,parameter=6.6260755e-34]
  • numerical_parameters/h_planck_erg [real,parameter=6.6260755e-27]
  • numerical_parameters/hecto [real,parameter=100.0]
  • numerical_parameters/hubble [real,parameter=0.6774]
  • numerical_parameters/hubble_cgs [real,parameter=2.1952879411478756e-18]
  • numerical_parameters/imf_max [real]
  • numerical_parameters/imf_min [real]
  • numerical_parameters/iso_fac [real,parameter=1.085736204755322]
  • numerical_parameters/j2erg [real,parameter=10000000.0]
  • numerical_parameters/k_boltzmann [real,parameter=1.3806503e-23]
  • numerical_parameters/k_boltzmann_erg [real,parameter=1.3806503e-16]
  • numerical_parameters/khorizon [real,parameter=0.00033356409519815205]
  • numerical_parameters/kilo [real,parameter=1000.0]
  • numerical_parameters/kind99 [integer,parameter=selected_real_kind(p,r=99)]
  • numerical_parameters/km2m [real,parameter=1000.0]
  • numerical_parameters/kms2mpcgyr [real,parameter=0.0010227121653079844]
  • numerical_parameters/kom [real,parameter=2.9434011554855882e-21]
  • numerical_parameters/kpc2cm [real,parameter=3.0856775807e+21]
  • numerical_parameters/kstrc_1 [real,parameter=0.6251543028]
  • numerical_parameters/kstrc_2 [real,parameter=1.191648617]
  • numerical_parameters/l_ab0 [real,parameter=43442114331434.14]
  • numerical_parameters/l_box [real]

    (linear box size)

  • numerical_parameters/ln10 [real,parameter=2.302585093]
  • numerical_parameters/ln2 [real,parameter=0.6931471806]
  • numerical_parameters/log4 [real,parameter=0.602059991]
  • numerical_parameters/log_lsun_erg [real,parameter=33.5848963]
  • numerical_parameters/logang2m [real,parameter=-10.0]
  • numerical_parameters/logc_light [real,parameter=8.4768207]
  • numerical_parameters/logc_light_angstroms [real,parameter=18.476820699999998]
  • numerical_parameters/logdeca [real,parameter=1.0]
  • numerical_parameters/logev2erg [real,parameter=-11.7952894176]
  • numerical_parameters/loggyr2s [real,parameter=16.499103967]
  • numerical_parameters/loghecto [real,parameter=2.0]
  • numerical_parameters/logj2erg [real,parameter=7.0]
  • numerical_parameters/logk_boltzmann [real,parameter=-22.859916308]
  • numerical_parameters/logkilo [real,parameter=3.0]
  • numerical_parameters/logl_ab0 [real,parameter=13.637910954000002]
  • numerical_parameters/loglsun [real,parameter=26.5848963]
  • numerical_parameters/loglsun_bc [real,parameter=26.5868122]
  • numerical_parameters/logm2cm [real,parameter=2.0]
  • numerical_parameters/logmega [real,parameter=6.0]
  • numerical_parameters/logmpc2asat10pc [real,parameter=10.31442513287]
  • numerical_parameters/logpc2m [real,parameter=16.489350545]
  • numerical_parameters/logpi [real,parameter=0.497149873]
  • numerical_parameters/lograd2as [real,parameter=5.31442513287]
  • numerical_parameters/long_bn [integer,parameter=selected_int_kind(16)]
  • numerical_parameters/lsun [real,parameter=3.845e+26]
  • numerical_parameters/lsun40_bc [real,parameter=3.862e-07]
  • numerical_parameters/lsun_bc [real,parameter=3.862e+26]
  • numerical_parameters/lsun_erg [real,parameter=3.8450000000000004e+33]
  • numerical_parameters/lw_mode [integer]
  • numerical_parameters/m2cm [real,parameter=100.0]
  • numerical_parameters/m32cm3 [real,parameter=1000000.0]
  • numerical_parameters/m8crit [real,parameter=595167835572453.2]
  • numerical_parameters/m_atomic [real,parameter=1.66053873e-27]
  • numerical_parameters/m_atomic_g [real,parameter=1.66053873e-24]
  • numerical_parameters/m_be [real,parameter=1000.0]
  • numerical_parameters/m_cha_out [real]
  • numerical_parameters/m_electron [real,parameter=9.10938188e-31]
  • numerical_parameters/mean_mol [real,parameter=1.2281617489023304]
  • numerical_parameters/mega [real,parameter=1000000.0]
  • numerical_parameters/mend [real,parameter=2.0]

    , Mcutoff !, MLongLive

  • numerical_parameters/milli [real,parameter=0.001]
  • numerical_parameters/mp_cgs [real,parameter=1.673e-24]
  • numerical_parameters/mpc2asat10pc [real,parameter=20626480624.64262]
  • numerical_parameters/mpc2cm [real,parameter=3.0856775807e+24]
  • numerical_parameters/mpc2m [real,parameter=3.0856775807e+22]
  • numerical_parameters/mpc_cgs [real,parameter=3.0857e+24]
  • numerical_parameters/mpckm2gyr [real,parameter=977.79222143002]
  • numerical_parameters/mres [real,parameter=500000.0]
  • numerical_parameters/msolar [real,parameter=1.9891e+30]
  • numerical_parameters/msolar_1030kg [real,parameter=1.9890999999999999]
  • numerical_parameters/msolar_g [real,parameter=1.9891e+33]
  • numerical_parameters/mstart [real,parameter=-2.0]
  • numerical_parameters/msun_cgs [real,parameter=1.989e+33]
  • numerical_parameters/msun_mpc3_cgs [real,parameter=6.769766377778569e-41]
  • numerical_parameters/mu_primordial [real,parameter=0.5847493469151946]
  • numerical_parameters/myr2s [real,parameter=31557600000000.0]
  • numerical_parameters/n_b_cgs [real,parameter=2.038744909036926e-07]
  • numerical_parameters/n_cloud_default [real,parameter=1000]

    ! dense, and diffuse density in the cold phase

  • numerical_parameters/n_cold_default [real,parameter=100.0]
  • numerical_parameters/n_h_cgs [real,parameter=1.5337477950684792e-07]
  • numerical_parameters/n_ionii [real,parameter=4000.0]
  • numerical_parameters/n_ioniii [real,parameter=90000.0]
  • numerical_parameters/n_spec [real,parameter=0.9667]
  • numerical_parameters/nextoutputid [integer]
  • numerical_parameters/nlev [integer]
  • numerical_parameters/nlev_max [integer,parameter=400]
  • numerical_parameters/nstarmassbinii [integer,parameter=64]
  • numerical_parameters/nstarmassbiniii [integer,parameter=64]
  • numerical_parameters/omega_b [real,parameter=0.0486]
  • numerical_parameters/omega_l [real,parameter=0.6911]
  • numerical_parameters/omega_m [real,parameter=0.3089]
  • numerical_parameters/pc2cm [real,parameter=3.0856775807e+18]
  • numerical_parameters/pc2m [real,parameter=3.0856775807e+16]
  • numerical_parameters/pi [real,parameter=3.1415926536]
  • numerical_parameters/pi4 [real,parameter=12.5663706144]
  • numerical_parameters/pio2 [real,parameter=1.5707963268]
  • numerical_parameters/pio4 [real,parameter=0.7853981634]
  • numerical_parameters/pisq [real,parameter=9.86960440115349]
  • numerical_parameters/popii_slope [real,parameter=2.3]
  • numerical_parameters/rad2as [real,parameter=206264.80624642622]
  • numerical_parameters/rad2degrees [real,parameter=57.29577951289617]
  • numerical_parameters/rdisk_half_scale [real,parameter=1.67834699]
  • numerical_parameters/rho_b_ast [real,parameter=6187863374.547873]
  • numerical_parameters/rho_b_cgs [real,parameter=4.189038942330162e-31]
  • numerical_parameters/rho_crit_cgs [real,parameter=8.619421692037372e-30]
  • numerical_parameters/rho_m_cgs [real,parameter=2.6625393606703443e-30]
  • numerical_parameters/rhocrit [real,parameter=277511434571.3454]
  • numerical_parameters/right_angle_degrees [real,parameter=90.0]
  • numerical_parameters/semi_circle_degrees [real,parameter=180.0]
  • numerical_parameters/sigma_8 [real,parameter=0.8159]
  • numerical_parameters/sigma_t [real,parameter=6.65e-25]
  • numerical_parameters/sigma_thomson [real,parameter=6.65245854e-29]
  • numerical_parameters/slope [real]
  • numerical_parameters/sqrt2 [real,parameter=1.4142135624]
  • numerical_parameters/sqrt2opi [real,parameter=0.7978845608]
  • numerical_parameters/sqrt2pi [real,parameter=2.5066282746]
  • numerical_parameters/sqrtatomic_mass_hydrogen [real,parameter=1.00396215]
  • numerical_parameters/sqrtdelta200 [real,parameter=14.142135623]
  • numerical_parameters/sqrtg [real,parameter=6.558435567945511e-05]
  • numerical_parameters/sqrtg_si [real,parameter=8.16859228998e-06]
  • numerical_parameters/sqrtk_boltzmann [real,parameter=3.715710295489e-12]
  • numerical_parameters/sqrtm_atomic [real,parameter=4.07497083425e-14]
  • numerical_parameters/sqrtmpc2m [real,parameter=175660968365.0]
  • numerical_parameters/sqrtmsolar [real,parameter=1410354565300000.0]
  • numerical_parameters/sqrtpi [real,parameter=1.7724538509]
  • numerical_parameters/sqrtrhocrit [real,parameter=526793.5407696816]
  • numerical_parameters/stdout_bn [integer,parameter=10]
  • numerical_parameters/surf_dens_norm_half [real,parameter=0.0297114]
  • numerical_parameters/t_crit [real,parameter=2200.0]
  • numerical_parameters/t_ion [real,parameter=10000.0]
  • numerical_parameters/t_univ [real,parameter=4.3549644e+17]
  • numerical_parameters/tree_mass [real]
  • numerical_parameters/v_com [real]
  • numerical_parameters/v_enrich_popiii [real,parameter=11000000.0]
  • numerical_parameters/v_out [real,parameter=11000000.0]
  • numerical_parameters/vbc [real]
  • numerical_parameters/x_hydrogen [real,parameter=0.778]
  • numerical_parameters/x_hydrogen_solar [real,parameter=0.7381]
  • numerical_parameters/xh [real,parameter=0.7523]
  • numerical_parameters/y_helium [real,parameter=0.222]
  • numerical_parameters/y_helium_solar [real,parameter=0.2485]
  • numerical_parameters/year_cgs [real,parameter=31560000.0]
  • numerical_parameters/yhe [real,parameter=0.2477]
  • numerical_parameters/z_crit [real,parameter=-5.0]
  • numerical_parameters/z_metals [real,parameter=5.36e-10]
  • numerical_parameters/z_metals_solar [real,parameter=0.0134]

resolution parameters

Quick access

Variables:

m_ccsn_max, m_ccsn_min, m_ion, m_pisn_max, m_pisn_min, m_survive, zmax, zmin

Variables

  • resolution_parameters/m_ccsn_max [real,parameter=40]

    Mass range for CCSNe

  • resolution_parameters/m_ccsn_min [real,parameter=10]
  • resolution_parameters/m_ion [real,parameter=5]

    Stellar mass above which we account for ionizing radition

  • resolution_parameters/m_pisn_max [real,parameter=260]

    Mass range for PISNe

  • resolution_parameters/m_pisn_min [real,parameter=140]
  • resolution_parameters/m_survive [real,parameter=0.81]

    Stellar mass with lifetime 13.8Gyr.

  • resolution_parameters/zmax [real,parameter=35.0]
  • resolution_parameters/zmin [real,parameter=0.0]

Precomputed time-dependent quantities

crit_masses

Quick access

Variables:

atomic_cooling_mass, mcrit

Routines:

get_atomic_cooling_mass(), get_crit_mass(), get_m_vir(), get_mcrit_schauer(), get_mcrit_stacy(), get_t_vir(), set_mcrit()

Needed modules

Variables

  • crit_masses/atomic_cooling_mass (nlev_max) [real,protected]
  • crit_masses/mcrit (nlev_max) [real,protected]

Subroutines and functions

subroutine  crit_masses/set_mcrit()
Called from:

asloth

Call to:

get_crit_mass(), get_atomic_cooling_mass()

function  crit_masses/get_crit_mass(z)

Determines critical mass for star formation at current redshift

Parameters:

z [real,in] :: redshift

Return:

get_crit_mass [real]

Called from:

set_mcrit()

Call to:

get_mcrit_schauer(), get_mcrit_stacy(), get_atomic_cooling_mass()

function  crit_masses/get_mcrit_stacy(z)

Modified to also check Jeans Mass calculated with modified speed of sound See Stacy et al. 2011

Parameters:

z [real,in] :: redshift

Return:

get_mcrit_stacy [real]

Called from:

get_crit_mass()

Call to:

get_m_vir()

function  crit_masses/get_mcrit_schauer(z)

based on Schauer+21, MNRAS, Volume 507, Issue 2, pp.1775-1787 Critical mass for halo collapse with LW feedback and baryonic streaming

Parameters:

z [real,in] :: redshift

Return:

get_mcrit_schauer [real]

Called from:

get_crit_mass()

function  crit_masses/get_atomic_cooling_mass(z)

the atomic cooling mass is the virial mass of 10^4 K

Parameters:

z [real,in] :: redshift

Return:

get_atomic_cooling_mass [real]

Called from:

set_mcrit(), get_crit_mass()

Call to:

get_m_vir()

function  crit_masses/get_m_vir(z, t)

function determines the halo mass corresponding to a virial temperature we use the definition of Hummel et al 2012

Parameters:

  • z [real,in] :: redshift

  • t [real,in] :: virial temperature

Return:

get_m_vir [real]

Called from:

get_mcrit_stacy(), get_atomic_cooling_mass()

function  crit_masses/get_t_vir(m, z)
Parameters:

  • m [real,in]

  • z [real,in]

Return:

t [real]

Called from:

sf_step()

cosmic_time

Quick access

Variables:

a3_inv, alev, tlev, zlev

Routines:

cosmic_dt(), distribute_cosmic_time(), dtdz(), read_scale_file(), set_cosmic_time()

Needed modules

Variables

  • cosmic_time/a3_inv (nlev_max) [real,protected]

    1/(a^3)

  • cosmic_time/alev (nlev_max) [real,protected]

    cosmic expansion factor

  • cosmic_time/tlev (nlev_max) [real,protected]

    time since big bang

  • cosmic_time/zlev (nlev_max) [real,protected]

    redshift

Subroutines and functions

subroutine  cosmic_time/set_cosmic_time()

subroutine to set redshift steps and compute expansion factors and cosmic time

Called from:

asloth

Call to:

read_scale_file(), distribute_cosmic_time(), cosmic_dt()

subroutine  cosmic_time/distribute_cosmic_time()

subroutine to distribute redshift steps (either linearly or logarithmically)

Called from:

set_cosmic_time()

subroutine  cosmic_time/read_scale_file()

subroutine to read the redshift steps from scale file

Called from:

set_cosmic_time()

function  cosmic_time/dtdz(a)

Subroutine to calculate absolute value of dt/dz in seconds for density field normalized at reference epoch a0=1.

Notation:

Parameters:

a [real,in] ::

expansion factor rel to a0=1

Uses

Return:

res [real]

Use :

numerical_parameters (hubble_cgs(), omega_l(), omega_m())

Called from:

cosmic_dt()

function  cosmic_time/cosmic_dt(z1, z2)

determines the time in seconds between redshift z1 and z2

Parameters:

  • z1 [real,in]

  • z2 [real,in]

Return:

res [real]

Called from:

set_cosmic_time()

Call to:

dtdz()

EPS wrapper

eps_wrapper

Description

module to interface with the Parkinson et al. 2008 tree generation routines

Quick access

Routines:

make_eps_tree()

Needed modules

Subroutines and functions

subroutine  eps_wrapper/make_eps_tree()
Called from:

asloth

Call to:

set_tree_params(), memory(), make_tree(), walk_tree()

Chemistry

chemistry

Quick access

Variables:

abhe, x_ion

Routines:

ch21(), ch22(), cl4(), cl64(), fh2_crit(), h2_cooling_rate(), h2_form_rate(), ph4(), t_hot2cold()

Needed modules

Variables

  • chemistry/abhe [real,parameter=y_helium/x_hydrogen/4.0]
  • chemistry/x_ion [real,parameter=0.0002]

    residual ionization fraction

Subroutines and functions

function  chemistry/t_hot2cold(n, t, z)
Parameters:

  • n [real,in]

  • t [real,in]

  • z [real,in]

Return:

res [real]

Called from:

timescaledetermination()

Call to:

fh2_crit(), h2_form_rate()

function  chemistry/h2_form_rate(n, t, z)
Parameters:

  • n [real,in]

  • t [real,in]

  • z [real,in]

Return:

res [real]

Called from:

t_hot2cold()

Call to:

ch21(), ph4(), ch22()

function  chemistry/fh2_crit(n, t, z)
Parameters:

  • n [real,in]

  • t [real,in]

  • z [real,in]

Return:

res [real]

Called from:

t_hot2cold()

Call to:

h2_cooling_rate()

function  chemistry/h2_cooling_rate(n, t)
Parameters:

  • n [real,in]

  • t [real,in]

Return:

res [real]

Called from:

fh2_crit()

Call to:

cl4(), cl64()

function  chemistry/ch22(temp)
Parameters:

temp [real,in]

Return:

res [real]

Called from:

h2_form_rate()

function  chemistry/ch21(temp)
Parameters:

temp [real,in]

Return:

res [real]

Called from:

h2_form_rate()

function  chemistry/ph4(z)
Parameters:

z [real,in]

Return:

res [real]

Called from:

h2_form_rate()

function  chemistry/cl64(temp)
Parameters:

temp [real,in]

Return:

res [real]

Called from:

h2_cooling_rate()

function  chemistry/cl4(temp)
Parameters:

temp [real,in]

Return:

res [real]

Called from:

h2_cooling_rate()

Input files

Description

List of input files

Quick access

Variables:

lines_lt_popii, lines_mmetals_popii, lines_qion_popii, mtree_file, output_string, pkinfile_dummy, popii_lt_file, popii_mmetals_file, popii_qion_file, popii_yields_file, popiii_yields_file, scale_file, tree_file, tree_name, tree_path

Variables

  • input_files/lines_lt_popii [integer]
  • input_files/lines_mmetals_popii [integer,parameter=27]
  • input_files/lines_qion_popii [integer]
  • input_files/mtree_file [character,parameter='tree_files/mass_estimates_lx13.dat']
  • input_files/output_string [character]
  • input_files/pkinfile_dummy [character,parameter='data/kvector_z1.txt']

    is later passed on to pkinfile

  • input_files/popii_lt_file [character]
  • input_files/popii_mmetals_file [character,parameter='data/kobayashi2006_z1e-3_mmetal.txt']
  • input_files/popii_qion_file [character]
  • input_files/popii_yields_file [character,parameter='data/kobayashi2006_z1e-3.txt']
  • input_files/popiii_yields_file [character,parameter='data/yields_isotopes_sum_tr.dat']
  • input_files/scale_file [character]
  • input_files/tree_file [character]
  • input_files/tree_name [character]
  • input_files/tree_path [character]

Memory modules

Tree_memory_arrays

Description

This module is use for managing the memory of the merger tree

Quick access

Variables:

mergertree_aux

Needed modules

Variables

  • tree_memory_arrays/mergertree_aux (*) [treenode,target/allocatable]

Tree_memory_arrays_passable

Description

This module is used to reference the merger tree

Quick access

Variables:

mergertree, number_of_nodes

Routines:

tree_index()

Needed modules

Variables

  • tree_memory_arrays_passable/mergertree (*) [treenode,pointer]
  • tree_memory_arrays_passable/number_of_nodes [integer]

Subroutines and functions

function  tree_memory_arrays_passable/tree_index(node)
Parameters:

node [treenode,pointer]

Return:

tree_index [integer,pure]

Called from:

associate_siblings()

Star formation

SF_routines

Quick access

Routines:

add_hmstars(), add_to_base_mstariisurv(), coldgasbindingenergy_cold(), coldgasbindingenergy_hot(), coldgasbindingenergy_nfw(), coldgasbindingenergy_stellar(), dm_concen(), get_n_ion(), hms_feedback(), hotgasbindingenergy_disk(), hotgasbindingenergy_hot(), hotgasbindingenergy_nfw(), imfsampling(), nfw_const(), rm_hmstars(), snedrivenoutflow(), survivingstarlist(), timescaledetermination(), writellstarstofileanddel()

Needed modules

Subroutines and functions

subroutine  sf_routines/add_hmstars(p_node, ii, timeborn, num_, zgas, pop)
Parameters:

  • p_node [treenode,inout]

  • ii [integer,in]

  • timeborn [real,in]

  • num [integer,in]

  • zgas (n_elements) [real,in]

  • pop [integer,in]

Called from:

imfsampling()

subroutine  sf_routines/rm_hmstars(p_node, highstar)
Parameters:
Called from:

hms_feedback()

function  sf_routines/get_n_ion(p_node, t_current, del_t)

This routine calculates the total number of ionizing photons emitted during this time-step

Parameters:

  • p_node [treenode,inout] :: the halo

  • t_current [real,in] :: beginning of the sub-step

  • del_t [real,in] :: time-step

Return:

n_ion_tot [real] :: number of total ionizing photons

Called from:

sf_step()

subroutine  sf_routines/hms_feedback(p_node, t_current, totalmetalejectathisstep, totalsneenergythisstep, m_heat, rel_del_t)

This routine calculates the cold gas mass that will be heated/ blown out from the halo and the explosion energy/ ejecta from SNe.

Parameters:

  • p_node [treenode,inout] :: the halo

  • t_current [real,in] :: beginnng of the sub-step

  • totalmetalejectathisstep (n_elements) [real,inout] :: metals produced in this time-step

  • totalsneenergythisstep [real,out] :: Msun cm^2 s^-2

  • m_heat [real,out] :: mass that will trasfer from cold to hot in this time-step

  • rel_del_t [real,in] :: time-step

Called from:

sf_step()

Call to:

add_to_m_metals(), rm_hmstars()

subroutine  sf_routines/add_to_base_mstariisurv(this_node)

Inherit M_star_II_surv to the base node at the final redshift

Parameters:

this_node [treenode,inout]

Called from:

asloth

subroutine  sf_routines/survivingstarlist(p_node, num, ii, z, timeborn, pop)

Non-standard routine if users are interested in stellar information of survival stars in satellites

Parameters:

  • p_node [treenode,inout] :: the halo

  • num [integer,in] :: number of stars in this IMF mass bin

  • ii [integer,in] :: index of the IMF mass bin

  • z (n_elements) [real,in] :: metallicity of the star

  • timeborn [real,in] :: time of star formation

  • pop [integer,in] :: stellar population

Called from:

imfsampling()

subroutine  sf_routines/writellstarstofileanddel(p_node)
Parameters:

p_node [treenode,inout]

Called from:

sf_step()

subroutine  sf_routines/imfsampling(p_node, totalmassnewstars, m_star_form, t_current, this_rng, zgas, pop)

This routine samples individual stars from the IMF and store massive stars in the list.

Parameters:

  • p_node [treenode,inout] :: the halo

  • totalmassnewstars [real,inout] :: total mass of newly formed stars

  • m_star_form [real,in] :: estimated total stellar mass based on cold gas mass and star formation efficiency

  • t_current [real,in] :: the time in the beginning of the time-step

  • this_rng [rng,inout] :: random number generator

  • zgas (n_elements) [real,in] :: gas metallicity

  • pop [population,in] :: stellar population

Use :

trace_mdf

Called from:

sf_step()

Call to:

get_mdf_bin(), poisson(), add_hmstars(), survivingstarlist(), add_to_node_mdf()

subroutine  sf_routines/dm_concen(m_halo, z, dm_concentration)

This routine determins the concentration of the NFW profile. We adopt the fitting curve by Correa et al. 2015 (doi:10.1093/mnras/stv1363) Planck cosmology (appendix B) Note that the fitting functions are valid for z=0-10 based on private conversation

Parameters:

  • m_halo [real,in] :: virial mass of the halo

  • z [real,in] :: current redshift

  • dm_concentration [real,out] :: (arbitrary) constant value, reference from Fig. 7 in Correa et al. 2015 for haloes > 10^5 Msun

Called from:

sf_step()

subroutine  sf_routines/timescaledetermination(p_node, v_dyn, t_dyn, ncold, n_hot, coldgas_tff, r_vir_current, dm_concentration, t_vir, t_cool)

This routines determines all the relative timescales that we need. Note that the time scales we eventually use are dimensionless.

Parameters:

  • p_node [treenode,inout] :: This also includes “dead” stars

  • v_dyn [real,inout] :: s; dynamical time of the halo center

  • t_dyn [real,inout] :: dynamical timescale of central part of the halo

  • ncold [real,inout] :: number density of cold gas

  • n_hot [real,inout] :: number density of hot gas

  • coldgas_tff [real,inout] :: free fall time of cold gas

  • r_vir_current [real,in] :: the current virial raius of the halo

  • dm_concentration [real,in] :: dark matter halo concentration

  • t_vir [real,in] :: virial temperature of the halo

  • t_cool [real,inout] :: cooling time of hot gas

Called from:

sf_step()

Call to:

t_hot2cold()

function  sf_routines/nfw_const(r_s, r_vir_current)
Parameters:

  • r_s [real,in]

  • r_vir_current [real,in]

Return:

a [real]

Called from:

hotgasbindingenergy_nfw(), coldgasbindingenergy_nfw()

function  sf_routines/hotgasbindingenergy_nfw(r_s, r_vir_current, m_hot, m_peak)
Parameters:

  • r_s [real,in]

  • r_vir_current [real,in]

  • m_hot [real,in]

  • m_peak [real,in]

Return:

be_hot_nfw [real]

Called from:

snedrivenoutflow()

Call to:

nfw_const()

function  sf_routines/hotgasbindingenergy_disk(r_s, r_vir_current, m_hot, m_disk)
Parameters:

  • r_s [real,in]

  • r_vir_current [real,in]

  • m_hot [real,in]

  • m_disk [real,in]

Return:

be_hot_disk [real]

Called from:

snedrivenoutflow()

function  sf_routines/hotgasbindingenergy_hot(r_vir_current, m_hot)
Parameters:

  • r_vir_current [real,in]

  • m_hot [real,in]

Return:

be_hot_hot [real]

Called from:

snedrivenoutflow()

function  sf_routines/coldgasbindingenergy_nfw(r_s, r_vir_current, m_cold, m_peak)
Parameters:

  • r_s [real,in]

  • r_vir_current [real,in]

  • m_cold [real,in]

  • m_peak [real,in]

Return:

be_cold_nfw [real]

Called from:

snedrivenoutflow()

Call to:

nfw_const()

function  sf_routines/coldgasbindingenergy_stellar(r_s, m_cold, m_stellar)
Parameters:

  • r_s [real,in]

  • m_cold [real,in]

  • m_stellar [real,in]

Return:

be_cold_stellar [real]

Called from:

snedrivenoutflow()

function  sf_routines/coldgasbindingenergy_hot(r_s, r_vir_current, m_cold, m_hot)
Parameters:

  • r_s [real,in]

  • r_vir_current [real,in]

  • m_cold [real,in]

  • m_hot [real,in]

Return:

be_cold_hot [real]

Called from:

snedrivenoutflow()

function  sf_routines/coldgasbindingenergy_cold(r_s, m_cold)
Parameters:

  • r_s [real,in]

  • m_cold [real,in]

Return:

be_cold_cold [real]

Called from:

snedrivenoutflow()

subroutine  sf_routines/snedrivenoutflow(p_node, r_vir_current, dm_concentration, totalsneenergythisstep, del_m_hot, totalmassnewstars, m_out_hot, m_out_cold)

This routines determines the binding energy of hot gas and cold gas. It then compare the SNe energy and the binding energy to determine how much gas SNe will blow out of the halo.

Parameters:

  • p_node [treenode,inout] :: the halo

  • r_vir_current [real,in] :: current virial radius of the halo

  • dm_concentration [real,in] :: dark matter halo’s concentration

  • totalsneenergythisstep [real,inout] :: sum of SNe energy that occur in this time-step

  • del_m_hot [real,in] :: smoothly accreted hot gas in this time-step

  • totalmassnewstars [real,in] :: total newly formed stellar mass in this time-step

  • m_out_hot [real,inout] :: hot gas mass that ends up in outflow in this time-step

  • m_out_cold [real,inout] :: cold gas mass that ends up in outflow in this time-step

Called from:

sf_step()

Call to:

hotgasbindingenergy_nfw(), hotgasbindingenergy_disk(), hotgasbindingenergy_hot(), coldgasbindingenergy_nfw(), coldgasbindingenergy_stellar(), coldgasbindingenergy_hot(), coldgasbindingenergy_cold()

Star_formation

Quick access

Routines:

add_parent_to_tree(), dt_adaptive(), m_star_dot(), prepare_node(), sf_step()

Needed modules

Subroutines and functions

subroutine  star_formation/prepare_node(p_node, m_star_max, is_ion)
Parameters:

  • p_node [treenode,inout] :: parent node

  • m_star_max [real,out]

  • is_ion [logical,out] :: set from most massive child

Called from:

sf_step()

Call to:

check_m_metals(), inherit_igm_z()

subroutine  star_formation/sf_step(node)
Parameters:

node [treenode,inout,target]

Called from:

asloth

Call to:

prepare_node(), rvir_mpc(), dm_concen(), get_t_vir(), check_m_metals(), get_metallicity(), calc_dz(), metalpoor(), timescaledetermination(), dt_adaptive(), m_star_dot(), imfsampling(), get_n_ion(), hms_feedback(), snedrivenoutflow(), assert(), v_plow_ad(), r_ion_ad(), writellstarstofileanddel()

subroutine  star_formation/add_parent_to_tree(node)
Parameters:

node [treenode,inout]

Called from:

asloth

Call to:

add_bubble_to_tree(), get_m_metals()

function  star_formation/dt_adaptive(p_node, coldgas_tff, t_dyn, del_m_hot, del_t, t_step, pop_now)

timestep in sec

Parameters:

  • p_node [treenode,in]

  • coldgas_tff [real,in]

  • t_dyn [real,in]

  • del_m_hot [real,in]

  • del_t [real,in]

  • t_step [real,in]

  • pop_now [integer,in]

Return:

dt_adaptive [real]

Called from:

sf_step()

Call to:

m_star_dot()

function  star_formation/m_star_dot(m_cold, coldgas_tff, pop_now)

star formation rate in this substep (Msun/s)

Parameters:

  • m_cold [real,in]

  • coldgas_tff [real,in]

  • pop_now [integer,in]

Return:

m_star_dot [real]

Use :

numerical_parameters (etaii(), etaiii())

Called from:

sf_step(), dt_adaptive()

Stellar feedback

feedback_routines

Quick access

Routines:

node_feedback(), self_enrichment(), statistical_feedback()

Needed modules

Subroutines and functions

subroutine  feedback_routines/node_feedback(this_node)
Parameters:

this_node [treenode,target]

Called from:

asloth

Call to:

self_enrichment(), statistical_feedback(), bubble_check(), ion_check()

subroutine  feedback_routines/self_enrichment(node)

This routine computes feedback exactly as the bubble-tree in the case that there is only self-enrichment

Parameters:

node [treenode]

Called from:

node_feedback()

subroutine  feedback_routines/statistical_feedback(this_node)
Parameters:

this_node [treenode,target]

Use :

random_object

Called from:

node_feedback()

Call to:

get_ran()

snowplow

Quick access

Routines:

r_ion_ad(), v_plow_ad()

Needed modules

Subroutines and functions

function  snowplow/v_plow_ad(m, z, r, r_v)
Parameters:

  • m [real,in]

  • z [real,in]

  • r [real,in]

  • r_v [real,in]

Return:

v [real]

Called from:

sf_step()

function  snowplow/r_ion_ad(nion, r_old, z, del_t, mhalo)
Parameters:

  • nion [real,in]

  • r_old [real,in]

  • z [real,in]

  • del_t [real,in]

  • mhalo [real,in]

Return:

r_new [real]

Called from:

sf_step()

Call to:

rvir_mpc(), rvir()

Feedback tracing

feedback_arrays

Quick access

Variables:

sfr, v_enriched, v_ionized

Routines:

add_to_feedback_volumes(), add_to_sfr()

Needed modules

Variables

  • feedback_arrays/sfr (2,nlev_max) [real,protected/optional/default=0]

    SFR: population, timestep

  • feedback_arrays/v_enriched (nlev_max) [real,protected/optional/default=0]
  • feedback_arrays/v_ionized (nlev_max) [real,protected/optional/default=0]

    (in physical Mpc^3)

Subroutines and functions

subroutine  feedback_arrays/add_to_feedback_volumes(node)
Parameters:

node [treenode]

Called from:

asloth

subroutine  feedback_arrays/add_to_sfr(node)
Parameters:

node [treenode]

Called from:

asloth

volume_fractions

Quick access

Routines:

get_fracs()

Needed modules

Subroutines and functions

subroutine  volume_fractions/get_fracs(jlevel)
Parameters:

jlevel [integer,in]

Use :

random_object

Called from:

asloth

Call to:

get_ran(), position_check()

Outputs

to_file

Quick access

Types:

file_specs

Routines:

add_timestamp(), close_outputs(), init_outputs(), outputs_satellites(), write_file(), write_files(), write_files_reali(), write_mw_properties()

Needed modules

Types

  • type  to_file/file_specs
    Type fields:

    • % file_string [character]

    • % head1 [character]

    • % head2 [character]

    • % out_dir [character]

    • % stuff (30,nlev_max) [real]

    • % var_name [character]

Variables

Subroutines and functions

subroutine  to_file/init_outputs()
Called from:

asloth

Call to:

add_timestamp()

subroutine  to_file/write_file(fs, var_name, ncoll)
Parameters:

  • fs [file_specs]

  • var_name [character]

  • ncoll [integer] :: loop writes one line

Called from:

write_files()

subroutine  to_file/outputs_satellites(count)
Parameters:

count [integer]

Use :

tree_memory_arrays_passable, trace_mdf

Called from:

asloth

Call to:

rvir_mpc(), add_to_node_mdf()

subroutine  to_file/write_files_reali()
Use :

tree_memory_arrays_passable, input_files, metal_functions, trace_mdf

Called from:

asloth

Call to:

rvir_mpc(), add_to_node_mdf()

subroutine  to_file/write_files(output)
Parameters:

output [file_specs,inout]

Called from:

asloth

Call to:

real_to_string(), write_file()

subroutine  to_file/close_outputs()

go through list of possible open file IDs

Called from:

asloth

subroutine  to_file/add_timestamp(output_path)
Parameters:

output_path [character,inout]

Called from:

init_outputs()

Call to:

int_to_string4(), int_to_string2()

subroutine  to_file/write_mw_properties()
Use :

tree_memory_arrays_passable, input_files

Called from:

asloth