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