change frac  = ionised fraction