Abstract
The Laplace transform technique applied as an initialisation and a time-integration scheme of a baroclinic primitive equations model is outlined. The numerical stability of this procedure and the selective damping of the high-frequency waves are improved by the introduction of an elliptic domain in the complex projection field. The time-step is doubled compared to a split-explicit integration scheme.