Replace ioneq with ionization_fraction and add setter

Replace ioneq with ionization_fraction and add setter #426