Merge pull request #201 from mklauser/fix189

A critical warning was added to the model indicating that no packet has left the simulation through the outer boundary.

tardis/model.py

@@ -359,6 +359,9 @@ def simulate(self, update_radiation_field=True, enable_virtual=False, initialize
             montecarlo.montecarlo_radial1d(self,
                                                      virtual_packet_flag=no_of_virtual_packets)
 
+        if np.sum(montecarlo_energies < 0) == len(montecarlo_energies):
+            logger.critical("No r-packet escaped through the outer boundary.")
+
         self.montecarlo_nu = montecarlo_nu * u.Hz
         self.montecarlo_luminosity = montecarlo_energies *  1 * u.erg / self.time_of_simulation