3
integer :: ndim1, ndim2, ndim3
4
real, dimension( :, :, : ), allocatable :: field
5
integer :: idim1, idim2, idim3
7
open( 5, file = 'filedim', status = 'unknown' )
8
read( 5, * ) ndim1, ndim2, ndim3
11
allocate( field( ndim1, ndim2, ndim3 ))
14
field( 1 : ndim1, 1, 1 ) = 1.
19
print*, field( idim1, idim2, idim3 )