program read_geos4_3hrly

c  This fortran program reads in the daily accumulated total  
c  precipitation (mm/day)

      real (kind=4) preacc(288,181)

      open(unit=10,file=
     $'preacc.daily.1998'
     $,form='UNFORMATTED',access='DIRECT',recl=4*288*181 )

      irec =1
      do nday = 1,365
       
       read(10,rec=irec) preacc
        irec=irec+1

      print *, nday
      enddo
     
      stop
      end