Membuat Peta Kontur Model IGRF

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

Arti dari perintah diatas adalah mengambil data IGRF untuk komponen Total (t) pada ketinggian 1 km tanggal 29 Oktober 1989 di bujur 98 dan lintang 3.5. Bagian -F digunakan untuk memilih komponen dan jenis model. Pilihan komponen yang lain adalah :

Deklinasi (d, derajat, positif searah jarum jam)

Inklinasi (i, derajat, positif ke atas)

Horizontal (h, nT)

komponen X (x, nT, positif ke Utara)

komponen Y (y, nT, positif ke Timur)

komponen Z (z, nT, positif ke bawah)

komponen Total (t, nT).

Bagian r pada -F artinya kita tampilkan koordinatnya.

Nilai 0 pada -F artinya kita menggunakan model IGRF. Pilihan yang lain adalah :

1 CM4 Core field
2 CM4 Lithospheric field
3 CM4 Primary Magnetospheric field
4 CM4 Induced Magnetospheric field
5 CM4 Primary ionospheric field
6 CM4 Induced ionospheric field
7 CM4 Toroidal field
9 Core field dari IGRF model lainnya dari CM4. Jangan gunakan 0 bersama dengan 9.

Bagian -A+a1+t1989-10-29T10:00:00 artinya kita menggunakan model IGRF pada ketinggian 1 km (+a1) tanggal 29 Oktober 1989 (+t1989-10-29) dan jam 10:00:00 UTC.

Selanjutnya kita akan menggabungkan perintah mgd77magref dan memplotnya di peta. Kita menggunakan perintah "for" untuk membuat data koordinat titik-titik gridnya. Berikut adalah scriptnya:

#File output
set F=igrf.ps
#Batas kiri
set lon_start=94
#Batas kanan
set lon_stop=141
#Batas bawah
set lat_start=-12
#Batas atas
set lat_stop=9
#Ketinggian referensi
set altitude=1
#Tanggal
set tanggal=1989-10-29
#Jam
set jam=10:00:00

@echo off
(FOR /L %%H IN (%lon_start%,1,%lon_stop%) DO (FOR /L %%G IN (%lat_start%,1,%lat_stop%) DO echo %%H %%G)) > datamagnet

@echo on
mgd77magref datamagnet -A+a%altitude%+t%tanggal%T%jam% -Frt/0 > igrf.dat

#Plot nilai medan magnet total (kolom ke-3)
gawk "{print $1, $2, $3}" igrf.dat | blockmean -R%lon_start%/%lon_stop%/%lat_start%/%lat_stop% -I0.1 | surface -R%lon_start%/%lon_stop%/%lat_start%/%lat_stop% -I0.1 -Gigrf.grd

makecpt -Z -Cglobe > elev.cpt

grdimage -JM9.5i -B4:."PETA INTENSITAS MEDAN MAGNET TOTAL MODEL IGRF":WSne -R%lon_start%/%lon_stop%/%lat_start%/%lat_stop% -K -Celev.cpt indo.nc > %F%

pscoast -K -W1 -Dh -JM -R -O >> %F%

#Plot kontur dengan interval 400nT dan label dengan interval 1200 nT
grdcontour igrf.grd -R -JM -O -C400 -A1200+f10+gwhite+u" nT" -W2,yellow >> %F%

#Konversi dari file postscript ke PNG
ps2raster %F% -A -TG -P

del datamagnet igrf.grd mag.cpt

Setelah dieksekusi hasilnya akan seperti berikut.

Sekian untuk kali ini. Selamat mencoba dengan komponen-komponen medan magnet lainnya. Salam.