added functionality to include zs array (only implemented for initialize so far)
???? Please see comment about rho weighting... Why are the weights based on the RA field?????