Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 5 additions & 3 deletions src/disk_instability.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -188,9 +188,11 @@ void DiskInstability::create_starburst(SubhaloPtr &subhalo, GalaxyPtr &galaxy, d
// Define accretion rate.
galaxy->smbh.macc_sb = delta_mbh/tdyn;

// Grow SMBH.
galaxy->smbh.mass += delta_mbh;
galaxy->smbh.mass_metals += delta_mzbh;
// Check if SMBH exists in this galaxy, then grow it.
if(galaxy->smbh.mass != 0){
galaxy->smbh.mass += delta_mbh;
galaxy->smbh.mass_metals += delta_mzbh;
}

// Reduce gas available for star formation due to black hole growth.
galaxy->bulge_gas.mass -= delta_mbh;
Expand Down
7 changes: 5 additions & 2 deletions src/galaxy_mergers.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -540,8 +540,11 @@ void GalaxyMergers::create_starbursts(HaloPtr &halo, double z, double delta_t){
physicalmodel->evolve_galaxy_starburst(*subhalo, *galaxy, z, delta_t, true);

// Grow SMBH after starbursts, as during it we need a realistical measurement of Ledd the BH had before the starburst.
galaxy->smbh.mass += delta_mbh;
galaxy->smbh.mass_metals += delta_mzbh;
// Check if SMBH exists in this galaxy, then grow it.
if(galaxy->smbh.mass != 0){
galaxy->smbh.mass += delta_mbh;
galaxy->smbh.mass_metals += delta_mzbh;
}

// Check for small gas reservoirs left in the bulge, in case mass is small, transfer to disk.
if(galaxy->bulge_gas.mass > 0 && galaxy->bulge_gas.mass < parameters.mass_min){
Expand Down