2 ! PR 15553 : the array used to be filled with garbage
3 ! this problem disappeared between 2004-05-20 and 2004-09-15
7 double precision x(10,10)
18 if (any(reshape(x(1:10,1:3), (/ 30 /)) &
19 /= (/ 11, 12, 13, 14, 15, 16, &
20 21, 22, 23, 24, 25, 26, &
21 31, 32, 33, 34, 35, 36, &
22 41, 42, 43, 44, 45, 46, &
23 51, 52, 53, 54, 55, 56 /))) STOP 1
27 subroutine pack (arr
, ni
, nj
)
28 integer, intent(in
) :: ni
, nj
29 double precision, intent(inout
) :: arr(:,:)
30 double precision :: tmp(ni
,nj
)
31 tmp(:,:) = arr(1:ni
, 1:nj
)
32 call copy (arr
, tmp
, ni
, nj
)
35 subroutine copy (dst
, src
, ni
, nj
)
36 integer, intent(in
) :: ni
, nj
37 double precision, intent(out
) :: dst(ni
, nj
)
38 double precision, intent(in
) :: src(ni
, nj
)