C COUROT    SOURCE    CHAT      05/01/12    22:25:39     5004           
      SUBROUTINE COUROT(REL,LRE,P,RP)                                           
C-----------------------------------------------------------------------        
C  ROUTINE DE PASSAGE D UNE MATRICE ELEMENTAIRE DES AXES LOCAUX AUX AXES        
C  GLOBAUX                                                                      
C
C  REL = TRANSPOSE(P).RP.P
C                                                                               
C  ENTREE                                                                       
C         LRE    NOMBRE DE LIGNES DU TABLEAU REL                                
C         P      MATRICE DE PASSAGE                                             
C         RP     MATRICE ELEMENTAIRE 10X10 EN AXES LOCAUX                       
C                                                                               
C  SORTIE                                                                       
C         REL    MATRICE ELEMENTAIRE EN AXES GLOBAUX                            
C-----------------------------------------------------------------------        
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)                                                  
      DIMENSION REL(LRE,LRE),P(LRE,LRE),RP(LRE,LRE),RPP(10,10)
C
      DO 1 I=1,100
1     RPP(I,1)=0.
C
      DO 2 I=1,10
         DO 2 J=1,10
            DO 3 K=1,10
3           RPP(I,J)=RPP(I,J)+RP(I,K)*P(K,J)                        
2     CONTINUE
C
      DO 4 I=1,10
         DO 4 J=1,10
            DO 5 K=1,10
5           REL(I,J)=REL(I,J)+P(K,I)*RPP(K,J)                        
4     CONTINUE
C
      RETURN
      END