Two-temperature treatments in magnetically arrested disk GRMHD simulations more accurately predict light curves of Sagittarius A*

Reproduction package for Salas, Liska, Markoff et al. 2025, MNRAS"Two-temperature treatments in magnetically arrested disk GRMHD simulations more accurately predict light curves of Sagittarius A*"Paper: https://doi.org/10.1093/mnras/staf240 The GRMHD simulation npz files in the folder "GRMHD", are split up into subfolders: MAD2TC_R0_27m20: simulation 2TC in table 1 of Salas et al. 2025MAD2TNC_R0_27m20: simulation 2T in table 1 of Salas et al. 2025MADNC: simulation 1T in table 1 of Salas et al. 2025 The python jupyter notebooks:  - BHOSS_raytracing_2TC.ipynb to read and plot the BHOSS postprocess data: synchrotron light curves at 43, 86, 214, 228, 230, 345 and 1360 GHz- 2TSgrA_pyhamrds.ipynb to read and plot the H-AMR portprocess data: GRMHD variables in spacetime The figures are saved in "Output". The notebook was built using Python 3.9.18. A description of the files are as follows: "Fluxes1T.npz": synchrotron fluxes of single temperature (1T) simulations "Fluxes2T.npz": synchrotron fluxes of two temperature (2T) simulations "postprocess.npz": GRMHD variables in spacetime