A unified model of multi-GNSS and multi‑frequency precise point positioning for the joint estimation of ionospheric TEC and time-varying receiver code bias