writing x,y,z,t grids using fortran

[Joy Romanski]

Hi GrADS users,

I need to write (x,y,z,t) grids using Fortran 77 such that GrADS can
read them.  When I use the code below,

       open(15,file='myfile',
     *   form='unformatted',access='direct',action='write',
     *   recl=4*nlon*nlat*nlev)

       do it=1,nt
         write(15,rec=it) (((data(ilon,ilat,ilev,it),ilon=1,nlon),
     *     ilat=1,nlat),ilev=1,nlev)
       enddo

the first grid (x, y varying, z=1,t=1) is correct but all of the
other grids are wrong.  What am I doing wrong?

Thanks!

Joy Romanski