Esope source files

   1 : C ELPDJ2    SOURCE    CHAT      05/01/12    23:36:47     5004           
   2 :       SUBROUTINE ELPDJ2(XLP0                                                    
   3 :      &                  ,XD,XNU,CB,ISTAT,NTRAP2                                 
   4 :      &                  ,CS1)                                                   
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-B,D-H,O-Z)                                              
   7 :       IMPLICIT COMPLEX*16(C)                                                    
   8 : ************************************************************************        
   9 : *                                                                               
  10 : *    PLAQUES EN DYN      PAR EQUATION INTEGRALE:                                
  11 : *                                                                               
  12 : *       CALCUL DES TERMES "DIAGONAUX"                                           
  13 : *        V  ..        V,N0   0                                                  
  14 : *        VN 0        VN,N0   ..                                                 
  15 : *        MN ..       MN,N0   0                                                  
  16 : *        KN 0        KN,N0   ****** voir elpdm1
  17 : *                                                                               
  18 : ************************************************************************        
  19 :           DIMENSION CS0 (8)                                                     
  20 :           DIMENSION CS1 (8)                                                     
  21 :           DIMENSION  S1 (8)                                                     
  22 : *                                                                               
  23 : *  1-- CALCUL ANALYTIQUE DES TERMES STATIQUES  (sauf kn,n0)
  24 : *
  24 :         CALL ELPI02(XLP0,XD,XNU,S1)                                             
  25 : *                                                                               
  26 : *  2-- CALCUL PAR TRAPEZE SUR op(w) -op(ws)
  27 :          IF ( ISTAT.NE.1) THEN                                                  
  27 :            CALL ELPDJ1(XLP0,XD,XNU,NTRAP2,CS0,CB)                               
  28 :          ENDIF                                                                  
  29 : *  3-- ON AJOUTE                                                                
  30 : *                                                                               
  31 :          CS1(1) = CS0(1) + S1(1)                                                
  32 :          CS1(2) = 0D0                                                           
  33 :          CS1(3) = CS0(3) + S1(3)                                                
  34 :          CS1(4) = 0D0                                                           
  35 :          CS1(5) = 0D0                                                           
  36 :          CS1(6) = CS0(6) + S1(6)                                                
  37 :          CS1(7) = 0D0                                                           
  38 :          CS1(8) = CS0(8) + S1(8)                                                
  39 : *                                                                               
  40 :        RETURN                                                                   
  41 :       END                                                                       
  42 : *                                                                               
  43 :