12 include
'geodesic.inc'
14 double precision a, f, lat1, lon1, azi1, lat2, lon2, azi2, s12,
27 read(*, *, end=90, err=90) lat1, lon1, azi1, s12
28 call
direct(a, f, lat1, lon1, azi1, s12, arcmod,
29 + lat2, lon2, azi2, omask, dummy, dummy, dummy, dummy, dummy)
30 print 20, lat2, lon2, azi2
31 20
format(f20.15, 1x, f20.15, 1x, f20.15)