12 subroutine direct(a, f, lat1, lon1, azi1, s12a12, flags,
13 + lat2, lon2, azi2, omask, a12s12, m12, MM12, MM21, SS12)
14 double precision,
intent(in) :: a, f, lat1, lon1, azi1, s12a12
15 integer,
intent(in) :: flags, omask
16 double precision,
intent(out) :: lat2, lon2, azi2
17 double precision,
intent(out) :: a12s12, m12, MM12, MM21, SS12
20 subroutine invers(a, f, lat1, lon1, lat2, lon2,
21 + s12, azi1, azi2, omask, a12, m12, MM12, MM21, SS12)
22 double precision,
intent(in) :: a, f, lat1, lon1, lat2, lon2
23 integer,
intent(in) :: omask
24 double precision,
intent(out) :: s12, azi1, azi2
25 double precision,
intent(out) :: a12, m12, MM12, MM21, SS12
28 subroutine area(a, f, lats, lons, n, AA, PP)
29 integer,
intent(in) :: n
30 double precision,
intent(in) :: a, f, lats(n), lons(n)
31 double precision,
intent(out) :: AA, PP