Inflationary models based on a single scalar field $\phi$ with a quadratic potential $V = \frac{1}{2} m^2 \phi^2$ are disfavoured by the recent Planck constraints on the scalar index, $n_s$, and the tensor-to-scalar ratio for cosmological density perturbations, $r_T$. In this paper we study how such a quadratic inflationary model can be rescued by postulating additional fields with quadratic potentials, such as might occur in sneutrino models, which might serve as either curvatons or supplementary inflatons. Introducing a second scalar field reduces but does not remove the pressure on quadratic inflation, but we find a sample of three-field models that are highly compatible with the Planck data on $n_s$ and $r_T$. We exhibit a specific three-sneutrino example that is also compatible with the data on neutrino mass difference and mixing angles.

We analyse the situation where the primordial curvature perturbations are produced by the joint effects of an inflaton field and two curvaton fields. We present general equations which allow the reader to obtain f_NL for several different scenarios which differ in the order in which fields decay into radiation after inflation. In order to investigate the physics of these equations we analyse some simplified situations where the fields are harmonic and both curvatons are frozen at the same expectation value during inflation. We find quite complex behaviour - for a given situation where the inflaton contributes a fixed amount to the total curvature perturbation there are situations where f_NL is maximised if both curvatons share equally in contributing the rest and situations where f_NL is maximised if only one of the two curvatons contributes the rest. We are unable therefore to make any completely general extrapolations about the expected non-gaussianity from N curvatons. We find that as the curvaton contribution to the overall perturbation is gradually increased, f_NL rises to a maximum before falling again and that a given f_NL can correspond to many different parameter sets for the two curvatons.