code: def dtraj_T100K_dt10_n(self, divides):
          
          disc = np.zeros(100, dtype=int)
          divides = np.concatenate([divides, [100]])
          for i in range(len(divides)-1):
              disc[divides[i]:divides[i+1]] = i+1
          return disc[self.dtraj_T100K_dt10]
docs: 100K frames trajectory at timestep 10, arbitrary n-state discretization.