DIFFERENTIAL ROTATION AND TURBULENT CONVECTION : A NEW REYNOLDS STRESS MODEL AND COMPARISON WITH SOLAR DATA