program read_geos4_3hrly

c  This fortran program reads in the 3-hour mean geos4 
c  surface temperature (K), and short wave radiation reaching the
c  ground (watts/m^2).

      real (kind=4) prec(288,181),t2m(288,181),sw(288,181)

      open(unit=10,file=
     $'geos4.surface.3hr.1998'
     $,form='UNFORMATTED',access='DIRECT',recl=4*288*181 )

      irec =1
      do nday = 1,365
      do nhr =1,8 
       
       read(10,rec=irec) t2m
        irec=irec+1
       read(10,rec=irec) sw
        irec=irec+1

      enddo
      print *, nday
      enddo
     
      stop
      end