PROGRAM MultiArrays

! Written by Kate T-C
! September 14, 2009

implicit none

integer :: i,j,k
integer :: lon, lat, time
integer, dimension(6) :: test1 = (/(i,i=1,6)/)
integer, dimension(3,2) :: test2
integer, dimension(24) :: test3d1 = (/(i,i=1,24)/)
integer, dimension(3,2,4) :: test3d


print *, 'One dimensional array:'
print *, test1

print *, 'Two dimensional array:'
test2 = reshape(test1, (/3,2/))
do i=1,3
	print *, (test2(i,j),j=1,2)
enddo

test3d = reshape(test3d1, (/3,2,4/))
print *, 'Three dimensional array straight up:'
do i=1,3
	do j=1,2
		write *, (test3d(i,j,k),k=1,4)
	enddo
	write *, aChar(10)
enddo

print *, '3d array with time on the outside:'
do time=1,4
	do lon=1,2
		write *, (test3d(lat,lon,time),lat=1,3)
	enddo
	write *, aChar(10)
enddo

END