Performance: avoid unnecessary heap allocations

[SouthEndMusic]

I have found, for instance looking at the sediment code, that a lot of computations are done by allocating arrays in update! instead of using pre-allocated caches. I haven't tested the performance myself, but it's good practice to avoid heap allocations in simulation code. This can be investigated rigorously using @profview allocs (see the documentation here.