Incorporating Matrix ionization scheme 

This pull request would incorporate the matrix scheme into dev.

Currently in beta, but @kslong wanted it merged in. If you still do Knox, then merge this pull request and consult my recent email. 

The recent reorganisation I did to the matrix scheme is described here https://github.com/agnwinds/python/pull/130