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