  QuaLiKiz, a model based on a local gyrokinetic eigenvalue solver is expanded
to include momentum flux modeling in addition to heat and particle fluxes.
Essential for accurate momentum flux predictions, the parallel asymmetrization
of the eigenfunctions is successfully recovered by an analytical fluid model.
This is tested against self-consistent gyrokinetic calculations and allows for
a correct prediction of the ExB shear impact on the saturated potential
amplitude by means of a mixing length rule. Hence, the effect of the ExB shear
is recovered on all the transport channels including the induced residual
stress. Including these additions, QuaLiKiz remains ~10 000 faster than
non-linear gyrokinetic codes allowing for comparisons with experiments without
resorting to high performance computing. The example is given of momentum pinch
calculations in NBI modulation experiments.
