zb_quadrat function

Divide a region into quadrats