12 subroutine direct(a, f, lat1, lon1, azi1, s12a12, arcmod,
13 + lat2, lon2, azi2, omask, a12s12, m12, mm12, mm21, ss12)
14 double precision,
intent(in) :: a, f, lat1, lon1, azi1, s12a12
15 logical,
intent(in) :: arcmod
16 integer,
intent(in) :: omask
17 double precision,
intent(out) :: lat2, lon2, azi2
18 double precision,
intent(out) :: a12s12, m12, mm12, mm21, ss12
21 subroutine invers(a, f, lat1, lon1, lat2, lon2,
22 + s12, azi1, azi2, omask, a12, m12, mm12, mm21, ss12)
23 double precision,
intent(in) :: a, f, lat1, lon1, lat2, lon2
24 integer,
intent(in) :: omask
25 double precision,
intent(out) :: s12, azi1, azi2
26 double precision,
intent(out) :: a12, m12, mm12, mm21, ss12
29 subroutine area(a, f, lats, lons, n, AA, PP)
30 integer,
intent(in) :: n
31 double precision,
intent(in) :: a, f, lats(n), lons(n)
32 double precision,
intent(out) :: aa, pp