Soil hydrology constants are soil type dependent - create functions to calculate across grid

We noticed that some constants are shared across soil and hydrology model and we stored them in different places, this will be addressed with #768

In the longer term, we want these constants NOT to be floats but rather soil type dependent arrays. The hydrology model is already set up to take arrays, for example for soil moisture capacity.

We would like a function to calculate these constants based on input soil maps. This could be part of the `init` of the model, need to discuss where to best place it (core or hydrology or soil).

Start by reviewing the shared-constant work in #768 and the model's init path. Compare the core, hydrology, and soil areas to determine where soil-map inputs and soil-type-dependent arrays belong. Done means the constants are calculated from soil maps and supplied as arrays to the hydrology model, with the integration location agreed.

Written by the indexing model from the issue text.