estimate within grain variability in data code
rework Zhang method a bit, extend data example