program read

	implicit none

	integer :: i,j,k,l

	real,dimension(150,150,15,13) :: in

	open(1,file='radar_grid_all-steps.unf',form='unformatted')
	open(2,file='radar_grid_all-steps.bin',form='binary',status='unknown')

!	Ferret unformatted dump yields X number of values in each record where X is length of x-axis.  
!	There will be Y,Z and T number of records depending on the dimensions of the array.

	do l=1,13
	do k=1,15
	do j=1,150

	read(1) (in(i,j,k,l), i=1,150)

	enddo
	enddo
	enddo

	write(2)in

	end program read