c      
c       This dataset is from NCAR achived 
c       ECMWF TOGA Global Supplementary Fields
c       http://dss.ucar.edu/datasets/ds111.3/

        program output
        parameter (lons=320,lats=160,kvar=6)
        real*4 data(lons,lats,kvar)

        open(8,file='data.file',
     &  form='unformatted',access='sequential')



        k=365     ! or 366
        do 110 idy=1,kdy*2
        do 120 ivar=1,kvar

        read(8) ((data(i,j,ivar),i=1,lons),j=1,lats)

c       for ivar= 
c                1 Surface latent heat flux   [Wm**-2 s]
C                2 Surface thermal radiation  [Wm**-2 s]
c                3 Surface solar radiation    [Wm**-2 s]
c                4 Surface sensible heat flux [Wm**-2 s]
c                5 East/West surface stress   [Nm**-2 s]
c                6 North/South surface stress [Nm**-2 s]

  120   continue
  110   continue
 
        stop
        end