C NDST      SOURCE    CHAT      05/01/13    01:56:29     5004           
      SUBROUTINE NDST(IGAU,XEL,SHPTOT,SHP,BGENE,DJAC)                           
C                                                                               
C ====================================================================          
C     FORME LA MATRICE N POUR L'ELEMENT DST                                     
C     UTILISEE POUR LA MATRICE DE MASSE                                         
C ====================================================================          
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
      PARAMETER (LRE=18,NST=6,NBNO=3)                                           
      DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)                 
      DATA  ZER/0.D0/                                                           
C                                                                               
      CALL ZERO(BGENE,NST,LRE)                                                  
      DO 101 NP=1,NBNO                                                          
      SHP(1,NP)=SHPTOT(1,NP,IGAU)                                               
      SHP(2,NP)=SHPTOT(2,NP,IGAU)                                               
      SHP(3,NP)=SHPTOT(3,NP,IGAU)                                               
101   CONTINUE                                                                  
      CALL JACOBI(XEL,SHP,2,NBNO,DJAC)                                          
C                                                                               
C      MATRICE  - N -                                                           
C                                                                               
      K=0                                                                       
      DO 102 NP=1,NBNO                                                          
      DO 103 I =1,NST                                                           
      BGENE(I,K+I)=SHP(1,NP)                                                    
103   CONTINUE                                                                  
      K=K+6                                                                     
102   CONTINUE                                                                  
      RETURN                                                                    
      END