program output
        parameter (lons=721,lats=361,mons=12)
        dimension dd(lons,lats,mons)

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



        do k=1,12
        read(8) ((dd(i,j,k), i=1,lons), j=1,lats)
        enddo

        stop
        end