Only printout and save every 5000 step to file. Choose which paramater to give by equation of state. Close to final version of program. Only simulates cirka 9% of mass before P_gas < 0.