Esope source files

   1 : C ELPDJ1    SOURCE    CHAT      05/01/12    23:36:43     5004           
   2 :       SUBROUTINE ELPDJ1(XLP0                                                    
   3 :      &                  ,XD,XNU                                                 
   4 :      &                  ,NTRAP2,CS1,CB)                                         
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-B,D-H,O-Z)                                              
   7 :       IMPLICIT COMPLEX*16(C)                                                    
   8 : ************************************************************************        
   9 : *                                                                               
  10 : *    PLAQUES PAR EQUATION INTEGRALE:                                            
  11 : *                                                                               
  12 : *       ON INTEGRE PAR LA METHODE DES TRAPEZES EN STATIQUE                      
  13 : *       ET EN DYNAMIQUE                                                         
  14 : *                                                                               
  15 : ************************************************************************        
  16 : *                                                                               
  17 : *                                                                               
  18 :           DIMENSION P0(2)                                                       
  19 :           DIMENSION QI(2)                                                       
  20 :           DIMENSION XOP(8)                                                      
  21 :           DIMENSION COPD(8)                                                     
  22 : *                                                                               
  23 :           DIMENSION XN (2)                                                      
  24 : *                                                                               
  25 :           DIMENSION XT (2)                                                      
  26 : *                                                                               
  27 : *                                                                               
  28 :           DIMENSION CS1 (8)                                                     
  29 : *                                                                               
  30 : * 1-   TAILLE DES TRAPEZES ( DEMI SEGMENT SEULEMENT)                            
  31 : *                                                                               
  32 :         X0 = XLP0/(4*NTRAP2)                                                    
  33 :         X1 = XLP0 /(2*NTRAP2)                                                   
  34 : *                                                                               
  35 : * 2-   POSITIONS DE REFERENCE                                                   
  36 : *                                                                               
  37 :         P0(1) = 0D0                                                             
  38 :         P0(2) = 0D0                                                             
  39 :         XN(1) = 1D0                                                             
  40 :         XN(2) = 0D0                                                             
  41 :         XT(1) = 0D0                                                             
  42 :         XT(2) = 1D0                                                             
  43 :         QI(1) = 0D0                                                             
  44 :         QI(2) = 0D0                                                             
  45 : *                                                                               
  46 : * 2-   INTEGRATION                                                              
  47 : *                                                                               
  48 :          DO 10 J=1,8                                                            
  49 :              CS1(J) = (0d0,0d0)
  50 : 10       CONTINUE                                                               
  51 :          DO 100 I=1,NTRAP2                                                      
  52 :            QI(2) = X0 + (I-1)*X1                                                
  53 :            istat = 0
  53 :            CALL ELPDOP(P0,XN,QI,XN,XT,XD,XNU,COPD,CB,istat)
  54 :                                                                                 
  54 :            CALL ELPOP1(P0,XN,QI,XN,XT,XD,XNU,XOP)                               
  55 :                                                                                 
  56 :                                                                                 
  57 :                                                                                 
  58 :                                                                                 
  59 :            DO 110 J=1,8                                                         
  60 :              ccop =xop(j)*(1d0,0d0)
  61 :              CS1(J) = CS1(J) + 2*( X1 * (COPD(J)- ccop))
  62 : 110        CONTINUE                                                             
  63 : 100      CONTINUE                                                               
  64 : *                                                                               
  65 :        RETURN                                                                   
  66 :       END                                                                       
  67 :