Magma dynamics using FD-PDE: a new, PETSc-based, finite-difference staggered-grid framework for solving partial differential equations