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.