Longitudinal variation with latitude (R script)

 

lat  = 0
long_dd_km <- 6378.145 * 2 * atan2(sqrt(cos(lat * pi/180)^2 * 7.615242e-05), sqrt(1 - cos(lat * pi/180)^2 * 7.615242e-05))