Panel (a) shows colour-filled contours of solar differential rotation parametrized from helioseismic observations by Charbonneau et al. (1997), steady-state meridional circulation patterns produced from our model for density function exponent m = 0.2, 0.5 and 1.5 are presented in panels (b), (c) and (d).