def get_verts(voxel, g):
    """return list (len=8) of point coordinates (x,y,z) that are vertices of the voxel (i,j,k)"""
    i, j, k = voxel
    dx, dy, dz = g['dx'], g['dy'], g['dz']
    v1_0, v1_1, v1_2 = g['xlo'] + i * dx, g['ylo'] + j * dy, g['zlo'] + k * dz
    vertices = [(v1_0, v1_1, v1_2), (v1_0 + dx, v1_1, v1_2), (v1_0 + dx, 
        v1_1 + dy, v1_2), (v1_0, v1_1 + dy, v1_2), (v1_0, v1_1, v1_2 + dz),
        (v1_0 + dx, v1_1, v1_2 + dz), (v1_0 + dx, v1_1 + dy, v1_2 + dz), (
        v1_0, v1_1 + dy, v1_2 + dz)]
    return vertices
