code: def run_mesh(self,
                   mesh=100.0,
                   shift=None,
                   is_time_reversal=True,
                   is_mesh_symmetry=True,
                   with_eigenvectors=False,
                   with_group_velocities=False,
                   is_gamma_center=False):
          
  
          self.init_mesh(mesh=mesh,
                         shift=shift,
                         is_time_reversal=is_time_reversal,
                         is_mesh_symmetry=is_mesh_symmetry,
                         with_eigenvectors=with_eigenvectors,
                         with_group_velocities=with_group_velocities,
                         is_gamma_center=is_gamma_center)
          self._mesh.run()
docs: Run mesh sampling phonon calculation.
  
          See the parameter details in Phonopy.init_mesh().