added ability to handle mean rho(z)