             subroutine addvec(a,b,c,D) 
             dimension a(1),b(1),c(1) 
             do i = 1,D 
               c(i) = a(i) + b(i) 
             enddo
             return 
             end  
