Sources Esope | Cast3M

elpdj2
C ELPDJ2    SOURCE    CHAT      05/01/12    23:36:47     5004
      SUBROUTINE ELPDJ2(XLP0
     &                  ,XD,XNU,CB,ISTAT,NTRAP2
     &                  ,CS1)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-B,D-H,O-Z)
      IMPLICIT COMPLEX*16(C)
************************************************************************
*
*    PLAQUES EN DYN      PAR EQUATION INTEGRALE:
*
*       CALCUL DES TERMES "DIAGONAUX"
*        V  ..        V,N0   0
*        VN 0        VN,N0   ..
*        MN ..       MN,N0   0
*        KN 0        KN,N0   ****** voir elpdm1
*
************************************************************************
          DIMENSION CS0 (8)
          DIMENSION CS1 (8)
          DIMENSION  S1 (8)
*
*  1-- CALCUL ANALYTIQUE DES TERMES STATIQUES  (sauf kn,n0)
*
        CALL ELPI02(XLP0,XD,XNU,S1)
*
*  2-- CALCUL PAR TRAPEZE SUR op(w) -op(ws)
         IF ( ISTAT.NE.1) THEN
           CALL ELPDJ1(XLP0,XD,XNU,NTRAP2,CS0,CB)
         ENDIF
*  3-- ON AJOUTE
*
         CS1(1) = CS0(1) + S1(1)
         CS1(2) = 0D0
         CS1(3) = CS0(3) + S1(3)
         CS1(4) = 0D0
         CS1(5) = 0D0
         CS1(6) = CS0(6) + S1(6)
         CS1(7) = 0D0
         CS1(8) = CS0(8) + S1(8)
*
       RETURN
      END
*