code: def calculate_angular_momentum(self):
          
          clibrebound.reb_tools_angular_momentum.restype = reb_vec3d
          L = clibrebound.reb_tools_angular_momentum(byref(self))
          return [L.x, L.y, L.z]
docs: Returns a list of the three (x,y,z) components of the total angular momentum of all particles in the simulation.