!------------------------------ ! dot prodict of two vectors ! of length 3 !------------------------------ subroutine dp(x,y,d) implicit none real :: d real, dimension(3) :: x, y d = sum(x*y) end subroutine dp