HL, CW, HH, and XL developed the model. HL and CW do the computational work. HL wrote the manuscript.