Numérotation des lignes :

C ELPDKP    SOURCE    CHAT      05/01/12    23:36:58     5004      SUBROUTINE ELPDKP(P0,XN0,QI,XN,XT,XD,XNU,COP,CB)      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-B,D-H,O-Z)      IMPLICIT COMPLEX*16(C)**************************************************************************    PLAQUES PAR EQUATION INTEGRALE:**       CALCUL DES OPERATEURS PLAQUE ET DE LEURS DERIVEES*       SUR WD - WS***************************************************************************          DIMENSION P0(2)          DIMENSION QI(2)          DIMENSION COP(8)*          DIMENSION XN (2)          DIMENSION XN0(2)*          DIMENSION XT (2)          DIMENSION XT0(2)          DIMENSION CG (5)          DIMENSION CG0(5)          DIMENSION XG (5)          DIMENSION XR (2)** 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R*         R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5         XR(1) = (QI(1) - P0(1)) / R         XR(2) = (QI(2) - P0(2)) / R***** 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R*       CALL ELPDGF(R,CB,CG0)       CALL ELPGFS(R,XG)        DO 10 I=1,5          CG(I) = CG0(I) - XG(I)10      CONTINUE** 3-   CALCUL DE V*       COP(1) = CG(1)** 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)*       CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,2,CAUX1)       CALL ELPDDR(CG,R,XR,XT,XT,XN,XN,2,CAUX2)       COP(3) = -1* XD* ( CAUX1 +  XNU*CAUX2)** 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)* ON NE LE CALCULE PAS*       CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,3,CAUX1)*       CALL ELPDDR(CG,R,XR,XN,XT,XT,XN,3,CAUX2)*       COP(4) = -1* XD* ( CAUX1 +  (2. - XNU) *CAUX2)******      TERMES DERIVES******* 8-   CALCUL DE V,NN0*       CALL ELPDDR(CG,R,XR,XN,XN0,XN,XN,2,CAUX1)       COP(6) =  - CAUX1** 10-  CALCUL DE KN,N0*       CALL ELPDDR(CG,R,XR,XN,XN,XN,XN0,4,CAUX1)       CALL ELPDDR(CG,R,XR,XN,XT,XT,XN0,4,CAUX2)       COP(8) = -1* XD* (-  CAUX1 -  (2. - XNU) *CAUX2)        RETURN      END

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