12 include
'geodesic.inc'
14 double precision a, f, lat1, lon1, azi1, lat2, lon2, azi2, s12,
26 read(*, *, end=90, err=90) lat1, lon1, azi1, s12
27 call
direct(a, f, lat1, lon1, azi1, s12, flags,
28 + lat2, lon2, azi2, omask, dummy, dummy, dummy, dummy, dummy)
29 print 20, lat2, lon2, azi2
30 20
format(f20.15, 1x, f20.15, 1x, f20.15)
program geoddirect
A simple program to solve the direct geodesic problem.