Different authors use sub-ionospheric latitude and longitude to to express latitudinal gradient of TEC, Das et al., 2014. How can we find the GPS tracks in terms of sub ionospheric latitude and longitude ?
lat is the latitude of ionospheric pierce point (IPP) in radians.
lon is the longitude of ionospheric pierce point in radians.
Oblique is the angle the ray path makes with the assumed ionospheric thin shell.
p is the angle subtended at earth center.
!PI/180 factors are used for degrees to radian conversion, since these lines of code come from my IDL scripts (which operate in radians). Remove these factors in all locations or convert to degrees if your coding language uses degrees.
Be careful, some authors will reference an equation for longitude that incorrectly uses the station latitude instead of the IPP latitude. This is due to a typo in one of the older texts in the field.
David R. Themens Could you confirm as to whether the Asin(a) factor in the equation for p stands for arcsin(a)? Also, the asin() factors in the equations for lat and lon (IPP latitude and longitude) too stands for arcsin() right? Its NOT a (= R*COS(Elev_use*!Pi/180.0)/(R+h)) multiplied with sin().