Experiences (continued)
The code generally vectorizes well, because much of it involves loops over all gridpoints (IJK) and loops over all droplets (N).
However, the coupling terms between flowfield and droplets vectorize only moderately because the droplets move relative to the grid and are not coincident with the grid points.
We chose an algorithm wherein the droplets are tied to the grid that they are physically located in, with the intention of minimizing communication in the parallel version.