Sources Esope | Cast3M

elpdgf
C ELPDGF    SOURCE    CHAT      05/01/12    23:36:29     5004
      SUBROUTINE ELPDGF(R,CB,CG)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-B,D-H,O-Z)
      IMPLICIT COMPLEX*16(C)
C
C =====================================================================
C     ELFE
C     CALCUL DE LA FONCTION DE GREEN ET DE DES DERIVEES
C     POUR UNE PLAQUE EN DYNAMIQUE
C
C =====================================================================
C
C
C
      DIMENSION CG(5)
      DIMENSION CZ(1)
      DIMENSION CGF0(1)
      DIMENSION CGF1(1)
      DIMENSION CGF2(1)
      DIMENSION CGF3(1)
      DIMENSION CGF4(1)
C
C ON UTILISE LES CGFI COMME TABLEAUX DE TRAVAIL
C POUR CALCULER LES H10,H11,H10(IBR),H11(IBR)
C
       np = 10
       xm    = 5d0
       CI = (0d0,1d0)
       CZ(1) = CB* R
       N=1
       CALL ELPDBE(CZ ,N,np,CGF0,CGF1,xm)
       DO 100 I = 1,N
         CZ(I) = CI*CZ(I)
100   CONTINUE
       CALL ELPDBE(CZ ,N,np,CGF2,CGF3,xm)
       DO 110 I = 1,N
         CZ(I) = CZ(I)/CI
110   CONTINUE
 
       DO 150 I = 1,N
        CH0M= CGF2(I) - CGF0(I)
        CH0P= CGF2(I) + CGF0(I)
        CH1M= CI * CGF3(I) - CGF1(I)
        CH1P= CI * CGF3(I) + CGF1(I)
        CGF0(I) = CH0M
        CGF1(I) =-1.* CH1M
        CGF2(I) = CH0P    + ( CH1M / CZ(I) )
        CGF3(I) = -1/CZ(I)       *CH0P
     &                     -1    *CH1P
     &           -2 /(CZ(I)** 2) *CH1M
        CGF4(I) =                 CH0M
     &           +3 /(CZ(I)** 2) *CH0P
     &           +2 / CZ(I)      *CH1P
     &           +6 /(CZ(I)** 3) *CH1M
150   CONTINUE
C
       CG(1)= -1*( CGF0(1))*CI/(8*CB*CB)
       CG(2)= -1*CGF1(1)*CI/(8*CB)
       CG(3)= -1*CGF2(1)*CI/8
       CG(4)= -1*CGF3(1)*CI*(CB)/8
       CG(5)= -1*CGF4(1)*CI*(CB*CB)/8
C
      RETURN
      END