anchor.diff:  
           select case(sys%system)
           case(hub_k, hub_real, heisenberg, chung_landau, ringium)
  -            if (ibasis /= nspatial) call stop_all('init_basis_fns','Not enough basis functions found.')
  +            if (ibasis /= nspatial) call stop_all(this,'Not enough basis functions found.')
           case(ueg)
           end select
   
  
positive.code: 
          type(sys_t), intent(inout) :: sys
          logical, intent(in), optional :: store_info
  
          logical :: t_store
  
  
negative.code:                         kp = (/ i, j, k /)
                          if (in_FBZ(sys%system, sys%lattice, kp(1:sys%lattice%ndim))) then
                              if (ibasis == nspatial) then
                                  call stop_all(this,'Too many basis functions found.')
                              else
                                  ! Have found an allowed wavevector/site.
                                  ! Add 2 spin orbitals to the set of the basis functions.
  