Text this: High-Order Compact Implicit Difference Methods For Parabolic Equations in Geodynamo Simulation