GNSS ECEF from SP3 to Eastward velocity in measurement plane for a given receiver location

[himanshu10140]

Flowchart

1. User Input of ecef x,y,z from sp3 file for GPS/GALILEO/GLONASS/BEIDOU/QZSS
2. Calculation of Azimuth and Elevation for given receiver location
3. Satellite velocities (NED) at Ionospheric Pierce Point(350 km) in propagation system
4. Eastward apparent velocity in measurement plane in geographic and magnetic coordinate.

I have implemented the above for geographic coordinate system with success but my velocities in magnetic coordinates are very high.
Hence, if you can implement this feature, it would be helpful for utilization in zonal drift velocity calculations using GNSS observations.

References

1. Appendix (A2) of 10.1002/2014RS005513 (Ghafoori and Skone)
2. Rino 2011 https://www.mathworks.com/matlabcentral/fileexchange/28888-satellite-orbit-computation?s_tid=prof_contriblnk
3. 10.1002/2015RS005864 (Carrano et. al. 2016)