Esope source files

   1 : C ELPDK2    SOURCE    CHAT      05/01/12    23:36:54     5004           
   2 :       SUBROUTINE ELPDK2(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 :           DO 10 I=1,8                                                           
  24 :            S1(I) = 0D0                                                          
  25 :            CS0(I) = 0D0                                                         
  26 :  10      CONTINUE                                                               
  27 : *                                                                               
  28 : *  1-- CALCUL ANALYTIQUE DES TERMES STATIQUES  sauf knn0
  29 : *                                                                               
  29 :         CALL ELPI02(XLP0,XD,XNU,S1)                                             
  30 : *                                                                               
  31 : *  2-- CALCUL PAR TRAPEZE SUR op(W - WS)
  32 : *                                                                               
  33 :        IF ( ISTAT.NE.1) THEN                                                    
  33 :           CALL ELPDK1(XLP0,XD,XNU,NTRAP2,CS0,CB)                                
  34 :        ENDIF                                                                    
  35 : *                                                                               
  36 : *  3-- ON AJOUTE                                                                
  37 : *                                                                               
  38 :          CS1(1) = CS0(1) + S1(1)                                                
  39 :          CS1(2) = 0D0                                                           
  40 :          CS1(3) = CS0(3) + S1(3)                                                
  41 :          CS1(4) = 0D0                                                           
  42 :          CS1(5) = 0D0                                                           
  43 :          CS1(6) = CS0(6) + S1(6)                                                
  44 :          CS1(7) = 0D0                                                           
  45 :          CS1(8) = CS0(8) + S1(8)                                                
  46 : *                                                                               
  47 :        RETURN                                                                   
  48 :       END                                                                       
  49 : *                                                                               
  50 :