Esope source files

   1 : C CERMAS    SOURCE    CHAT      05/01/12    21:52:26     5004           
   2 :       SUBROUTINE CERMAS(REL,LRE,ROS,COORD)
   3 : C-----------------------------------------------------------------------
   4 : C  ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE
   5 : C  CERCE
   6 : C     
   7 : C  ENTREE
   8 : C          LRE     DIMENSION DU TABLEAU REL
   9 : C          ROS     RHO * SECT
  10 : C          COORD   COORDONNEES DU NOEUD
  11 : C     
  12 : C  SORTIE
  13 : C          REL     MATRICE DE MASSE ELEMENTAIRE
  14 : C-----------------------------------------------------------------------
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)
  17 : -INC CCREEL
  18 : -INC CCOPTIO
  19 :       DIMENSION REL(LRE,*),COORD(3,*)
  20 : *     
  21 :       R = COORD(1,1)
  22 :       IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN
  23 :          FAC=ROS*R*2*XPI
  24 :       ELSE
  25 :          FAC=ROS*R*XPI
  26 :       ENDIF
  27 : *     
  28 :       DO 1 I=1,LRE
  29 :          DO 2 J=1,LRE
  30 :             REL(I,J)=0.D0
  31 :  2       CONTINUE
  32 :  1    CONTINUE
  33 : *     
  34 :       REL(1,1)=FAC
  35 :       REL(2,2)=FAC
  36 :       IF(IFOUR.EQ.1) THEN
  37 :          REL(3,3)=FAC
  38 :       ENDIF
  39 :       RETURN
  40 :       END
  41 : 
  42 : 
  43 : 
  44 :  
  45 :  
  46 :