using System;
using NETGeographicLib;

namespace example_Geocentric
{
    class Program
    {
        static void Main(string[] args)
        {
            try {
                Geocentric earth = new Geocentric( Constants.WGS84.EquatorialRadius,
                                                   Constants.WGS84.Flattening);
                // Alternatively: Geocentric earth = new Geocentric();
                {
                    // Sample forward calculation
                    double lat = 27.99, lon = 86.93, h = 8820; // Mt Everest
                    double X, Y, Z;
                    earth.Forward(lat, lon, h, out X, out Y, out Z);
                    Console.WriteLine( String.Format( "{0} {1} {2}",
                        Math.Floor(X / 1000 + 0.5),
                        Math.Floor(Y / 1000 + 0.5),
                        Math.Floor(Z / 1000 + 0.5) ) );
                }
                {
                    // Sample reverse calculation
                    double X = 302e3, Y = 5636e3, Z = 2980e3;
                    double lat, lon, h;
                    earth.Reverse(X, Y, Z, out lat, out lon, out h);
                    Console.WriteLine(String.Format("{0} {1} {2}", lat, lon, h));
                }
            }
            catch (GeographicErr e) {
                Console.WriteLine( String.Format( "Caught exception: {0}", e.Message ) );
            }
        }
    }
}
