Move the content of the loop over scales in a separate function, so it can be called with fresh arguments every time to check parallelism and dependencies.

- turn C4 from a structure to a class (to do the above), move its initialisation to its constructor
- eliminate "allocate_vectors()", redundant and including a (wrong) double allocation