Quad tree indexing
ix = qnode(d,bx,by)
depth in splitting
box coordinates at that depth
index of that node in tree
ix = 4^(d) + by + bx*2^(d)