geoddirect.for File Reference

Go to the source code of this file.

Functions/Subroutines
*file geoddirect for *brief A test program for	direct ()* > A simple program to solve the direct geodesic problem.*!*!This program reads in lines with lat1
*file geoddirect for *brief A test program for s12 and prints out *!lines with	azi2 (for the WGS84 ellipsoid).program geoddirect implicit none include 'geodesic.inc'double precision a
*file geoddirect for *brief A test program for s12 and prints out *!lines with dummy logical arcmod integer omask *WGS84 values s12 call	direct (a, f, lat1, lon1, azi1, s12, arcmod,+lat2, lon2, azi2, omask, dummy, dummy, dummy, dummy, dummy) print 20
Variables
*file geoddirect for *brief A test program for	lon1
*file geoddirect for *brief A test program for	azi1 = 0&deg
*file geoddirect for *brief A test program for s12 and prints out *!lines with	lat2
*file geoddirect for *brief A test program for s12 and prints out *!lines with	lon2
*file geoddirect for *brief A test program for s12 and prints out *!lines with	f
*file geoddirect for *brief A test program for s12 and prints out *!lines with	lat1
*file geoddirect for *brief A test program for s12 and prints out *!lines with	azi2
*file geoddirect for *brief A test program for s12 and prints out *!lines with	s12
*file geoddirect for *brief A test program for s12 and prints out *!lines with dummy logical arcmod integer omask *WGS84 values	a

---

Function Documentation

* file geoddirect for* brief A test program for direct

(

)

* file geoddirect for* brief A test program for s12 and prints out* ! lines with azi2

(

for the WGS84

ellipsoid

)

Type Constraints

* file geoddirect for* brief A test program for s12 and prints out* ! lines with dummy logical arcmod integer omask* WGS84 values s12 call direct	(	a	,
		f	,
		lat1	,
		lon1	,
		azi1	,
		s12	,
		arcmod	,
		+	lat2,
		lon2	,
		azi2	,
		omask	,
		dummy	,
		dummy	,
		dummy	,
		dummy	,
		dummy
	)

Type Constraints

---

Variable Documentation

*file geodinverse for *brief A test program for lat2 and prints *out lines with dummy integer omask *WGS84 values lon1

Definition at line 6 of file geoddirect.for.

*file geodinverse for *brief A test program for lat2 and prints *out lines with dummy integer omask *WGS84 values lon2 call azi1 = 0&deg

Definition at line 6 of file geoddirect.for.

*file geodinverse for *brief A test program for lat2 and prints *out lines with dummy integer omask *WGS84 values lat2

Definition at line 6 of file geoddirect.for.

*file geodinverse for *brief A test program for lat2 and prints *out lines with lon2

Definition at line 6 of file geoddirect.for.

* file geoddirect for* brief A test program for s12 and prints out* ! lines with f

Definition at line 14 of file geoddirect.for.

* file geoddirect for* brief A test program for s12 and prints out* ! lines with lat1

Definition at line 14 of file geoddirect.for.

*file geodinverse for *brief A test program for lat2 and prints *out lines with dummy integer omask *WGS84 values lon2 call azi2

Initial value:

*!   [\e azi1, \e azi2] + [\e d, &minus

Definition at line 14 of file geoddirect.for.

*differential and integral properties of geodesics are computed **The shortest path between two points on the ellipsoid e lon2 is called the geodesic Its length is *e s12 and the geodesic from point to point has forward azimuths *e azi1 and e azi2 at the two end points **Traditionally two geodesic problems are e e and e *determine e e and e azi2 This is solved by the *subroutine e e e *determine e s12

Definition at line 14 of file geoddirect.for.

* file geoddirect for* brief A test program for s12 and prints out* ! lines with dummy logical arcmod integer omask* WGS84 values a

Initial value:

 6378137d0
      f = 1/298.257223563d0

      arcmod = .false.
      omask = 0

 10   continue
      read(*, *, end=90, err=90) lat1

Definition at line 20 of file geoddirect.for.