integer, parameter :: nx=1152,ny=36,nz=32 real, dimension (nx,ny,nz,3) :: bb real, dimension (nx) :: x real, dimension (ny) :: y real, dimension (nz) :: z real :: t ! ! open and read ! open(1,file='data/proc10/bb.dat',form='unformatted') read(1) bb,x,y,z,t close(1) ! ! check ! print*,'max(bb)=',maxval(bb) print*,'min(bb)=',minval(bb) print*,'z=',z print*,'x=',x print*,'t=',t END