We complete the computation of the gravitational waveform amplitude for
non-spinning compact binaries to the third and a half post-Newtonian (PN) order
in the quasi-circular case. This is done by performing a spin-weighted
spherical harmonics decomposition of the amplitude. This computation is
achieved using the post-Newtonian-Multipolar-post-Minkowskian formalism. The
amplitude modes are written in a PN expanded way as well as a factorized form
suitable for EOB template building. Combining this work with previous ones, we
provide in the supplemental materials the expressions of the modes including
all non-spinning and spinning effects to the 3.5PN order for non-precessing
quasi-circular binaries.