Sources Esope | Cast3M

elpdk2
C ELPDK2    SOURCE    CHAT      05/01/12    23:36:54     5004
      SUBROUTINE ELPDK2(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)
*
          DO 10 I=1,8
           S1(I) = 0D0
           CS0(I) = 0D0
 10      CONTINUE
*
*  1-- CALCUL ANALYTIQUE DES TERMES STATIQUES  sauf knn0
*
        CALL ELPI02(XLP0,XD,XNU,S1)
*
*  2-- CALCUL PAR TRAPEZE SUR op(W - WS)
*
       IF ( ISTAT.NE.1) THEN
          CALL ELPDK1(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
*