Changed eq_state.f90 to speed up computation by doing Saha only
in species with significant aboundances (gt 1e-6 times H)
Going from 1e-4 to 1-7 doesn't change results or time significantly