Sources Esope | Cast3M

poujac
C POUJAC    SOURCE    PV        18/06/18    21:15:27     9860           
      SUBROUTINE POUJAC(COORD,XL)
C-----------------------------------------------------------------------
C  ROUTINE DE CALCUL DU JACOBIEN DE LA POUTRE DROITE
C
C  ENTREE
C          COORD   COORDONNEES DES 2 NOEUDS
C
C  SORTIE
C          XL      JACOBIEN ( = LONGUEUR - CTE DANS L'ELEMENT )
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION COORD(3,*),XE(2),YE(2),ZE(2)
C
      DO 2 I=1,2
      XE(I)=COORD(1,I)
      YE(I)=COORD(2,I)
  2   ZE(I)=COORD(3,I)
      XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2
      XL=SQRT(XL)
      RETURN
      END