I am only adding the functions for the drift/kick/timestep and not the structure around them. I am also implementing a new set of min/max functions because we will need them when computing the time step of a cell (min/max for gparts, parts, sparts, bparts and sinks -> 5 values).