C DEVFXA    SOURCE    BP208322  18/12/20    21:15:31     10048          
*      SUBROUTINE DEVFXA(FEXA,IFEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS,
      SUBROUTINE DEVFXA(FEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS,
     &                  FTEXB,FEXB,NPLB,IDIMB,RIGIDE)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*--------------------------------------------------------------------*
*                                                                    *
*     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *
*     ________________________________________________               *
*                                                                    *
*     Ajout des termes de forces ext{rieures libres.                 *
*                                                                    *
*     Param}tres:                                                    *
*                                                                    *
* e   FEXA    Evolution des forces extérieures, en base A            *
* e   IFEXA   Numéro du mode correspondant au point de chargement    *
*             (supprime le 2018-12-14 par bp)
* es  FTOTA   Forces extérieures libres totalisées, en base A        *
* e   NPFEXA  Nombre de points de chargement                         *
* e   NA1     Nombre total d'inconnues dans la base A                *
* e   NPC1    Nombre de pas de calcul                                *
* e   NPAS    Num{ro du pas courant                                  *
* e   FTEXB   Evolution des forces extérieures sur base B, pour le   *
*             calcul des  moments                                    *
* es  FEXB    Forces extérieures sur base B, pour le calcul des      *
*             moments                                                *                                               *
*                                                                    *
*     Auteur, date de création:                                      *
*                                                                    *
*     Denis ROBERT-MOUGIN, le 22 mai 1989.                           *
*                                                                    *
*--------------------------------------------------------------------*
      REAL*8 FEXA(NPFEXA,NPC1,*),FTOTA(NA1,*)
      REAL*8 FEXB(NPLB,2,*),FTEXB(NPLB,NPC1,2,*)
c       INTEGER IFEXA(*)
*
      LOGICAL RIGIDE
*
*     Boucle sur les points de chargement:
*
      NPA = NPAS + 1
      DO 100 I=1,NPFEXA
         FTOTA(I,1) = FEXA(I,NPA,1)
         FTOTA(I,2) = FEXA(I,NPA,2)
 100     CONTINUE
*     END DO
*
* Dans le cas des corps rigides, boucle sur les chargements Base B
      IF (RIGIDE) THEN
         IF (IDIMB.EQ.6) THEN
             IDIM = 3
         ELSE
             IDIM = 2
         ENDIF
         DO 200 IP=1,NPLB
             DO 210 ID=1,IDIM
                FEXB(IP,1,ID) = FTEXB(IP,NPA,1,ID)
                FEXB(IP,2,ID) = FTEXB(IP,NPA,2,ID)
 210         CONTINUE
 200     CONTINUE
      ENDIF
      RETURN
      END