Mathematical modeling of COVID-19 in British Columbia: An age-structured model with time-dependent contact rates