Modelling of coupled normal modes of the Earth: the spectral method