Télécharger la plaquette Cast3M


fr en

Télécharger elpdi1.eso

Retour à la liste

Numérotation des lignes : 

elpdi1

  1. C ELPDI1    SOURCE    CHAT      05/01/12    23:36:35     5004
  2.       SUBROUTINE ELPDI1(P0,XN0,XT0
  3.      &                  ,Q,XN,XT,A,B,OAB,TETA,R,XLQ
  4.      &                  ,XD,XNU
  5.      &                  ,NTRAP,CS1,CB,ISTAT)
  6.       IMPLICIT INTEGER(I-N)
  7.       IMPLICIT REAL*8(A-B,D-H,O-Z)
  8.       IMPLICIT COMPLEX*16(C)
  9. ************************************************************************
  10. *
  11. *    PLAQUES PAR EQUATION INTEGRALE:
  12. *
  13. *       ON INTEGRE PAR LA METHODE DES TRAPEZES POUR CALCULER
  14. *       LES COEFFICIENTS NON  DIAGONAUX
  15. *
  16. ************************************************************************
  17. *
  18. *
  19.           DIMENSION P0(2)
  20.           DIMENSION A (2)
  21.           DIMENSION B (2)
  22.           DIMENSION OAB(2)
  23.           DIMENSION Q (2)
  24.           DIMENSION QI(2)
  25.           DIMENSION COP(8)
  26. *
  27.           DIMENSION XN (2)
  28.           DIMENSION XNQI(2)
  29.           DIMENSION XN0(2)
  30. *
  31.           DIMENSION XT (2)
  32.           DIMENSION XTQI(2)
  33.           DIMENSION XT0(2)
  34. *
  35. *
  36.           DIMENSION CS1 (8)
  37. *
  38. * 1-   TAILLE DES TRAPEZES  ARC DE CERCLE
  39. *
  40.         X0 = XLQ/(2*NTRAP)
  41.         X1 = (2*TETA*R) / NTRAP
  42. *
  43. * 2-   INTEGRATION
  44. *
  45.          DO 10 J=1,8
  46.              CS1(J) = 0D0
  47. 10       CONTINUE
  48.          DO 100 I=1,NTRAP
  49.            TETAI = ( (I-1)*2 + 1) * TETA /NTRAP
  50.  
  51.             QI(1) = OAB(1) + COS(TETAI) *( A(1) -OAB(1))
  52.      &                     - SIN(TETAI) *( A(2) -OAB(2))
  53.             QI(2) = OAB(2) + SIN(TETAI) *( A(1) -OAB(1))
  54.      &                     + COS(TETAI) *( A(2) -OAB(2))
  55. *         -- NORMALE
  56.             XNQI(1) = ( QI(1) - OAB(1) ) / R
  57.             XNQI(2) = ( QI(2) - OAB(2) ) / R
  58.             PS = XNQI(1) * XN(1) + XNQI(2) * XN(2)
  59.             XNQI(1) = ( PS / ABS(PS)) * XNQI(1)
  60.             XNQI(2) = ( PS / ABS(PS)) * XNQI(2)
  61.  
  62. *         -- TANGENTE
  63.             XTQI(1) = -1*XNQI(2)
  64.             XTQI(2) =    XNQI(1)
  65.  
  66.            CALL ELPDOP(P0,XN0,QI,XNQI,XTQI,XD,XNU,COP,CB,ISTAT)
  67.            DO 110 J=1,8
  68.              CS1(J) = CS1(J) + ( X1 * COP(J))
  69. 110        CONTINUE
  70. 100      CONTINUE
  71. *
  72.        RETURN
  73.       END
  74.  
  75.  

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