Sekian lama tidak diupdate, kali ini saya akan berbagi cara membuat peta kontur intensitas medan magnet dari model IGRF (International Geomagnetic Reference Field). Terinspirasi dari pertanyaan kawan yang menanyakan cara mengambil data IGRF di beberapa lokasi sekaligus, saya kembangkan menjadi cara membuat peta kontur intensitas medan magnetnya. Perlu diketahui bahwa model IGRF diupdate setiap 5 tahun sekali. Di GMT 5 (GMT 5.1.0) telah tersedia data IGRF dari tahun 1900-2015. Sehingga data yang bisa digunakan hanya sampai tanggal 1 Januari 2015. Di versi terbaru dari GMT 5 (GMT 5.1.2) mungkin hal ini telah diupdate. Saya sendiri belum mencoba GMT 5.1.2.
Disini kita menggunakan perintah "mgd77magref" untuk mengambil data IGRF. Contoh syntax modul ini adalah:
echo 98 3.5 | mgd77magref -A+a1+t1989-10-29T10:00:00 -Frt/0
98 3.5 41464.6232039