[PL/CT] Store chemical system id at the gauss integration point
As suggested by @endJunction,
chemical_system_id is moved and stored at the gauss integration point. This change brings added convenience to chemical solver interface.
- Feature description was added to the changelog
- Tests covering your feature were added?
- Any new feature or behavior change was documented?