c    TRMM 3B42 V6 three hourly merged rainfall
c    ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B42_V6


c    1) description of the data
c        rr(i,j):  Hourly Rain Rate (mm/hr)
c        er(i,j):  relative error (mm/hr)
c    2) coverage:
c        179.875S => 179.875E
c        49.875S  => 49.875N
c      January 1,1998 to June 30, 2006
c    3) resolution:
c        0.25 deg lat x 0.25 deg lon
c    4) missing values: -9999.0
c    5) example program: read one day data from 00z,03z,....21z.

      parameter (lons=1440,lats=400, ihr=8) 
      real*4  rr(lons,lats),er(lons,lats)

      open(21,file='1998/3B42.980130.dat',
     &  form='unformatted',access='sequential')

      do 180  k=1,ihr
      read(21) ((rr(i,j),i=1,lons),j=1,lats)
      read(21) ((er(i,j),i=1,lons),j=1,lats)

      write(6,*) rr(720,200)
      write(6,*) er(720,200)
 180  continue

       stop
       end

c       Contact:
c       Yongsheng Zhang: yshzhang@hawaii.edu
c       Augut 12, 2006