Added correct density profile, from Piro and Nakar 2014, can now re-create Appendix B plots from Piro 2013. Only with small t_peak