Sources Esope | Cast3M

elpdj1
C ELPDJ1    SOURCE    CHAT      05/01/12    23:36:43     5004
      SUBROUTINE ELPDJ1(XLP0
     &                  ,XD,XNU
     &                  ,NTRAP2,CS1,CB)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-B,D-H,O-Z)
      IMPLICIT COMPLEX*16(C)
************************************************************************
*
*    PLAQUES PAR EQUATION INTEGRALE:
*
*       ON INTEGRE PAR LA METHODE DES TRAPEZES EN STATIQUE
*       ET EN DYNAMIQUE
*
************************************************************************
*
*
          DIMENSION P0(2)
          DIMENSION QI(2)
          DIMENSION XOP(8)
          DIMENSION COPD(8)
*
          DIMENSION XN (2)
*
          DIMENSION XT (2)
*
*
          DIMENSION CS1 (8)
*
* 1-   TAILLE DES TRAPEZES ( DEMI SEGMENT SEULEMENT)
*
        X0 = XLP0/(4*NTRAP2)
        X1 = XLP0 /(2*NTRAP2)
*
* 2-   POSITIONS DE REFERENCE
*
        P0(1) = 0D0
        P0(2) = 0D0
        XN(1) = 1D0
        XN(2) = 0D0
        XT(1) = 0D0
        XT(2) = 1D0
        QI(1) = 0D0
        QI(2) = 0D0
*
* 2-   INTEGRATION
*
         DO 10 J=1,8
             CS1(J) = (0d0,0d0)
10       CONTINUE
         DO 100 I=1,NTRAP2
           QI(2) = X0 + (I-1)*X1
           istat = 0
           CALL ELPDOP(P0,XN,QI,XN,XT,XD,XNU,COPD,CB,istat)
 
           CALL ELPOP1(P0,XN,QI,XN,XT,XD,XNU,XOP)
 
 
 
 
           DO 110 J=1,8
             ccop =xop(j)*(1d0,0d0)
             CS1(J) = CS1(J) + 2*( X1 * (COPD(J)- ccop))
110        CONTINUE
100      CONTINUE
*
       RETURN
      END