Esope source files

   1 : C ELPDKP    SOURCE    CHAT      05/01/12    23:36:58     5004           
   2 :       SUBROUTINE ELPDKP(P0,XN0,QI,XN,XT,XD,XNU,COP,CB)                          
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-B,D-H,O-Z)                                              
   5 :       IMPLICIT COMPLEX*16(C)                                                    
   6 : ************************************************************************        
   7 : *                                                                               
   8 : *    PLAQUES PAR EQUATION INTEGRALE:                                            
   9 : *                                                                               
  10 : *       CALCUL DES OPERATEURS PLAQUE ET DE LEURS DERIVEES                       
  11 : *       SUR WD - WS                                                             
  12 : *                                                                               
  13 : ************************************************************************        
  14 : *                                                                               
  15 : *                                                                               
  16 :           DIMENSION P0(2)                                                       
  17 :           DIMENSION QI(2)                                                       
  18 :           DIMENSION COP(8)                                                      
  19 : *                                                                               
  20 :           DIMENSION XN (2)                                                      
  21 :           DIMENSION XN0(2)                                                      
  22 : *                                                                               
  23 :           DIMENSION XT (2)                                                      
  24 :           DIMENSION XT0(2)                                                      
  25 :           DIMENSION CG (5)                                                      
  26 :           DIMENSION CG0(5)                                                      
  27 :           DIMENSION XG (5)                                                      
  28 :           DIMENSION XR (2)                                                      
  29 : *                                                                               
  30 : * 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R                               
  31 : *                                                                               
  32 :          R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5             
  33 :          XR(1) = (QI(1) - P0(1)) / R                                            
  34 :          XR(2) = (QI(2) - P0(2)) / R                                            
  35 : *                                                                               
  36 : *                                                                               
  37 : *                                                                               
  38 : *                                                                               
  39 : * 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R                      
  40 : *                                                                               
  40 :        CALL ELPDGF(R,CB,CG0)                                                    
  40 :        CALL ELPGFS(R,XG)                                                        
  41 :         DO 10 I=1,5                                                             
  42 :           CG(I) = CG0(I) - XG(I)                                                
  43 : 10      CONTINUE                                                                
  44 : *                                                                               
  45 : * 3-   CALCUL DE V                                                              
  46 : *                                                                               
  47 :        COP(1) = CG(1)                                                           
  48 : *                                                                               
  49 : * 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)                                     
  50 : *                                                                               
  50 :        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,2,CAUX1)                                 
  50 :        CALL ELPDDR(CG,R,XR,XT,XT,XN,XN,2,CAUX2)                                 
  51 :        COP(3) = -1* XD* ( CAUX1 +  XNU*CAUX2)                                   
  52 : *                                                                               
  53 : * 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)                             
  54 : * ON NE LE CALCULE PAS
  54 : *       CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,3,CAUX1)                                
  54 : *       CALL ELPDDR(CG,R,XR,XN,XT,XT,XN,3,CAUX2)                                
  55 : *       COP(4) = -1* XD* ( CAUX1 +  (2. - XNU) *CAUX2)                          
  56 : *****                                                                           
  57 : *      TERMES DERIVES                                                           
  58 : *****                                                                           
  59 : *                                                                               
  60 : * 8-   CALCUL DE V,NN0                                                          
  61 : *                                                                               
  61 :        CALL ELPDDR(CG,R,XR,XN,XN0,XN,XN,2,CAUX1)                                
  62 :        COP(6) =  - CAUX1                                                        
  63 : *                                                                               
  64 : * 10-  CALCUL DE KN,N0                                                          
  65 : *                                                                               
  65 :        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN0,4,CAUX1)                                
  65 :        CALL ELPDDR(CG,R,XR,XN,XT,XT,XN0,4,CAUX2)                                
  66 :        COP(8) = -1* XD* (-  CAUX1 -  (2. - XNU) *CAUX2)                         
  67 :                                                                                 
  68 :        RETURN                                                                   
  69 :       END                                                                       
  70 :