!
!  set_a_0.f90
!
      subroutine set_a_0(a,n)
      integer :: n
      real :: a(n)
!      
      do i=1,n
        a(i) = 0.
      enddo
!
      end subroutine set_a_0      