TODO Find a suitable way to obtain the initial guess Handle solid and species concentration in a single array Reduce problem for species distribution (beta to beta prime conversion)