Adaptive numerical integration on spherical triangles