Name

---

      ELL_GEO_SIDE

Purpose

---

      Find which side of a geodesic a given point located on.

Category

---

Calling Sequence

---

      side = ell_geo_side( p1, p2, pt)

Inputs

---

      p1 = Start point of the geodesic {lon:lon1, lat:lat1}. in
      p2 = End point of the geodesic {lon:lon2, lat:lat2}. in
      pt = Point to check {lon:lont, lat:latt}. in
          If lont and latt are arrays side will be an array also.
      side = Side code (pt as seen from p1):
        -1 means pt is left of p2.
          0 means pt is aligned with p2.
        +1 means pt is right of p2.

Keyword Parameters

---

Outputs

---

Common Blocks

---

Notes

---

Modification History

---

      R. Sterner, 2006 Sep 25
      R. Sterner, 2009 Jan 15 --- Allowed arrays in pt.
  Copyright (C) 2006, Johns Hopkins University/Applied Physics Laboratory
  This software may be used, copied, or redistributed as long as it is not
  sold and this copyright notice is reproduced on each copy made. This
  routine is provided as is without any express or implied warranties
  whatsoever. Other limitations apply as described in the file disclaimer.txt.