are computed directly from the Lagrangian interface and then distributed onto the Eulerian grid using the discrete delta function and the immersed boundary method of [@Pes1977]. A detailed description of the procedure for calculating $\textbf{F}$, $\textbf{N}$ and $I_{\rm num}$ can be found in [@ShJu07; @ShJu09], where in particular we demonstrate that this method of calculating the surface tension force reduces any parasitic currents in the standard static drop test case to a level of $O(10^{-7})$ for fluids with properties similar to those we use here.

The Lagrangian interface is advected by integrating $\mathrm{d}\textbf{x}_f/\mathrm{d}t = \textbf{V}$ with a second-order Runge-Kutta method where the interface velocity, $\textbf{V}$, is interpolated from the Eulerian velocity.

The parallelization of the code is based on algebraic domain decomposition, where the velocity field is solved by a parallel generalized minimum residual (GMRES) method for the implicit viscous terms and the pressure by a parallel multigrid method motivated by the algorithm of [@Kwak-InterScience-2004]. Communication across process threads is handled by message passing interface (MPI) procedures.

The code contains a module for the definition of immersed solid objects and their interaction with the flow, which we have used to simulate Faraday waves in a spherical container. In order to 