Model 4A from K. Arredondo and M. I. Billen. The effects of phase transitions and compositional layering in two- dimensional kinematic models of subduction. Journal of Geodynamics, 100:159–174, 2016.