int bond_off = 0;

  int angle_off = 0;

  for (i = 0; i < modify->nfix; i++)

    if (strcmp(modify->fix[i]->style,"shake") == 0)

    if ((strcmp(modify->fix[i]->style,"shake") == 0)

        || (strcmp(modify->fix[i]->style,"rattle") == 0))

      bond_off = angle_off = 1;

  if (force->bond && force->bond_match("quartic")) bond_off = 1;

  // extract count from buffer that has cycled back to me

  // nspecial[i][1] = # of 1-3 neighbors of atom i

  j = 0;

<  j = 0;

  for (i = 0; i < nlocal; i++) {

    nspecial[i][1] = buf[j];

    j += 2 + nspecial[i][0];