Esope source files

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