anchor.diff:                  write (6,'()')
               end if
   
  -        end if
  +        end select
   
       end subroutine init_symmetry
   
  
positive.code:         ! Construct the symmetry tables.
  
          use basis, only: nbasis, basis_fns, write_basis_fn
          use system, only: ndim, system_type, hub_real
          use kpoints, only: is_reciprocal_lattice_vector
          use parallel, only: parent
          use utils, only: int_fmt
  
negative.code:         integer :: ksum(ndim)
          character(4) :: fmt1
  
          select case(system_type)
          case(hub_real)
  
              ! Symmetry currently not implemented.
              nsym = 1
  
          case(hub_k)
  
              ! Each wavevector corresponds to its own irreducible representation.
              ! The maximum system size considered will be quite small (ie <200)
              ! and the sum of any two wavevectors is another wavevector in the
              ! basis (up to a primitive reciprocal lattice vector).
              ! It is thus feasible to store the nsym^2 product table.
  
              nsym = nbasis/2
              allocate(sym_table(nsym, nsym), stat=ierr)
  