tried to implement the cosmology

choosing inittype = 4, leads to NAN. not sure why, but also not sure
what whether the stuff I did makes makes much sense…