Esope source files

   1 : C NDST      SOURCE    CHAT      05/01/13    01:56:29     5004           
   2 :       SUBROUTINE NDST(IGAU,XEL,SHPTOT,SHP,BGENE,DJAC)                           
   3 : C                                                                               
   4 : C ====================================================================          
   5 : C     FORME LA MATRICE N POUR L'ELEMENT DST                                     
   6 : C     UTILISEE POUR LA MATRICE DE MASSE                                         
   7 : C ====================================================================          
   8 :       IMPLICIT INTEGER(I-N)
   9 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  10 :       PARAMETER (LRE=18,NST=6,NBNO=3)                                           
  11 :       DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)                 
  12 :       DATA  ZER/0.D0/                                                           
  13 : C                                                                               
  13 :       CALL ZERO(BGENE,NST,LRE)                                                  
  14 :       DO 101 NP=1,NBNO                                                          
  15 :       SHP(1,NP)=SHPTOT(1,NP,IGAU)                                               
  16 :       SHP(2,NP)=SHPTOT(2,NP,IGAU)                                               
  17 :       SHP(3,NP)=SHPTOT(3,NP,IGAU)                                               
  18 : 101   CONTINUE                                                                  
  18 :       CALL JACOBI(XEL,SHP,2,NBNO,DJAC)                                          
  19 : C                                                                               
  20 : C      MATRICE  - N -                                                           
  21 : C                                                                               
  22 :       K=0                                                                       
  23 :       DO 102 NP=1,NBNO                                                          
  24 :       DO 103 I =1,NST                                                           
  25 :       BGENE(I,K+I)=SHP(1,NP)                                                    
  26 : 103   CONTINUE                                                                  
  27 :       K=K+6                                                                     
  28 : 102   CONTINUE                                                                  
  29 :       RETURN                                                                    
  30 :       END                                                                       
  31 :