1      subroutine ifill(n,val,a,ia)
2C$Id$
3      implicit none
4      integer n, val, a(*), ia, i
5c
6c     initialise integer precision array to scalar value
7c
8      if (ia.eq.1) then
9         do 10 i = 1, n
10            a(i) = val
11 10      continue
12      else
13         do 20 i = 1,(n-1)*ia+1,ia
14            a(i) = val
15 20      continue
16      endif
17c
18      end
19