Add to paradox_sim docs; remove comments; return burned ts now
Re document
Add crude axes to filled_contour()
Fix paradox_sim_R to match C++ version
A couple minor changes to match new function names
Add paradox_sim_reps which runs sim over multiple reps easily

And in parallel if desired
Add function to setup parallel
Add early version of crude diagnostic plot
Add some diagnostic options to paradox_sim.cpp
Re-document