C DDOTPV    SOURCE    CHAT      05/01/12    22:38:16     5004           
      function ddotpv(lon,val1,val2)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      dimension val1(*),val2(*)
      resu1=0.d0
      resu2=0.d0
      resu3=0.d0
      resu4=0.d0
      if (lon.gt.0) then
      do 10 i=0,lon-4,4
       resu1=resu1+val1(i+1)*val2(i+1)
       resu2=resu2+val1(i+2)*val2(i+2)
       resu3=resu3+val1(i+3)*val2(i+3)
       resu4=resu4+val1(i+4)*val2(i+4)
  10  continue
      do 20 i=4*(lon/4)+1,lon
       resu1=resu1+val1(i)*val2(i)
  20  continue
      endif
      ddotpv=resu1+resu2+resu3+resu4
      return
      end