function llDistance ( $lat1, $lon1, $lat2, $lon2, $unit='M' )
        {
#-k     geo
#-      this routine calculates the distance between two points (given the
#-      latitude/longitude of those points).
#-
#-       unit:  'M' is statute miles (default)
#-              'K' is kilometers
#-              'N' is nautical miles

        $theta = $lon1 - $lon2;
        $dist = sin(deg2rad($lat1)) * sin(deg2rad($lat2)) +  cos(deg2rad($lat1)) * cos(deg2rad($lat2)) * cos(deg2rad($theta));
        $dist = acos($dist);
        $dist = rad2deg($dist);
        $miles = $dist * 60 * 1.1515;
        $unit = strtoupper($unit);

        if ($unit == "K")       return ($miles * 1.609344);
        elseif ($unit == "N")   return ($miles * 0.8684);
        else                    return $miles;

        }

?>