Disini saya akan menjelaskan cara menampilkan Wind Barbs dengan GMT 5. Wind Barbs adalah diagram yang menunjukkan arah dan kecepatan angin. Untuk menampilkan diagram ini, kita harus membuat file definisi simbol Wind Barbs. File ini saya dapat dari sebuah forum di website GMT dan saya modifikasi agar mampu menampilkan data hingga kecepatan lebih dari 200 knot. File definisi tersebut dapat didownload disini. Setelah didownload, simpan di C:\programs\gmt5\share\custom. Download contoh data angin disini dan file cpt disini.
Selanjutnya buatlah script berikut :
>> set ps=windbarb.ps
>>
>> pscoast -JM6i -R90/110/-13/10 -Gwheat -Sazure1 -Dh -Wthin,black -Bafg -BWSen -K > %ps%
>> gawk "{if($3>=0 && $3 <=50) print $1,$2,$3,$4}" angin.dat | psxy -J -R -W0.7p,blue -Skwindbarb/0.7i -O -K -Gblue >> %ps%
>> gawk "{if($3>50 && $3 <=100) print $1,$2,$3,$4}" angin.dat | psxy -J -R -W0.7p,green -Skwindbarb/0.7i -O -K -Ggreen >> %ps%
>> gawk "{if($3>100 && $3 <=150) print $1,$2,$3,$4}" angin.dat | psxy -J -R -W0.7p,orange -Skwindbarb/0.7i -O -K -Gorange >> %ps%
>> gawk "{if($3>150 && $3 <=200) print $1,$2,$3,$4}" angin.dat | psxy -J -R -W0.7p,red -Skwindbarb/0.7i -O -K -Gred >> %ps%
>> gawk "{if($3>200) print $1,$2,$3,$4}" angin.dat | psxy -J -R -W0.7p,purple -Skwindbarb/0.7i -O -K -Gpurple >> %ps%
>>
>> psscale -Ct.cpt -Dx8c/2.5c+w12c/0.5c+jTC+h -Bxaf+l"knot" -O >> %ps%
>> psconvert -A -Tg -P %ps% -F%ps% -E300
Simpan dengan ekstensi .bat dan jalankan. Simbol pada Wind Barbs disini memerlukan 4 kolom. Kolom tersebut adalah bujur, lintang, kecepatan dalam knot, dan arah dalam derajat.
Selamat mencoba.
0 komentar:
Post a Comment