end, we selected the model proposed by Wonham and colleagues [16] to compute the dynamics of the two populations. This model is based on 8 differential equations (Figure 12a ) which can compute over time the evolution of the different categories of individuals (called 'compartments'): susceptible,infected, recovered and dead birds, the larvae of mosquitoes and the susceptible, exposed and infected adult mosquitoes. We proposed some modifications in order to correct some discrepancies that we found in the model. We also included climate effects in the model using the work of Madder and colleagues [33] as a starting point. This was a difficult task because the model was no longer in equilibrium and this required several modifications to the differential equations [18] . Figure  12b presents an overview of the new equations of the proposed model. We notice that the birds equations (dS BJ /dt, dI BJ /dt, dR BJ /dt, dX BJ /dt) have an index j which represents a different bird species that we want to include in the simulation. Climate effects are computed using another set of equations that is not presented in this paper. These equations modify certain parameters in the differential equations of Figure 12b . For example v xL in the dL M /dt equation represents the rate of progress of larvae toward the state of nymphae, this rate depends on temperature conditions defined in a different set of equations. The adjusted model gives satisfactory results in terms of quality (e.g. distribution of the mosquitoes' generations). Indeed, the pace of the established curves reflects the biological behaviours of the studied species if we refer to the specialized literature. However, the quantitative results (e.g. the number of larvae, eggs, emerged Culex, dead crows, etc.) had not been conclusive with the first results of the simulation. We corrected this problem with the calibration of the system as presented in Section 2.4.