Menampilkan Batas-Batas Lempeng Tektonik Dengan GMT

Pada artikel kali ini, saya akan berbagi cara menampilkan batas-batas lempeng tektonik dengan GMT. Sebelumnya, download dahulu data-data batas lempeng tersebut disini. Anda juga dapat menambahkan data fracture zone disini. Setelah didownload, pastikan data-data tersebut berekstensi ".gmt". Selanjutnya downloadlah data topografi disini. 

Kemudian buatlah file warna seperti berikut :

  -10000    153    0    255    -9500    153    0    255
-9500    153    0    255    -9000    153    0    255
-9000    153    0    255    -8500    153    0    255
-8500    136    17    255    -8000    136    17    255
-8000    119    34    255    -7500    119    34    255
-7500    102    51    255    -7000    102    51    255
-7000    85    68    255    -6500    85    68    255
-6500    68    85    255    -6000    68    85    255
-6000    51    102    255    -5500    51    102    255
-5500    34    119    255    -5000    34    119    255
-5000    17    136    255    -4500    17    136    255
-4500    0    153    255    -4000    0    153    255
-4000    27    164    255    -3500    27    164    255
-3500    54    175    255    -3000    54    175    255
-3000    81    186    255    -2500    81    186    255
-2500    108    197    255    -2000    108    197    255
-2000    134    208    255    -1500    134    208    255
-1500    161    219    255    -1000    161    219    255
-1000    188    230    255    -500    188    230    255
-500    215    241    255    -200    215    241    255
-200    241    252    255    0    241    252    255
  0      0       97      71      50      16      122     47
  50     16      122     47      500     232     215     125
  500    232     215     125     1200    161     67      0      
  1200   161     67      0       1700    130     30      30
  1700   130     30      30      2800    110     110     110
  2800   110     110     110     4000    255     255     255
  4000   255     255     255     6000    255     255     255

Lalu simpan dengan nama "elev.cpt".

Setelah itu tulislah baris-baris perintah seperti berikut :

set F=maps.ps
psbasemap -JM25 -R85/150/-18/14 -K -Ba10g10WSne:: > %F%
grdimage indo.nc -R -JM -Celev.cpt -K -O -Ba10g10wsne:: >> %F%
pscoast -JM -R -K -O -Dh -Na/thickest,brown -Ir/0.1p,blue -Lf88/-15.2/17/500+l+jr --LABEL_FONT_SIZE=13p -Tf145/10/0.75i/2 --HEADER_FONT_SIZE=13p --HEADER_OFFSET=0.05i >>%F%
psxy -R -JM -W5 -Sf0.4i/0.05ilt -Gblack -m -O -K trench.gmt>> %F%
psxy -R -JM -Wthicker -O -K transform.gmt -m >> %F%
psxy -R -JM -Wthinner,- -m -O -K GSFML_SF_FZ_KM.gmt >> %F%
psxy -R -JM -Wthinnest,- -m -O GSFML_SF_ER_KM.gmt >> %F%
ps2raster maps.ps -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Ftectonic

Simpan dengan ekstensi ".bat". Setelah dieksekusi, hasilnya akan seperti berikut :

Pada baris perintah diatas, saya menggunakan "set" untuk menyingkat script. Dan selanjutnya anda harus memberikan tanda '%' diantara variabel yang telah ditetapkan pada "set". Pada perintah pscoast, saya tambahkan :

"-Na/thickest,brown" = untuk menampilkan batas negara dengan warna coklat.

"-Ir/0.1p,blue" = untuk menampilkan aliran sungai dengan ukuran 0.1p dan berwarna biru.

"-Lf88/-15.2/17/500+l+jr" = untuk menampilkan skala peta pada koordinat 88 BT/-15.2 LS dengan panjang 500 km dan diberi label dengan perataan ke kanan
"-Tf145/10/0.75i/2" = untuk menampilkan arah mata angin pada koordinat 145 BT/10 LU dengan ukuran 0.75 inch.

"--LABEL_FONT_SIZE=13p , --HEADER_FONT_SIZE=13p dan --HEADER_OFFSET=0.05i" untuk mengatur ukuran font.

Untuk ps2raster ada baiknya anda membaca postingan saya sebelumnya disini.

Selamat mencoba dan semoga bermanfaat. Salam Orang Indonesia...