Алгоритмы вычисления параметров магнитного поля Земли с испо

Существуют три наиболее используемые модели геомагнитного поля Земли:
1. IGRF (International Geomagnetic Reference Field) – разработано Международной Ассоциацией Геомагнетизма и аэрономии (The International Association of Geomagnetism and Aeronomy, IAGA, http://www.ngdc.noaa.gov/IAGA/vmod/igrf.html)
2. WMM (World Magnetic Model) – всемирная модель геомагнитного поля Земли - разработана Национальным агентством США по обработке пространственных геоданных (The U.S. National Geospatial-Intelligence Agency, NGA), Национальным центром геофизических данных США (The U.S. National Geophysical Data Center, NGDC, http://www.ngdc.noaa.gov/ngdc.html) и Британским геологоразведочным обществом (the British Geological Survey, BGS, http://www.geomag.bgs.ac.uk/)
3. Глобальная модель геомагнитного поля - BGGM (The BGS Global Geomagnetic Model – разработана Британским геологоразведочным обществом (the British Geological Survey, BGS, http://www.geomag.bgs.ac.uk/).

На просторах Интернета можно без труда вычислить параметры магнитного поля Земли в любой точке, используя ту или иную модель (IGRF, WMM, BGGM) как в режиме on-line (например, http://www.geomag.bgs.ac.uk/gifs/igrf_form.shtml), так и используя различные программы (GeoMagix) без выхода в сеть.

Также в сети можно найти коэффициенты моделей, по которым рассчитываются искомые параметры на интересующую нас эпоху (5 лет) и математику расчета. К сожалению, практически вся информация, мною найденная на английском с использованием довольно трудных формул для моего понимания.
Хотел бы запрограммировать оф-лайн программу на C или Pascale по расчету параметров магнитного поля Земли с использованием глобальных моделей геомагнитного поля Земли. Для этого мне хотелось бы узнать:
1. Одиннаковые ли алгоритмы вычисления параметров в различных моделях. Если алгоритм один и тот же, то просто подставляем нужные коэффициенты. И все.
2. На пальцах объяснить этот алгоритм. Мое представление – есть некий грид, по вершинам которого расставлены коэффициенты, производим мат функцию внутри грида и получаем результат.
Вообще я нашел и исходный код на С и Fortran для расчета, но он только для IGRF модели ( http://www.ngdc.noaa.gov/IAGA/vmod/igrf.html). Как для WMM и BGGM? На вышеуказанных сайтах тоже много информации. Может пообсуждаем и я услышу дельные советы?

 

kostyarico, получилась программа?