1 : C ELPDGF SOURCE CHAT 05/01/12 23:36:29 5004 2 : SUBROUTINE ELPDGF(R,CB,CG) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-B,D-H,O-Z) 5 : IMPLICIT COMPLEX*16(C) 6 : C 7 : C ===================================================================== 8 : C ELFE 9 : C CALCUL DE LA FONCTION DE GREEN ET DE DES DERIVEES 10 : C POUR UNE PLAQUE EN DYNAMIQUE 11 : C 12 : C ===================================================================== 13 : C 14 : C 15 : C 16 : DIMENSION CG(5) 17 : DIMENSION CZ(1) 18 : DIMENSION CGF0(1) 19 : DIMENSION CGF1(1) 20 : DIMENSION CGF2(1) 21 : DIMENSION CGF3(1) 22 : DIMENSION CGF4(1) 23 : C 24 : C ON UTILISE LES CGFI COMME TABLEAUX DE TRAVAIL 25 : C POUR CALCULER LES H10,H11,H10(IBR),H11(IBR) 26 : C 27 : np = 10 28 : xm = 5d0 29 : CI = (0d0,1d0) 30 : CZ(1) = CB* R 31 : N=1 31 : CALL ELPDBE(CZ ,N,np,CGF0,CGF1,xm) 32 : DO 100 I = 1,N 33 : CZ(I) = CI*CZ(I) 34 : 100 CONTINUE 34 : CALL ELPDBE(CZ ,N,np,CGF2,CGF3,xm) 35 : DO 110 I = 1,N 36 : CZ(I) = CZ(I)/CI 37 : 110 CONTINUE 38 : 39 : DO 150 I = 1,N 40 : CH0M= CGF2(I) - CGF0(I) 41 : CH0P= CGF2(I) + CGF0(I) 42 : CH1M= CI * CGF3(I) - CGF1(I) 43 : CH1P= CI * CGF3(I) + CGF1(I) 44 : CGF0(I) = CH0M 45 : CGF1(I) =-1.* CH1M 46 : CGF2(I) = CH0P + ( CH1M / CZ(I) ) 47 : CGF3(I) = -1/CZ(I) *CH0P 48 : & -1 *CH1P 49 : & -2 /(CZ(I)** 2) *CH1M 50 : CGF4(I) = CH0M 51 : & +3 /(CZ(I)** 2) *CH0P 52 : & +2 / CZ(I) *CH1P 53 : & +6 /(CZ(I)** 3) *CH1M 54 : 150 CONTINUE 55 : C 56 : CG(1)= -1*( CGF0(1))*CI/(8*CB*CB) 57 : CG(2)= -1*CGF1(1)*CI/(8*CB) 58 : CG(3)= -1*CGF2(1)*CI/8 59 : CG(4)= -1*CGF3(1)*CI*(CB)/8 60 : CG(5)= -1*CGF4(1)*CI*(CB*CB)/8 61 : C 62 : RETURN 63 : END 64 :
© Cast3M 2003 - All rights reserved.
Disclaimer