Data GSMaP yang sudah didownload masih dalam format *.dat.gz, contoh file tanggal 01 Januari 2020 : gsmap_nrt.20200101.0.1d.daily.00Z-23Z.dat.gz. Untuk itu file tersebut perlu diekstrak terlebih dahulu dengan menggunakan unzip sehingga nama filenya menjadi gsmap_nrt.20200101.0.1d.daily.00Z-23Z.dat.
Untuk plotting data diatas, dibutuhkan file deskriptor yang lazim juga disebut dengan kontrol file sehingga ekstensinya adalah *.ctl sebagai berikut :
DSET ^gsmap_nrt.%y4%m2%d2.0.1d.daily.00Z-23Z.dat
TITLE GSMaP_NRT 0.1deg Daily (00:00Z-23:59Z)
UNDEF -999.9
OPTIONS YREV LITTLE_ENDIAN TEMPLATE
XDEF 3600 LINEAR 0.05 0.1
YDEF 1200 LINEAR -59.95 0.1
zdef 1 levels 1013
tdef 365 linear 00:00z1jan2020 1dy
VARS 1
precip 0 99 daily averaged precip(mm/hr)
ENDVARS
untuk penjelasan dari tiap item dalam file deskriptor (*.ctl) diatas dapat dibaca : disini
Copy semua baris file *.ctl diatas dan Paste di Notepad, kemudian simpan dengan nama gsmap.ctl pada folder yang sama dengan folder penyimpanan data.
Selanjutnya copy dan paste baris perintah GrADS berikut ke notepad dan beri nama filenya dengan nama plotgsmap.gs :
'reinit'
Untuk plotting data diatas, dibutuhkan file deskriptor yang lazim juga disebut dengan kontrol file sehingga ekstensinya adalah *.ctl sebagai berikut :
DSET ^gsmap_nrt.%y4%m2%d2.0.1d.daily.00Z-23Z.dat
TITLE GSMaP_NRT 0.1deg Daily (00:00Z-23:59Z)
UNDEF -999.9
OPTIONS YREV LITTLE_ENDIAN TEMPLATE
XDEF 3600 LINEAR 0.05 0.1
YDEF 1200 LINEAR -59.95 0.1
zdef 1 levels 1013
tdef 365 linear 00:00z1jan2020 1dy
VARS 1
precip 0 99 daily averaged precip(mm/hr)
ENDVARS
untuk penjelasan dari tiap item dalam file deskriptor (*.ctl) diatas dapat dibaca : disini
Copy semua baris file *.ctl diatas dan Paste di Notepad, kemudian simpan dengan nama gsmap.ctl pada folder yang sama dengan folder penyimpanan data.
Selanjutnya copy dan paste baris perintah GrADS berikut ke notepad dan beri nama filenya dengan nama plotgsmap.gs :
'reinit'
tanggal = 01 Januari 2020
'open PATH/gsmapctl.ctl'
'set t 1'
'set lon 95 141'
'set lat -15 15'
'set mpdset hires'
'set csmooth on'
'set gxout shaded'
'set clevs 0 5 20 50 100 150'
'd precip*24'
'cbarn'
'draw title Distribusi Curah Hujan GSMaP 01 Januari 2020'
'printim PATH/01/01Januari2020.png white'
'set gxout print'
'set prnopts %6.2f 1 1'
write('PATH/'tanggal'.txt', 'X Y LON LAT Precip')
'q dims'
xline=sublin(result,2) ;* 2nd line
yline=sublin(result,3) ;* 3rd line
xmax=subwrd(xline,13) ;*13th word on xline
ymax=subwrd(yline,13) ;*13th word on yline
say 'X grid-points: 'xmax
say 'Y grid-points: 'ymax
y=570.5
ymax=595.5
while(y<=ymax)
x=1006.5
xmax=1045.5
while(x<=xmax)
'set x 'x
'set y 'y
'd precip*24'
* NOTE: It may be useful to test this to find out where the data is contained with in the result
* It just so happens that in this case, the data is the 1st word of the 2nd line, this is not always true
precip=sublin(result,2)
precip=subwrd(precip,1)
* Get Lat/Lon Data
'q dims'
lons=sublin(result,2)
lats=sublin(result,3)
lon=subwrd(lons,6)
lat=subwrd(lats,6)
* Save data to file
* Note the "append", so to add to the file instead of overwriting it
write('PATH/'tanggal'.txt', x' 'y' 'lon' 'lat' 'precip,append)
x=x+1
endwhile
y=y+1
endwhile
####
Keterangan : PATH adalah alamat folder target penyimpanan file data dan file .ctl serta file .gs yang dibatasi oleh garis miring.
untuk menjalankan script GrADS diatas, masuk dulu ke aplikasi GrADS dengan klik dua kali pada icon GrADS. Maka akan tampil sebagai berikut :
Tekan ENTER maka akan tampil
selanjutnya pada GrADS prompt, ketik run PATH/plotgsmap.gs maka akan tampil plot curah hujan tanggal 01 Januari 2020 wilayah Indonesia sekaligus menyimpan data tiap grid dalam format .txt kedalam folder target.
Anda bisa juga melakukan plotting GSMaP yang telah tersimpan dalam format .txt dengan menggunakan aplikasi ArcGIS... https://arifmarufi.blogspot.com/2020/01/cara-plot-data-curah-hujan-harian-gsmap.html
Mau tanya, GSMaP kan datanya adanya harian, jika mau memetakan curah hujan bulanan pakai GSMaP sebaiknya bagaimana? Terima kasih
BalasHapusAda 2 cara jika ingin mendapatkan data bulanan :
Hapus1. Bisa dengan menjumlahkan data harian dalam bulan yang sama
2. Bisa juga dengan melalui alamat : ftp://hokusai.eorc.jaxa.jp/realtime/monthly/
salam
semoga berhasil
Pak bagaimana ya jika saat di plot ke grads, hasilnya all undefined value?
BalasHapusoh harusnya bisa.. mungkin bisa lebih detil bisa didiskusikan via email saya marufi.arif@yahoo.com
Hapussama saya juga seperti ini pak, izin email pak
HapusIjin bertanya Pak, kenapa saat display itu precip*24? kenapa tidak d precip saja? apakah pengaruhnya karena plotting akumulasi curah hujan harian ya Pak? terima kasih Pak.
BalasHapusitu saya kalikan 24 karena di metadata gsmap nya ditulis bahwa datanya average daily yang menurut saya adalah rata-rata harian. sehingga ketika kita akan melihat data total harian harus dikalikan 24 jam...
Hapuspernah saya coba langsung d precip saja tapi nilainya sangat kecil sekali... terima kasih diskusinya... coba cek lagi metadata gsmap nya... untuk memastikan salam..
bisa dibaca di file README.first.txt di alamat : ftp://hokusai.eorc.jaxa.jp/realtime/
HapusSore mas Arif.
BalasHapussaya jalankan plotgsmap nya, ada error
"Unable to locate ENDWHILE statement for the WHILE statement at line 32"
knapa ya??
perlu diperhatikan identasi dari WHILE dan ENDWHILE nya, identasi adalah posisi pengetikan yang menjorok kedalam.. jika salah biasanya akan error...
Hapuskalau masih error, boleh tinggalkan alamat emailnya nanti saya kirimkan file gs yang sudah benar identasinya... salam
HapusTerima Kasih pak Arif Ma'rufi.
HapusSetelah mentoring dengan bapak kurang lebih 2 jam secara online, tugas saya jadi lebih cepat. Sekaligus mendapatkan ilmu baru, terutama untuk alternatif mendapatkan data hujan di Indonesia.
Salam hormat dan tetap jaga kesehatan.
sama-sama pak, semoga bermanfaat. dan semoga kita semua selalu dilindungi Allah SWT.
HapusSalam Sehat
kalau saya mau filter daerah daratn aja gimana ya caranya?
BalasHapusada dua alternatif :
Hapus1. Menggunakan perintah wrseries, baca : https://arifmarufi.blogspot.com/2019/11/cara-download-dan-ekstrak-data-chirps.html
2. Hasil ekstrak diatas dibuka di ArcGIS dan dilakukan proses clip pada geoprocessing wizard
semoga bermanfaat
salam
Izin bertanya pak, ketika saya mencoba menjalankan script plotgsmap di atas, pada tampilan grads tertera 'unknown command : 1' dan setelah saya cek file txt nya ternyata kosong pak. Ini bagaimana pak? Mohon arahannya pak
BalasHapuskemungkinan domain wilayahnya salah sehingga tidak berhasil diekstrak...
Hapuskalau boleh tau, itu mau ekstrak wilayah mana ?
Izin pak, saya rencananya akan mengekstrak wilayah Kabupaten Lebak, banten pak
Hapusmungkin bisa dijelaskan di email saja mbak...
Hapusemail saya marufi.arif@yahoo.com
salam
siap pak, terima kasih
HapusMohon maaf pak, untuk download GrADS nya dimana ya? saya mencoba download di internet tidak ada
BalasHapushttp://opengrads.org/ coba disini mbak
HapusIzin bertanya Kak, kalau untuk per jamnya curah hujannya
BalasHapusscritnya seperti apa Kak?
prinsipnya sama saja sih.. tinggal download file nya yang jam2an.. dan di bagian d precip*24 diubah d precip saja
HapusTerima kasih mas Arif Makrufi.. blognya sangat bermanfaat.. terima kasih sudah memberi inspirasi..
BalasHapusterima kasih supportnya pak, semoga bermanfaat..
Hapussalam sehat
Terima kasih tutorialnya. Jika ingin menghitung curah hujan selama 10 hari bagaimana? Apakah cukup dengan d ave(precip,t=1,t=10)*24?
BalasHapusKalau total 10 hari seharusnya menggunakan fungsi sum (jumlah) bukan ave (rata2)
Hapusmohon izin bertanya pak, kok saya ada tulisan Open Error: Missing or invalid dimension starting value
BalasHapus--> The invalid description file record is:
--> TDEF 365 LINEAR 00:00Z9okt2022 1dy
mksudnya gmana ya pak?
TDEF 365 LINEAR 00:00Z9okt2022 1dy artinya ada 365 t (time), dimulai dari tanggal 09 oktober 2022.. tp menurut saya yang salah adalah penulisan singkatan bulan.. bukan okt tapi oct..
Hapusjadi penulisan singkatan bulan harus dalam bahasa inggris
semoga bermanfaat
Komentar ini telah dihapus oleh pengarang.
BalasHapusPak, untuk nilai
BalasHapusy=570.5
ymax=595.5
while(y<=ymax)
x=1006.5
xmax=1045.5
while(x<=xmax)
didapatkan dari mana ya? apakah area? klo saya cuma mau plot pulau sumatra saya, brp nilai yang harus sy inputkan?
Terima kasih
sepertinya itu hanya untuk area Jambi. itu adalah nomor urut grid X dan Y. Memang membatasi wilayah atau area. untuk menghitung suatu wilayah/area kita lihat dulu file ctl nya, disitu terlihat bahwa XDEF 0.05 0.1 yang berarti start di 0.05 derajat dan ukuran grid nya 0.1 derajat sehingga total gridnya menjadi 3600 grid (karena bujur data yang tersedia 360 derajat). sedangkan YDEF -59.95 0.1 yang berarti start di -59.95 derajat dengan ukuran grid 0.1 derajat sehingga total grid nya menjadi 1200 (karena total lintang yang tersedia dalam data tsb adalah 120 derajat).
Hapusmisalnya kita ingin mengambil wilayah dengan koordinat 90 derajat sd 100 derajat, maka kita hitung mulai dari 0.05 derajat dengan penambahan 0.1 derajat sampai ketemu 90 derajat, nomor barisnya adalah Xawal dan Xakhir adalah nomor baris pada 100 derajat.
demikian pula cara menghitung nomor grid Y nya.. tp start dari -59.95 + 0.1 sesuai dengan panduan file ctl.
maaf pak, saya juga terkendala untuk mencari nilai x y nya karena masih belum paham
Hapusdibuat di excel untuk menghitung nomor grid x dan nomor grid y nya... diurutkan ke bawah, dimana x start dari 0.05 dengan penjumlahan 0.1 hingga 3600 grid. sedangkan y start -59.5 dengan penjumlahan 0.1 hingga 1200 grid. setelah itu tinggal menyesuaikan koordinat mana yang akan kita ambil.. samakan dengan nomor grid nya...
Hapuskalau masih kesulitan bisa tinggalkan emailnya nanti saya share caranya...
Hapuscara mendapatkan x dan y itu bagaimana ya pak
BalasHapusya dihitung seperti saya jelaskan diatas
Hapussiang pak, izin bertanya saya juga mengolah data gsmap mau diubah dulu ke nc dan mau coba olah pakai python, tapi saat udat dat ke nc saya mengalami kendala, tertulis file notfound bagaimana cara memperbaikinya ya pak?
BalasHapuswah kalau itu perlu diurut dari awal lagi.. saran saya didownload ulang saja file nya
HapusTerima kasih Pak... blog tutorial ini sangat membantu penelitian saya saat harus download data GSMaP.
BalasHapusTerima kasih juga kepada Bapak Arif Ma'rufi atas tutorialnya. semoga Bapak diberikan kelancaran dan selesai tepat waktu dalam menempuh program S3...aamiin.
sama2 pak.. terima kasih doanya...
Hapusterima kasih... blog tutorial ini sangat membantu penelitian saya saat harus download data GSMaP.
BalasHapusterima kasih juga kepada Bapak Arif Ma'rufi atas tutorialnya, semoga dalam menempuh program S3 diberikan kelancaran dan dapat menyelesaikan tepat waktu...aamiin.
sama2 pak Djoko Suryanto, semoga bermanfaat bagi kita semua... terima kasih atas doanya....
HapusIzin pak, saya sudah mengikuti semua langkah termasuk file yang didonwload pun sama namun ketika menjalankan perintah run pada grads tidak berjalan dengan notifikasi 'Error opening script file:'. Solusinya bagaimana ya pak?
BalasHapussdh diemailkan
Hapus