modify skirtor module so that the central source is isotropic

pcigale/sed_modules/fritz2006.py

@@ -125,9 +125,7 @@ class Fritz2006(SedModule):
                                   x=self.fritz2006.wave)
 
         # Apply polar-dust obscuration
-        # We define various constants necessary to compute the model. For
-        # consistency, we define the speed of light in nm s¯¹ rather than in
-        # m s¯¹.
+        # We define various constants necessary to compute the model
         self.c = cst.c * 1e9
         lambda_0 = 200e3
         # Calculate the extinction (SMC)

pcigale/sed_modules/skirtor2016.py

@@ -101,6 +101,15 @@ class SKIRTOR2016(SedModule):
                                                     self.oa, self.R, self.Mcl,
                                                     self.i)
 
+        # Modify the model so that the central source is isotropic
+        # Calculate the coverting factors disk emission
+        # (only for type-1 case, since disk emission in type-2 is zero)
+        # (since dust emission is normalized to unity; it remains same)
+        if self.i<=(90-self.oa):
+            self.SKIRTOR2016.disk*= (4*np.sin(self.oa*np.pi/180)**2 + 3*np.sin(self.oa*np.pi/180)) / \
+                                    (12*np.cos(self.i*np.pi/180)**2 + 6*np.cos(self.i *np.pi/180))
+
+
     def process(self, sed):
         """Add the IR re-emission contributions