Comment (by Lorenzo Sala):

Thank you so much for pointing this out, this has fixed everything! Fixing those two parameters, rho and temperature do not get stuck and the plot perfectly overlaps the non-MHD one.

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2272/grhydro-c2p-inconsistency-with-and-without