I implemented a solver for non-ideal MHD including super-time-stepping in Arepo. I use this solver to perform 3D, global simulations of protoplanetary discs to understand the interaction of gravitational instability, MRI and non-ideal MHD.