https://github.com/cran/SoDA
Tip revision: dc000f2a5f006d137f66716b086025d618bf8306 authored by John M Chambers on 14 July 2008, 00:00:00 UTC
version 1.0-5
version 1.0-5
Tip revision: dc000f2
geodistv.f
subroutine geodistv(lat1,lon1,lat2,lon2,dist, n)
double precision lat1(1), lon1(1), lat2(1), lon2(1), dist(1)
integer n
double precision faz, baz
double precision a,f
integer i
a = 6378137.00
f = 1./298.257223563
do 200 i=1,n
call geoddist(lat1(i), lon1(i), lat2(i), lon2(i), a, f,
1 faz, baz, dist(i))
200 continue
end