Sources Esope | Cast3M

ndst
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