Improved understanding of the impact of crystal mush rheology on the response of magma chambers to magmatic events is critical for better understanding crustal igneous systems with abundant crystals. This project is an extension of Liao et al (2018), which models the response of a partially-mushy magma chamber with poroviscoelastic mush under injection of new magma. Analytical solutions are realized numerically with the codes provided.