print IIS if solver returns status infeasible

scripts/solve_network.py

@@ -839,6 +839,9 @@ def solve_network(n, config, solving, opts="", **kwargs):
             f"Solving status '{status}' with termination condition '{condition}'"
         )
     if "infeasible" in condition:
+        labels = n.model.compute_infeasibilities()
+        logger.info("Labels:\n" + labels)
+        n.model.print_infeasibilities()
         raise RuntimeError("Solving status 'infeasible'")
 
     return n