Added additional debug output for zero time steps. More attempts at symmetrising the HLLC Riemann solver. Removed output message for every vacuum Riemann solver call in Exact Riemann solver. Made temperature and density axis logarithmic in Zeldovich plot script and removed analytic solution for redshifts below 1.