Esope source files

   1 : C MTOGY2    SOURCE    DC        05/02/15    21:15:50     5030           
   2 :       SUBROUTINE MTOGY2(LRE,NDDL,REL,VROT,RELB)                                
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C  Calcul de la matrice du couplage de Coriolis à partir de la matrice de masse   
   6 : C  pour les éléments massifs en Fourier 
   7 : C
   8 : C     Didier COMBESCURE Janvier 2005                                             
   9 : C                                                                               
  10 : C=======================================================================        
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  13 :       DIMENSION REL(LRE,1),RELB(2*LRE,1),VROT(*),ROME(6,6)
  14 : C
  15 :       NNOEU = LRE/NDDL
  16 : C
  17 : C
  18 :       DO 10 I=1,6
  19 :        DO 10 J=1,6
  20 :         ROME(I,J)=0.D0
  21 :   10  CONTINUE
  22 : C
  23 : C Signe choisi pour etre coherent avec Genta....
  24 : C
  25 :       ROME(1,5) =VROT(2)
  26 :       ROME(5,1) =VROT(2)
  27 :       ROME(2,6) = (-1.)*VROT(2)
  28 :       ROME(6,2) = (-1.)*VROT(2)
  29 : C
  30 : C                                 
  31 :       DO   100 IN=1,NNOEU
  32 : C
  33 :         DO   100 IM=1,NNOEU
  34 : C
  35 :           DO   100 IDD1=1,NDDL
  36 : C
  37 :             DO   100 IDD2=1,NDDL
  38 : C
  39 :             RELB((6*IN)-7+2*IDD1,(6*IM)-7+2*IDD2)=
  40 :      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2-1) 
  41 :             RELB((6*IN)-6+2*IDD1,(6*IM)-6+2*IDD2)=
  42 :      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2) 
  43 :             RELB((6*IN)-6+2*IDD1,(6*IM)-7+2*IDD2)=
  44 :      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1,2*IDD2-1) 
  45 :             RELB((6*IN)-7+2*IDD1,(6*IM)-6+2*IDD2)=
  46 :      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(2*IDD1-1,2*IDD2) 
  47 : C 
  48 :   100  CONTINUE
  49 : C                                   
  50 :       RETURN                                                                    
  51 :       END    
  52 :                                         
  53 :  
  54 :