*     *************************************************************************

      subroutine isetcd( n, ind, x, alpha )

*     *************************************************************************

*     Purpose :
*     ---------

*     This routine sets the vector X to the integer constant
*     ALPHA.
*     The components of the vector X that are considered
*     are the components of indices IND(i) for i=1 up to
*     i=N.

*     Parameters :
*     ------------

*     n       ( int )

*       input  : number of components of X to set to ALPHA.
*                It is also the length of array IND.
*       output : unmodified.

*     ind     ( int )

*       input  : array of length N, containing the indices
*                of the components of X to set to ALPHA.
*       output : unmodified.

*     x       ( dble )

*       input  : a vector containing the componenets to
*                set to ALPHA.
*       output : identical to the input value, except for
*                the components of indices IND(i) for i=1
*                up to i=N, that are set to the constant
*                ALPHA.

*     alpha   ( dble )

*       input  : the integer value to assign to the selected
*                components of X.
*       output : unmodified.

*     Routines used :
*     ---------------

*     None.

*     Programming :
*     -------------

*     D. Tuyttens

*     ========================================================================

*     Routine parameters

      integer  n, ind(*), x(*), alpha

*     Internal variables

      integer  i

      do 10 i = 1 , n
        x(ind(i)) = alpha
 10   continue
*
      return
      end