This is based on the following paper: Sośnica, K., Bury, G., Zajdel, R. et al. General relativistic effects acting on the orbits of Galileo satellites. Celest Mech Dyn Astr 133, 14 (2021). https://doi.org/10.1007/s10569-021-10014-y
If you happen to find any implementation mistakes, please let me know.