onemols = &atom->molecules[imol];

  nset = atom->molecules[imol]->nset;

  // error check on molecule template fields

  for (int i = 0; i < nset; i++) {

    if (onemols[i]->typeflag == 0)

      error->all(FLERR,"Atom style template molecule must have atom types");

  // set bonds_allow,angles_allow,etc based on the molecules in template set

  // similar to how atom_style bond,angle,full set it