Esope source files

   1 : C ELPDI3    SOURCE    CHAT      05/01/12    23:36:39     5004           
   2 :       SUBROUTINE ELPDI3(P0,XN0                                                  
   3 :      &                  ,A,XN1,XN2,XT1,XT2,XLL                                  
   4 :      &                  ,XD,XNU                                                 
   5 :      &                  ,CS2,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 EN DYNAMIQUE PAR EQUATION INTEGRALE:                               
  12 : *                                                                               
  13 : *       CALCUL DU TERME DE COIN : 12 COEFFICIENTS A PLACER DANS LES             
  14 : *       RELATIONS                                                               
  15 : *                                                                               
  16 : ************************************************************************        
  17 : *                                                                               
  18 : *                                                                               
  19 :           DIMENSION P0(2)                                                       
  20 :           DIMENSION A (2)                                                       
  21 : *                                                                               
  22 :           DIMENSION XR (2)                                                      
  23 :           DIMENSION XN0(2)                                                      
  24 :           DIMENSION XN1(2)                                                      
  25 :           DIMENSION XN2(2)                                                      
  26 : *                                                                               
  27 :           DIMENSION XT1(2)                                                      
  28 :           DIMENSION XT2(2)                                                      
  29 : *                                                                               
  30 :           DIMENSION XLL(4)                                                      
  31 : *                                                                               
  32 :           DIMENSION XG (5 )                                                     
  33 :           DIMENSION CG (5 )                                                     
  34 :           DIMENSION CS2 (12)                                                    
  35 : *                                                                               
  36 : * 1-   CALCUL DE GREEN AU COIN                                                  
  37 : *                                                                               
  38 :         R = ( (A(1)-  P0(1))**2 + ( A(2) - P0(2)) ** 2 ) ** .5                  
  39 :         XR(1)=(A(1)-  P0(1) )/ R                                                
  40 :         XR(2)=(A(2)-  P0(2) )/ R                                                
  41 :         IF ( ISTAT .EQ. 1)THEN                                                  
  41 :           CALL ELPGFS(R,XG)                                                     
  42 :           DO 10 I=1,5                                                           
  43 :             CG(I) = XG(I)                                                       
  44 :  10       CONTINUE                                                              
  45 :         ELSE                                                                    
  45 :           CALL ELPDGF(R,CB,CG)                                                  
  46 :         ENDIF                                                                   
  47 : *                                                                               
  48 :                                                                                 
  49 : * 2-   COEFFICIENTS SUR LA PREMIERE LIGNE D'EQUATION                            
  50 : *                                                                               
  51 :          CAUX1 = -1* XD* (1. - XNU ) * CG(1)                                    
  52 :          CS2(1) =   CAUX1 / ( ( XLL(3) + XLL(4) ) / 2)                          
  53 :          CS2(2) = - CS2(1)                                                      
  54 :          CS2(3) = - CAUX1 / ( ( XLL(1) + XLL(2) ) / 2)                          
  55 :          CS2(4) = - CS2(3)                                                      
  56 : *                                                                               
  57 : *       CALCUL DE MNT EN A                                                      
  58 : *                                                                               
  58 :          CALL ELPDDR(CG,R,XR,XN1,XT1,XT1,XT1,2,CVA1)                            
  58 :          CALL ELPDDR(CG,R,XR,XN2,XT2,XT1,XT1,2,CVA2)                            
  59 :          CMNT1 = - 1* ( XD* (1. - XNU ) ) * CVA1                                
  60 :          CMNT2 = - 1* ( XD* (1. - XNU ) ) * CVA2                                
  61 :            AX =   XLL(2) / (XLL(2) + XLL(3))                                    
  62 :          CS2(5) =  AX * ( - CMNT2 + CMNT1)                                      
  63 :          CS2(6) = (1 -  AX) * ( - CMNT2 + CMNT1)                                
  64 :                                                                                 
  65 : *                                                                               
  66 : * 3-   COEFFICIENTS SUR LA DEUXIEME LIGNE                                       
  67 : *                                                                               
  67 :          CALL ELPDDR(CG,R,XR,XN0,XN0,XN0,XN0,1,CVA )                            
  68 :          CAUX1 = -1* XD* (1. - XNU ) * (- CVA)                                  
  69 :          CS2(7) =   CAUX1 / ( ( XLL(3) + XLL(4) ) / 2)                          
  70 :          CS2(8) = - CS2(7)                                                      
  71 :          CS2(9) = - CAUX1 / ( ( XLL(1) + XLL(2) ) / 2)                          
  72 :          CS2(10) = - CS2(9)                                                     
  73 : *                                                                               
  74 : *       CALCUL DE MNT,N0 EN A                                                   
  75 : *                                                                               
  75 :          CALL ELPDDR(CG,R,XR,XN1,XT1,XN0,XN0,3,CVA1)                            
  75 :          CALL ELPDDR(CG,R,XR,XN2,XT2,XN0,XN0,3,CVA2)                            
  76 :          CMNT1  = - 1* ( XD* (1. - XNU ) ) * (- CVA1)                           
  77 :          CMNT2  = - 1* ( XD* (1. - XNU ) ) * (- CVA2  )                         
  78 :          CS2(11) =  AX        * ( - CMNT2 + CMNT1)                              
  79 :          CS2(12) = (1. -  AX) * ( - CMNT2 + CMNT1)                              
  80 : *                                                                               
  81 :        RETURN                                                                   
  82 :       END                                                                       
  83 :