code: def shift_com(elements, coordinates, com_adjust=np.zeros(3)):
      
      com = center_of_mass(elements, coordinates)
      com = np.array([com - com_adjust] * coordinates.shape[0])
      return coordinates - com
docs: Return coordinates translated by some vector.
  
      Parameters
      ----------
      elements : numpy.ndarray
          An array of all elements (type: str) in a molecule.
  
      coordinates : numpy.ndarray
          An array containing molecule's coordinates.
  
      com_adjust : numpy.ndarray (default = [0, 0, 0])
  
      Returns
      -------
      numpy.ndarray
          Translated array of molecule's coordinates.