  We study the two-field warm inflation models with a double quadratic
potential and a linear temperature dependent dissipative coefficient. We
derived the evolution equation of all kinds of perturbations without assuming
slow-roll approximation, and obtained the curvature power spectrum at the end
of inflation with a fully numerical method. Then we compute the scalar spectral
index $n_s$, tensor-to-scalar ratio $r$ for several representative potentials,
and compare our results with observational data. At last, we use Planck data to
constrain the parameters in our models. This work is a natural extension of
single-field warm inflation, and the aim of this work is to present some
features of multi-field warm inflation using a simple two-field model.
