program read_geos4_3hrly

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

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

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

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

      print *, nday
      enddo
     
      stop
      end