Cara Mengekstrak Data TRMM

Data TRMM yang sudah didownload dalam format NC (NetCDF) yang merupakan salah satu format standar data kebumian. Untuk membuka, mengekstrak, membuat display dan menganalisa data  format NC dapat menggunakan aplikasi GrADS (Grid Analysis and Display System). Aplikasi GrADS ini merupakan program Open Source yang dapat didownload gratis di  http://cola.gmu.edu/grads/downloads.php. 

Berikut Script GrADS yang digunakan dalam mengekstrak data TRMM Harian menjadi format .txt dan untuk ruang lingkup Provinsi Jambi :

'reinit'

#Misalnya data yang dioalah adalah data TRMM tanggal 25 Juli 2017
tanggal = 25_Juli_2017

#Membuka file TRMM sesuai target folder tempat data tersebut disimpan, melakukan setting lintang bujur sesuai keinginan (Provinsi Jambi) dan menginisiasi parameter.

'sdfopen PATH\3B42RT_Daily.20170725.7.nc4.nc'
'set dfile'
'set lon 101 105'
'set lat -3.0 -0.5'
'set t 1'
'prec1=precipitation'


'set gxout shaded'
#'set clevs 0 20 50 100 150 200 300 400 500'
'd prec1'
'cbarn'
'set display color white'
'draw title Analisis Curah Hujan Prop.Jambi Tgl. 'tanggal
'printim PATH\TRMM_'tanggal'.gif gif x800 y600 white'

'set gxout print'
'set prnopts %6.2f 1 1'


write('PATH\25072017.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=228
ymax=240
while(y<=ymax)

  x=1121
  xmax=1141
  while(x<=xmax)

    'set x 'x
    'set y 'y
    'd prec1'

*    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

     precipitation=sublin(result,2)
     precipitation=subwrd(precipitation,1)

*    Menentukan Lat/Lon Data

     'q dims'
     lons=sublin(result,2)
     lats=sublin(result,3)
     lon=subwrd(lons,6)
     lat=subwrd(lats,6)

*    Menyimpan Data dalam format .txt
*    Note the "append", so to add to the file instead of overwriting it


write('PATH\25072017.txt', x'    'y'    'lon'    'lat'    'precipitation,append)

     x=x+1
   endwhile
  y=y+1
endwhile

***