program output parameter (lons=72,lats=73,mons=12) real*4 dd(lons,lats,mons) open(8,file='data.dat', & form='unformatted',access='sequential') do mn=1,12 read(8) ((dd(i,j,mn), i=1,lons), j=1,lats) enddo stop end