C MTOKCE    SOURCE    BP208322  12/03/08    21:15:04     7294
c       SUBROUTINE MTOKCE(LRE,NDDL,REL,VROT,RELB)
      SUBROUTINE MTOKCE(LRE,NDDL,REL,ROME,RELB)
C=======================================================================
C
C  Calcul de la matrice de raideur centrifuge à partir dela matrice de masse
C  pour les éléments massifs en 3D
C
C     Didier COMBESCURE Mars 2003
C     BP on sort le remplissage de ROME (cst) de la boucle (cf KCENT3.eso)
C
C=======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
c       DIMENSION REL(LRE,1),RELB(LRE,1),VROT(*),ROME(3,3)
      DIMENSION REL(LRE,1),RELB(LRE,1),ROME(3,3)
C
      NNOEU = LRE/NDDL
C
c       ROME(1,1) = (-1.)*((VROT(2)**2) + (VROT(3)**2))
c       ROME(2,2) = (-1.)*((VROT(1)**2) + (VROT(3)**2))
c       ROME(3,3) = (-1.)*((VROT(1)**2) + (VROT(2)**2))
c       ROME(1,2) = VROT(1)*VROT(2)
c       ROME(1,3) = VROT(1)*VROT(3)
c       ROME(2,3) = VROT(2)*VROT(3)
c       ROME(2,1) = ROME(1,2)
c       ROME(3,1) = ROME(1,3)
c       ROME(3,2) = ROME(2,3)
C
C
      DO   100 IN=1,NNOEU
C
        DO   100 IM=1,NNOEU
C
          DO   100 IDD1=1,NDDL
C
            DO   100 IDD2=1,NDDL
C
            RELB((3*IN)-3+IDD1,(3*IM)-3+IDD2)=
     . REL((3*IN)-3+IDD1,(3*IM)-3+IDD1)*ROME(IDD1,IDD2)
C
  100      CONTINUE
C
      RETURN
      END