C POUROG    SOURCE    CHAT      05/01/13    02:19:50     5004           
      SUBROUTINE POUROG(REL,LRE,P,RP)                                           
C-----------------------------------------------------------------------        
C  ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES        
C  GLOBAUX                                                                      
C  POUR UNE MATRICE ANTISYMETRIQUE                                                              
C                                                                               
C  ENTREE                                                                       
C         LRE    NOMBRE DE LIGNES DU TABLEAU REL                                
C         P      MATRICE DE PASSAGE                                             
C         RP     MATRICE ELEMENTAIRE 12X12 EN AXES LOCAUX                       
C                                                                               
C  SORTIE                                                                       
C         REL    MATRICE ELEMENTAIRE EN AXES GLOBAUX                            
C                                                                               
C 
C     Didier COMBESCURE Mars 2003                                                                
C-----------------------------------------------------------------------        
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
      DIMENSION REL(12,12),P(3,3),RP(12,12)                                     
      DIMENSION NI(12),NJ(12)                                                   
C                                                                               
      DATA NI/1,1,1,4,4,4,7,7,7,10,10,10/                                       
      DATA NJ/1,2,3,1,2,3,1,2,3,1,2,3/ 
C
       DO 10 I=1,12  
        DO 10 J=1,12
         REL(I,J)=  0.D0                
C                                        
         DO 10 K=1,3  
          DO 10 L=1,3
C
           I1 = NI(I) + K - 1
           J1 = NI(J) + L - 1
           I2 = NJ(I)
           J2 = NJ(J)
C
           REL(I,J)=REL(I,J) + (P(K,I2)*RP(I1,J1)*P(L,J2))
C   
   10     CONTINUE    
C
C                                                               
      RETURN                                                                    
      END