The different engines behave differently when it comes to evolving porosity and surface area.
Check whether it's implemented for PFLOTRAN. Implement.
Check whether it works correctly for CrunchFlow, add test examples.
Does Alquimia need to provide information to driver that porosity and surface area is being evolved?