Membuat Profil Kedalaman dan Ketinggian Dengan GMT

Kali ini saya akan berbagi cara membuat profil kedalaman laut dan ketinggian daratan dengan GMT. Data yang saya pakai disini adalah data GEBCO dengan resolusi 30 detik.

Untuk membuatnya, buatlah script berikut :

set F=tsun.ps
set R=105/115/-12/-5
set azimuth=100
set elevation=20

#Membuat file warna
makecpt -Z -Cglobe > elev.cpt

#Menampilkan peta 3D dan wilayah yang di slice
#Wilayah yang dislice disini adalah 110BT-5LS & 110BT-12LS
#Pada bagian echo 110 -8.5, nilai -8.5 adalah hasil trial and error agar
#garis berada tepat ditengah
grdview indo.nc -Jx1.5c -Jy1.5c -R%R%/-11000/4000 -Celev.cpt -JZ1.5c -Qi  -E%azimuth%/%elevation% -K -B2::/2::/1000::wESnz+ -Y11 -X1.5> %F%
echo 110 -8.5 4000 > sliceh
echo 110 -8.5 2000 >> sliceh
echo 110 -8.5 0 >> sliceh
echo 110 -8.5 -2000 >> sliceh
echo 110 -8.5 -4000 >> sliceh
echo 110 -8.5 -6000 >> sliceh
echo 110 -8.5 -8000 >> sliceh
echo 110 -8.5 -11000 >> sliceh
psxyz sliceh -J -JZ -R -O -K -E%azimuth%/%elevation% -Sy10.5 -W2 >> %F%
echo 110 -12 -4000 > slicev
echo 110 -5 -4000 >> slicev
psxyz slicev -J -JZ -R -O -K -E%azimuth%/%elevation% -Sz1.5 -W2 >> %F%

#Menampilkan peta 2D dan wilayah yang di slice
grdview indo.nc -Jx1c -Jy1c -R%R%/-11000/4000 -Celev.cpt -Qi -K -B2::/2::wESn -X16 -O >> %F%
echo 110 -5 > slice
echo 110 -12 >> slice
psxy slice -J -R -O -K -W6 >> %F%

#Menampilkan hasil slice
#Pada bagian "project", -G1 artinya nilai sampling 1 km
project -C110/-5 -E110/-12 -Q -G1 > track
grdtrack track -Gindo.nc > tracked
gawk "{print $2, $4}" tracked | psxy -JX8i/2.5i -R-12/-5/-8000/4000 -Wthick -O -K -Y-8.5 -X-12 -B2:lintang:/3000:meter:WSne >> %F%
echo -12 0 > nol
echo -5 0 >> nol
psxy nol -J -R -O -Wthick,blue >> %F%

#Print peta
ps2raster %F% -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fslice

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