change AGN disk SED shape, to avoid UV overestimation

pcigale/sed_modules/skirtor2016.py

@@ -120,6 +120,29 @@ class SKIRTOR2016(SedModule):
                                                     self.oa, self.R, self.Mcl,
                                                     self.i)
 
+        # Modify the disk SED to Feltre et al. (2012), following the method in
+        # https://sites.google.com/site/skirtorus/sed-library
+        # Set the broken power law shape
+        # 1st part: < 50 nm
+        disk_new = self.SKIRTOR2016.wave.copy()
+        # 2nd part: 50-125 nm
+        wave_idxs = (self.SKIRTOR2016.wave>=50) & (self.SKIRTOR2016.wave<125)
+        # Here, the factor 50**0.2 is to re-nomalize for continuity
+        norm_fac = 50**0.2
+        disk_new[wave_idxs] = disk_new[wave_idxs]**0.8 * norm_fac
+        # 3rd part: 125-1e4 nm
+        norm_fac *= 125**2.3
+        wave_idxs = (self.SKIRTOR2016.wave>=125) & (self.SKIRTOR2016.wave<1e4)
+        disk_new[wave_idxs] = disk_new[wave_idxs]**-1.5 * norm_fac
+        # 4th part: >1e4 nm
+        norm_fac *= 1e4**2.5
+        wave_idxs = self.SKIRTOR2016.wave>=1e4
+        disk_new[wave_idxs] = disk_new[wave_idxs]**-4 * norm_fac
+        # Re-normalize the entire disk SED to keep energy conservation
+        disk_new *= np.trapz(self.SKIRTOR2016.disk, x=self.SKIRTOR2016.wave) / \
+                    np.trapz(disk_new, x=self.SKIRTOR2016.wave)
+        self.SKIRTOR2016.disk = disk_new
+
         # Apply wavelength cut to avoid X-ray wavelength
         lam_cut = 5.
         lam_idxs = self.SKIRTOR2016.wave>=lam_cut