  In this work we present IMRPhenomTP, a time domain phenomenological model for
the dominant $l=2$, $m=|2|$ modes of coalescing black hole binary systems and
its extension to describe general precessing systems within the "twisting up"
approximation. The underlying non-precessing model is calibrated to the new
release of Numerical Relativity simulations of the SXS Collaboration and its
accuracy is comparable to the state-of-the-art non-precessing dominant mode
models as IMRPhenomX and SEOBNRv4. The precessing extension allows for
flexibility choosing the Euler angles of the time-dependent rotation between
the co-precessing and the inertial reference systems, including the single spin
NNLO and the double spin MSA PN descriptions present in other models, numerical
integration of the orbit averaged spin evolution equations, different choices
for the evolution of the orbital angular momentum norm and a simple
approximation to the ringdown behaviour.
