Numérotation des lignes :

elpdi1
C ELPDI1    SOURCE    CHAT      05/01/12    23:36:35     5004      SUBROUTINE ELPDI1(P0,XN0,XT0     &                  ,Q,XN,XT,A,B,OAB,TETA,R,XLQ     &                  ,XD,XNU     &                  ,NTRAP,CS1,CB,ISTAT)      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-B,D-H,O-Z)      IMPLICIT COMPLEX*16(C)**************************************************************************    PLAQUES PAR EQUATION INTEGRALE:**       ON INTEGRE PAR LA METHODE DES TRAPEZES POUR CALCULER*       LES COEFFICIENTS NON  DIAGONAUX***************************************************************************          DIMENSION P0(2)          DIMENSION A (2)          DIMENSION B (2)          DIMENSION OAB(2)          DIMENSION Q (2)          DIMENSION QI(2)          DIMENSION COP(8)*          DIMENSION XN (2)          DIMENSION XNQI(2)          DIMENSION XN0(2)*          DIMENSION XT (2)          DIMENSION XTQI(2)          DIMENSION XT0(2)**          DIMENSION CS1 (8)** 1-   TAILLE DES TRAPEZES  ARC DE CERCLE*        X0 = XLQ/(2*NTRAP)        X1 = (2*TETA*R) / NTRAP** 2-   INTEGRATION*         DO 10 J=1,8             CS1(J) = 0D010       CONTINUE         DO 100 I=1,NTRAP           TETAI = ( (I-1)*2 + 1) * TETA /NTRAP             QI(1) = OAB(1) + COS(TETAI) *( A(1) -OAB(1))     &                     - SIN(TETAI) *( A(2) -OAB(2))            QI(2) = OAB(2) + SIN(TETAI) *( A(1) -OAB(1))     &                     + COS(TETAI) *( A(2) -OAB(2))*         -- NORMALE            XNQI(1) = ( QI(1) - OAB(1) ) / R            XNQI(2) = ( QI(2) - OAB(2) ) / R            PS = XNQI(1) * XN(1) + XNQI(2) * XN(2)            XNQI(1) = ( PS / ABS(PS)) * XNQI(1)            XNQI(2) = ( PS / ABS(PS)) * XNQI(2) *         -- TANGENTE            XTQI(1) = -1*XNQI(2)            XTQI(2) =    XNQI(1)            CALL ELPDOP(P0,XN0,QI,XNQI,XTQI,XD,XNU,COP,CB,ISTAT)           DO 110 J=1,8             CS1(J) = CS1(J) + ( X1 * COP(J))110        CONTINUE100      CONTINUE*       RETURN      END  

© Cast3M 2003 - Tous droits réservés.
Mentions légales