!==========================================
! dotprod.f90
! prompt for two real vectors of length 3
! calculate dot product using do loop
! print result
!==========================================
program dotprod
  implicit none
  integer :: i
  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 = 0.0
  do i = 1, 3
     c = c + a(i)*b(i)
  enddo

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

end program dotprod