def azidist(v1, v2, zenith=True):
    if (v1.shape[0] == 3):
        v1 = vec2sph(v1, zenith=zenith)
    if (v2.shape[0] == 3):
        v2 = vec2sph(v2, zenith=zenith)
    d = ((((v1[1] - v2[1]) + np.pi) % (2 * np.pi)) - np.pi)
    return np.absolute(d)
def azidist(v1, v2, zenith=True):
    if (v1.shape[0] == 3):
        v1 = vec2sph(v1, zenith=zenith)
    if (v2.shape[0] == 3):
        v2 = vec2sph(v2, zenith=zenith)
    d = ((((v1[1] - v2[1]) + np.pi) % (2 * np.pi)) - np.pi)
    return np.absolute(d)