Esope source files

   1 : C POUJAC    SOURCE    CHAT      05/01/13    02:18:57     5004           
   2 :       SUBROUTINE POUJAC(COORD,XL)                                               
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE CALCUL DU JACOBIEN DE LA POUTRE DROITE                            
   5 : C                                                                               
   6 : C  ENTREE                                                                       
   7 : C          COORD   COORDONNEES DES 2 NOEUDS                                     
   8 : C                                                                               
   9 : C  SORTIE                                                                       
  10 : C          XL      JACOBIEN ( = LONGUEUR - CTE DANS L'ELEMENT )                 
  11 : C-----------------------------------------------------------------------        
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  14 :       DIMENSION COORD(3,1),XE(2),YE(2),ZE(2)                                    
  15 : C                                                                               
  16 :       DO 2 I=1,2                                                                
  17 :       XE(I)=COORD(1,I)                                                          
  18 :       YE(I)=COORD(2,I)                                                          
  19 :   2   ZE(I)=COORD(3,I)                                                          
  20 :       XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2                     
  21 :       XL=SQRT(XL)                                                               
  22 :       RETURN                                                                    
  23 :       END                                                                       
  24 :