Christopher, you should use implicitFreeSurface=.TRUE., This will considerably reduce solver iteration for a given accuracy. After first few time steps, iterations should not exceed order 100 for 1e-7 residual. Dimitris