Our model is an age-dependent compartmental model. The model flowchart is given in Fig. 1 , where the subscript i denotes the i th age group. The model variables are described as follows, with the time unit t in days: