Add neutral H condition when iterating on ionization state

Proposed fix for #89