The solar differential rotation and meridional flow as computed with the new scheme. Top left: surface rotation rate versus colatitude. Top right: rotation rate as function of radius at 0, 15, 30, 45, 60, 75, and 90° latitude, respectively, from top to bottom. Bottom left: streamlines of the meridional flow. The flow is counter-clockwise, i.e. directed towards the pole at the surface and towards the equator at the bottom. Bottom right: flow speed as a function of latitude at the top (solid blue) and bottom (dashed red) of the convection zone. Positive values indicate a flow away from the north pole.