example-MagneticCircle.cs 1.57 KB
Newer Older
Valentin Platzgummer's avatar
Valentin Platzgummer committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
using System;
using NETGeographicLib;

namespace example_MagneticCircle
{
    class Program
    {
        static void Main(string[] args)
        {
            try {
                MagneticModel mag = new MagneticModel("wmm2010","");
                double lat = 27.99, lon0 = 86.93, h = 8820, t = 2012; // Mt Everest
                {
                    // Slow method of evaluating the values at several points on a circle of
                    // latitude.
                    for (int i = -5; i <= 5; ++i) {
                        double lon = lon0 + i * 0.2;
                        double Bx, By, Bz;
                        mag.Field(t, lat, lon, h, out Bx, out By, out Bz);
                        Console.WriteLine(String.Format("{0} {1} {2} {3}", lon, Bx, By, Bz));
                    }
                }
                {
                    // Fast method of evaluating the values at several points on a circle of
                    // latitude using MagneticCircle.
                    MagneticCircle circ = mag.Circle(t, lat, h);
                    for (int i = -5; i <= 5; ++i) {
                        double lon = lon0 + i * 0.2;
                        double Bx, By, Bz;
                        circ.Field(lon, out Bx, out By, out Bz);
                        Console.WriteLine(String.Format("{0} {1} {2} {3}", lon, Bx, By, Bz));
                    }
                }
            }
            catch (GeographicErr e) {
                Console.WriteLine(String.Format("Caught exception: {0}", e.Message));
            }
        }
    }
}