[variance] r/sigma_z/variance/g

and compute real variance in variance_mean()