alb_v[ind] += dzv[ind] * (1.0 - cosz[ind])