Array of species  moved to Reaction record.

- move record for array of species to abstract reaction class
- mapping from local indices of substances to their global indices done in abstract reaction class
- clear a bit the reaction and sorption classes
Dual porostity model works, further cleaning of reaction interface done.

- test 23 added, old and new dual porosity compared successfully
- TOS - simplified reaction constructing but further development required
- reaction interface simplified again
- other small changes in sorptions - mainly simplifiing and cleaning
Reaction interface edited, new sorptions enabled.

- methods for setting the sorbed concentration hidden from reaction interface
- methods modify_reaction_matrix hidden from linear reaction interface
- test20 - new sorption under reaction interface enabled and tested
Improve reaction interface and change sorption record name.

- reaction interfaced cleaned
- reaction: ini_from_input hidden in protected
- sorption record name changed to 'Sorption'