enable to calculate AGN accretion power

pcigale/sed_modules/fritz2006.py

@@ -234,7 +234,7 @@ class Fritz2006(SedModule):
         sed.add_info('agn.scatt_luminosity', l_agn_scatt, True)
         sed.add_info('agn.disk_luminosity', l_agn_agn, True)
         sed.add_info('agn.luminosity', l_agn_total, True)
-        sed.add_info('agn.intrin_disk_luminosity', l_agn_intrin_agn, True)
+        sed.add_info('agn.accretion_power', l_agn_intrin_agn, True)
         sed.add_info('agn.intrin_Lnu_2500A', l_agn_2500A, True)
 
         sed.add_contribution('agn.fritz2006_therm', self.fritz2006.wave,

pcigale/sed_modules/skirtor2016.py

@@ -161,10 +161,10 @@ class SKIRTOR2016(SedModule):
 
         # Integrate AGN luminosity for different components
         self.lumin_disk = np.trapz(self.SKIRTOR2016.disk, x=self.SKIRTOR2016.wave)
-        # Intrinsic (de-reddened) AGN luminosity from the central source
+        # Intrinsic (de-reddened) AGN luminosity from the central source at theta=30 deg
         self.lumin_intrin_disk = np.trapz(self.SKIRTOR2016.intrin_disk,
                                           x=self.SKIRTOR2016.wave)
-        # Calculate L_lam(2500A)
+        # Calculate L_lam(2500A) at theta=30 deg
         self.l_agn_2500A = np.interp(250, self.SKIRTOR2016.wave, self.SKIRTOR2016.intrin_disk)
         # Convert L_lam to L_nu
         self.l_agn_2500A *= 250**2/self.c
@@ -201,7 +201,11 @@ class SKIRTOR2016(SedModule):
             agn_power = luminosity * (1./(1.-self.fracAGN) - 1.)
             lumin_dust = agn_power
             lumin_disk = agn_power * self.lumin_disk
-            lumin_intrin_disk = agn_power * self.lumin_intrin_disk
+            # power_accretion means the intrinsic disk luminosity
+            # integrated over 4pi solid angles
+            # The factor (0.493) comes from the fact that lumin_intrin_disk
+            # is calculated at viewing angle = 30 deg
+            power_accretion = agn_power * self.lumin_intrin_disk * 0.493
             lumin = lumin_dust + lumin_disk
             l_agn_2500A = agn_power * self.l_agn_2500A
         else:
@@ -211,7 +215,7 @@ class SKIRTOR2016(SedModule):
         sed.add_info('agn.dust_luminosity', lumin_dust, True)
         sed.add_info('agn.disk_luminosity', lumin_disk, True)
         sed.add_info('agn.luminosity', lumin, True)
-        sed.add_info('agn.intrin_disk_luminosity', lumin_intrin_disk, True)
+        sed.add_info('agn.accretion_power', power_accretion, True)
         sed.add_info('agn.intrin_Lnu_2500A', l_agn_2500A, True)
 
         sed.add_contribution('agn.SKIRTOR2016_dust', self.SKIRTOR2016.wave,