We show our work to refine the process of evolutions in general relativistic magnetohydrodynamics. We investigate several areas in order to improve the overall accuracy of our results. We test several versions of conversion methodologies between different sets of variables. We compare both single equation and two equations solvers to do the conversion. We find no significant improvement for multiple equation conversion solvers when compared to single equation solvers. We also investigate the construction of initial data and the conversion of coordinate systems between initial data code and evolution code. In addition to the conversion work, we have improved some methodologies to ensure data integrity when moving data from the initial data code to the evolution code. Additionally we add into the system of MHD equations a new field to help control the no monopole constraint. We perform a characteristic decomposition of the system of equations in order to derive the associated boundary condition for this new field. Finally, we implement a WENO (weighted non-oscillatory) system. This is done so we can evolve and track shocks that are generated during an evolution of our GRMHD equations.