1  SELECT CASE (lp)
2  CASE (0)
3     CALL integrate_core_0(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
4                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
5  CASE (1)
6     CALL integrate_core_1(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
7                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
8  CASE (2)
9     CALL integrate_core_2(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
10                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
11  CASE (3)
12     CALL integrate_core_3(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
13                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
14  CASE (4)
15     CALL integrate_core_4(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
16                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
17  CASE (5)
18     CALL integrate_core_5(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
19                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
20  CASE (6)
21     CALL integrate_core_6(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
22                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
23  CASE (7)
24     CALL integrate_core_7(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
25                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
26  CASE (8)
27     CALL integrate_core_8(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
28                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
29  CASE (9)
30     CALL integrate_core_9(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
31                           map(-cmax, 1), sphere_bounds(1), cmax, gridbounds(1, 1))
32  CASE DEFAULT
33     CALL integrate_core_default(grid(1, 1, 1), coef_xyz(1), pol_x(0, -cmax), pol_y(1, 0, -cmax), pol_z(1, 0, -cmax), &
34                                 map(-cmax, 1), sphere_bounds(1), lp, cmax, gridbounds(1, 1))
35  END SELECT
36