!==========================================
! dotprod.f90
! prompt for two real vectors of length 3
! calculate dot product
! print result
!==========================================
program dotprod
  implicit none
  real :: c
  real, dimension(3) :: a, b

  !-------------------------------
  ! enter data
  !-------------------------------
  print*,'Enter first vector'
  read*, a
  print*,'Enter second vector'
  read*, b

  !-------------------------------
  ! calculate dot product
  !-------------------------------
  c = a(1)*b(1) + a(2)*b(2) + a(3)*b(3)

  !-------------------------------
  ! print result
  !-------------------------------
  print*,'Dot product = ', c

end program dotprod