If a pdep network causes problems, improve logging for debugging.

If a pdep network is going to cause the job to crash with an
InvalidMicrocanonicalRateError, then report a summary of 
the troublesome network and draw a PDF of it before
crashing, to make it easier to figure out what's wrong.

Author: rwest

rmgpy/pdep/network.py

@@ -408,6 +408,7 @@ def set_conditions(self, T, P, ymB=None):
                             logging.error("Increasing number of grains did not decrease error enough "
                                           "(Current badness: {0:.1f}, previous {1:.1f}). Something must be wrong with "
                                           "network {2}".format(badness, previous_error.badness(), self.label))
+                            self.log_summary()
                             raise error
                     previous_error = error
                     success = False

rmgpy/rmg/pdep.py

@@ -44,7 +44,7 @@
 import rmgpy.reaction
 from rmgpy.constants import R
 from rmgpy.data.kinetics.library import LibraryReaction
-from rmgpy.exceptions import PressureDependenceError, NetworkError
+from rmgpy.exceptions import PressureDependenceError, NetworkError, InvalidMicrocanonicalRateError
 from rmgpy.pdep import Configuration
 from rmgpy.rmg.react import react_species
 from rmgpy.statmech import Conformer
@@ -876,7 +876,13 @@ def update(self, reaction_model, pdep_settings, requires_rms=False):
         # Calculate the rate coefficients
         self.initialize(Tmin, Tmax, Pmin, Pmax, maximum_grain_size, minimum_grain_count, active_j_rotor, active_k_rotor,
                         rmgmode)
-        K = self.calculate_rate_coefficients(Tlist, Plist, method)
+        try:
+            K = self.calculate_rate_coefficients(Tlist, Plist, method)
+        except InvalidMicrocanonicalRateError:
+            if output_directory:
+                job.draw(output_directory, file_format='pdf')
+                logging.info(f"Network {self.index} has been drawn and saved as a pdf in {output_directory}.")
+            raise
 
         # Generate PDepReaction objects
         configurations = []