Télécharger la plaquette Cast3M


fr en

Télécharger elpdop.eso

Retour à la liste

Numérotation des lignes : 

elpdop

  1. C ELPDOP    SOURCE    KK2000    14/04/09    21:15:25     8027
  2.       SUBROUTINE ELPDOP(P0,XN0,QI,XN,XT,XD,XNU,COP,CB,ISTAT)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-B,D-H,O-Z)
  5.       IMPLICIT COMPLEX*16(C)
  6. ************************************************************************
  7. *
  8. *    PLAQUES PAR EQUATION INTEGRALE:
  9. *
  10. *       CALCUL DES OPERATEURS PLAQUE ET DE LEURS DERIVEES
  11. *       POUR UN COUPLE DE PTS P0 , QI
  12. *
  13. ************************************************************************
  14. *
  15. *
  16.           DIMENSION P0(2)
  17.           DIMENSION QI(2)
  18.           DIMENSION COP(8)
  19. *
  20.           DIMENSION XN (2)
  21.           DIMENSION XN0(2)
  22. *
  23.           DIMENSION XT (2)
  24. C          DIMENSION XT0(2)
  25.           DIMENSION CG (5)
  26.           DIMENSION XG (5)
  27.           DIMENSION XR (2)
  28. *
  29. * 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R
  30. *
  31.          R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5
  32.          XR(1) = (QI(1) - P0(1)) / R
  33.          XR(2) = (QI(2) - P0(2)) / R
  34. *
  35. *
  36. *
  37. *
  38. * 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R
  39. *
  40.        IF ( ISTAT.EQ.1) THEN
  41.          CALL ELPGFS(R,XG)
  42.          DO  10 I=1,5
  43.            CG(I) = XG(I)*(1d0,0d0)
  44.  10      CONTINUE
  45.        ELSE
  46.          CALL ELPDGF(R,CB,CG)
  47.        ENDIF
  48. *
  49. * 3-   CALCUL DE V
  50. *
  51.        COP(1) = CG(1)
  52. *
  53. * 4-   CALCUL DE V,N
  54. *
  55.        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,1,CAUX1)
  56.        COP(2) = CAUX1
  57. *
  58. * 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)
  59. *
  60.        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,2,CAUX1)
  61.        CALL ELPDDR(CG,R,XR,XT,XT,XN,XN,2,CAUX2)
  62.        COP(3) = CMPLX(-1d0)* XD* ( CAUX1 +  XNU*CAUX2)
  63. *
  64. * 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)
  65. *
  66.        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,3,CAUX1)
  67.        CALL ELPDDR(CG,R,XR,XN,XT,XT,XN,3,CAUX2)
  68.        COP(4) = CMPLX(-1d0)* XD* ( CAUX1 +  (CMPLX(2d0) - XNU) *CAUX2)
  69. *****
  70. *      TERMES DERIVES
  71. *****
  72. * 7-   CALCUL DE V,N0
  73. *
  74. *
  75. *
  76. *
  77.        CALL ELPDDR(CG,R,XR,XN0,XN,XN,XN,1,CAUX1)
  78.        COP(5) = - CAUX1
  79. *
  80. * 8-   CALCUL DE V,NN0
  81. *
  82.        CALL ELPDDR(CG,R,XR,XN,XN0,XN,XN,2,CAUX1)
  83.        COP(6) =  - CAUX1
  84. *
  85. * 9-   CALCUL DE MN,N0
  86. *
  87.        CALL ELPDDR(CG,R,XR,XN,XN,XN0,XN,3,CAUX1)
  88.        CALL ELPDDR(CG,R,XR,XT,XT,XN0,XN,3,CAUX2)
  89.        COP(7) = CMPLX(-1d0)* XD* ( - CAUX1 -  XNU*CAUX2)
  90. *
  91. * 9-   CALCUL DE KN,N0
  92. *
  93.        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN0,4,CAUX1)
  94.        CALL ELPDDR(CG,R,XR,XN,XT,XT,XN0,4,CAUX2)
  95.        COP(8) = CMPLX(-1d0)* XD* (-  CAUX1 -  (CMPLX(2d0) - XNU) *CAUX2)
  96.  
  97.        RETURN
  98.       END
  99.  
  100.  
  101.  

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