c    TRMM 3G68 hourly rainfall rate

c    coverage:
c        39.75S  --> 39.75N;
c        179.75W --> 179.75E
c       8 Dec 1997 --> 28 Feb 2006
c    resolution:
c        0.5 deg lat x 0.5 deg lon
c    missing values: -9999.0

c     example program

      parameter (lons=720,lats=160, khr=24)
      real*4  tmi(lons,lats),pr(lons,lats),comb(lons,lats)

      open(21,file='2006/2006-02-28.out',
     &       form='unformatted',access='sequential')

      do 180  k=1,khr
      read(21) ((tmi(i,j), i=1,lons),j=1,lats)
      read(21) ((pr(i,j),  i=1,lons),j=1,lats)
      read(21) ((comb(i,j),i=1,lons),j=1,lats)

c      tmi(i,j):  TMI mean rain rate (mm/hour)
c      pr(i,j):   PR mean rain rate  (mm/hour)
c      comb(i,j): PR and TMI combined rain rate (mm/hour)

      write(6,*) tmi(240,81)
      write(6,*) pr(240,81)
      write(6,*) comb(240,81)
 180  continue

       stop
       end

c       Contact:
c       Yongsheng Zhang: yshzhang@hawaii.edu
c       March 20, 2006