Add a final checkpoint

src/analysis.py

@@ -557,6 +557,11 @@ with MPIPool() as pool:
     for it_final, res in enumerate(sampler.add_live_points()):
         pass
 
+    # Create a final checkpoint and set of plots
+    write_current_state(sampler, resume_file, sampling_time)
+    if input_args.no_plot is False:
+        plot_current_state(sampler, sampling_keys, outdir, label)
+
     sampling_time += (datetime.datetime.now() - t0).total_seconds()
 
     out = sampler.results