real function dprodc(a1,a2)

      implicit none
      real,dimension(3)::a1,a2

      dprodc=a1(1)*a2(1)+a1(2)*a2(2)+a1(3)*a2(3)

      end function dprodc