Hi,
For the following code, some complilers write out 1.0 and some write out
0.0 for op(1:3,1:3,1). Which one is correct, or standard conforming?
First time to post a question. Thanks.
S Yuan
program main
implicit none
real :: op(3,4,2), s
op=0.0
s=1.0
call adds(op(1:3,1:3,1),op(1:3,1:3,1),s)
write(*,"(3(3f8.4,/))")op(1:3,1:3,1)
stop
contains
subroutine adds(b,a,s)
real, intent(in out) :: b(3,3),a(3,3),s
real :: c(3,3)
c=a+s
b=c
! write(*,"(3(3f8.4,/))")b
return
end subroutine adds
end program main
|