finished counting intersitials in get_vacancy.m