void FillLocal(int *iclos,int itri)
{
   int k,res1,res2;

// FillLocal fills iclos with atoms that belong to the two residues
// that define the rotational axes.

   if (DCKMETH == 0) {

	res1 = inttyp[TriArr[itri]->prot[0]]->res;
	res2 = inttyp[TriArr[itri]->prot[1]]->res;

	for (k=0; k<*xyzp->iatoms; k++) {
	   iclos[k] = 0;
	}
	for (k=0; k<*xyzp->iatoms; k++) {
	   if (xyzp->iresid[k] == res1 || xyzp->iresid[k] == res2) {
	      iclos[k] = 1;
	   }
	}
   }

}

