C DEVLK0    SOURCE    CHAT      05/01/12    22:46:43     5004           
C DEVLK0    SOURCE    PLAF      89/08/01    21:13:53                            
      SUBROUTINE DEVLK0(Q1,XK,XM,FTOTA,NA1,NPAS,IND)                            
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
*--------------------------------------------------------------------*          
*                                                                    *          
*     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *          
*     ________________________________________________               *          
*                                                                    *          
*     Actualisation des forces ext{rieures totalis{es sur la base A  *          
*     pour le premier demi pas de temps.                             *          
*                                                                    *          
*     Param}tres:                                                    *          
*                                                                    *          
* e   Q1(,)   D{placements g{n{ralis{s                               *          
* e   XK      Matrice de raideur g{n{ralis{e                         *          
* e   XM      Matrice de masse g{n{ralis{e                           *          
* es  FTOTA   Forces ext{rieures totalis{es, sur la base A           *          
* e   NA1     Nombre total d'inconnues dans la base A                *          
* e   NPAS    Num{ro du pas de calcul courant                        *          
* e   IND     Indice du pas                                          *          
*                                                                    *          
*     Auteur, date de cr{ation:                                      *          
*                                                                    *          
*     Denis ROBERT-MOUGIN, le 22 mai 1989.                           *          
*                                                                    *          
*                                                                    *          
*--------------------------------------------------------------------*          
      REAL*8 Q1(NA1,*),XK(*),XM(*),FTOTA(NA1,*)                                 
*                                                                               
*     Boucle sur les inconnues de la base A                                     
*                                                                               
      DO 100 I=1,NA1                                                            
         FTOTA(I,IND) = FTOTA(I,IND) - XK(I)*Q1(I,IND)                          
 100     CONTINUE                                                               
*     end do                                                                    
*                                                                               
      END