Esope source files

   1 : C ELPDI1    SOURCE    CHAT      05/01/12    23:36:35     5004           
   2 :       SUBROUTINE ELPDI1(P0,XN0,XT0                                              
   3 :      &                  ,Q,XN,XT,A,B,OAB,TETA,R,XLQ                             
   4 :      &                  ,XD,XNU                                                 
   5 :      &                  ,NTRAP,CS1,CB,ISTAT)                                    
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-B,D-H,O-Z)                                              
   8 :       IMPLICIT COMPLEX*16(C)                                                    
   9 : ************************************************************************        
  10 : *                                                                               
  11 : *    PLAQUES PAR EQUATION INTEGRALE:                                            
  12 : *                                                                               
  13 : *       ON INTEGRE PAR LA METHODE DES TRAPEZES POUR CALCULER                    
  14 : *       LES COEFFICIENTS NON  DIAGONAUX                                         
  15 : *                                                                               
  16 : ************************************************************************        
  17 : *                                                                               
  18 : *                                                                               
  19 :           DIMENSION P0(2)                                                       
  20 :           DIMENSION A (2)                                                       
  21 :           DIMENSION B (2)                                                       
  22 :           DIMENSION OAB(2)                                                      
  23 :           DIMENSION Q (2)                                                       
  24 :           DIMENSION QI(2)                                                       
  25 :           DIMENSION COP(8)                                                      
  26 : *                                                                               
  27 :           DIMENSION XN (2)                                                      
  28 :           DIMENSION XNQI(2)                                                     
  29 :           DIMENSION XN0(2)                                                      
  30 : *                                                                               
  31 :           DIMENSION XT (2)                                                      
  32 :           DIMENSION XTQI(2)                                                     
  33 :           DIMENSION XT0(2)                                                      
  34 : *                                                                               
  35 : *                                                                               
  36 :           DIMENSION CS1 (8)                                                     
  37 : *                                                                               
  38 : * 1-   TAILLE DES TRAPEZES  ARC DE CERCLE                                       
  39 : *                                                                               
  40 :         X0 = XLQ/(2*NTRAP)                                                      
  41 :         X1 = (2*TETA*R) / NTRAP                                                 
  42 : *                                                                               
  43 : * 2-   INTEGRATION                                                              
  44 : *                                                                               
  45 :          DO 10 J=1,8                                                            
  46 :              CS1(J) = 0D0                                                       
  47 : 10       CONTINUE                                                               
  48 :          DO 100 I=1,NTRAP                                                       
  49 :            TETAI = ( (I-1)*2 + 1) * TETA /NTRAP                                 
  50 :                                                                                 
  51 :             QI(1) = OAB(1) + COS(TETAI) *( A(1) -OAB(1))                        
  52 :      &                     - SIN(TETAI) *( A(2) -OAB(2))                        
  53 :             QI(2) = OAB(2) + SIN(TETAI) *( A(1) -OAB(1))                        
  54 :      &                     + COS(TETAI) *( A(2) -OAB(2))                        
  55 : *         -- NORMALE                                                            
  56 :             XNQI(1) = ( QI(1) - OAB(1) ) / R                                    
  57 :             XNQI(2) = ( QI(2) - OAB(2) ) / R                                    
  58 :             PS = XNQI(1) * XN(1) + XNQI(2) * XN(2)                              
  59 :             XNQI(1) = ( PS / ABS(PS)) * XNQI(1)                                 
  60 :             XNQI(2) = ( PS / ABS(PS)) * XNQI(2)                                 
  61 :                                                                                 
  62 : *         -- TANGENTE                                                           
  63 :             XTQI(1) = -1*XNQI(2)                                                
  64 :             XTQI(2) =    XNQI(1)                                                
  65 :                                                                                 
  65 :            CALL ELPDOP(P0,XN0,QI,XNQI,XTQI,XD,XNU,COP,CB,ISTAT)                 
  66 :            DO 110 J=1,8                                                         
  67 :              CS1(J) = CS1(J) + ( X1 * COP(J))                                   
  68 : 110        CONTINUE                                                             
  69 : 100      CONTINUE                                                               
  70 : *                                                                               
  71 :        RETURN                                                                   
  72 :       END                                                                       
  73 :