We propose a statistical model for homogeneous turbulence undergoing distortions, which improves and extends the MCS model by Mons, Cambon & Sagaut (J. Fluid Mech., vol. 788, 2016, 147–182). The spectral tensor of two-point second-order velocity correlations is predicted in the presence of arbitrary mean-velocity gradients and in a rotating frame. For this, we numerically solve coupled equations for the angle-dependent energy spectrum ${\mathcal{E}}(\boldsymbol{k},t)$ that includes directional anisotropy, and for the deviatoric pseudo-scalar $Z(\boldsymbol{k},t)$, that underlies polarization anisotropy ($\boldsymbol{k}$ is the wavevector, $t$ the time). These equations include two parts: (i) exact linear terms representing the viscous spectral linear theory (SLT) when considered alone; (ii) generalized transfer terms mediated by two-point third-order correlations. In contrast with MCS, our model retains the complete angular dependence of the linear terms, whereas the nonlinear transfer terms are closed by a reduced anisotropic eddy damped quasi-normal Markovian (EDQNM) technique similar to MCS, based on truncated angular harmonics expansions. And in contrast with most spectral approaches based on characteristic methods to represent mean-velocity gradient terms, we use high-order finite-difference schemes (FDSs). The resulting model is applied to homogeneous rotating turbulent shear flow with several Coriolis parameters and constant mean shear rate. First, we assess the validity of the model in the linear limit. We observe satisfactory agreement with existing numerical SLT results and with theoretical results for flows without rotation. Second, fully nonlinear results are obtained, which compare well to existing direct numerical simulation (DNS) results. In both regimes, the new model improves significantly the MCS model predictions. However, in the non-rotating shear case, the expected exponential growth of turbulent kinetic energy is found only with a hybrid model for nonlinear terms combining the anisotropic EDQNM closure and Weinstock’s return-to-isotropy model.