Esope source files

   1 : C ELPDOP    SOURCE    CHAT      05/01/12    23:37:11     5004           
   2 :       SUBROUTINE ELPDOP(P0,XN0,QI,XN,XT,XD,XNU,COP,CB,ISTAT)                    
   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 : *       POUR UN COUPLE DE PTS P0 , QI                                           
  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 XG (5)                                                      
  27 :           DIMENSION XR (2)                                                      
  28 : *                                                                               
  29 : * 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R                               
  30 : *                                                                               
  31 :          R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5             
  32 :          XR(1) = (QI(1) - P0(1)) / R                                            
  33 :          XR(2) = (QI(2) - P0(2)) / R                                            
  34 : *                                                                               
  35 : *                                                                               
  36 : *                                                                               
  37 : *                                                                               
  38 : * 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R                      
  39 : *                                                                               
  40 :        IF ( ISTAT.EQ.1) THEN                                                    
  40 :          CALL ELPGFS(R,XG)                                                      
  41 :          DO  10 I=1,5                                                           
  42 :            CG(I) = XG(I)*(1d0,0d0)
  43 :  10      CONTINUE                                                               
  44 :        ELSE                                                                     
  44 :          CALL ELPDGF(R,CB,CG)                                                   
  45 :        ENDIF                                                                    
  46 : *                                                                               
  47 : * 3-   CALCUL DE V                                                              
  48 : *                                                                               
  49 :        COP(1) = CG(1)                                                           
  50 : *                                                                               
  51 : * 4-   CALCUL DE V,N                                                            
  52 : *                                                                               
  52 :        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,1,CAUX1)                                 
  53 :        COP(2) = CAUX1                                                           
  54 : *                                                                               
  55 : * 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)                                     
  56 : *                                                                               
  56 :        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,2,CAUX1)                                 
  56 :        CALL ELPDDR(CG,R,XR,XT,XT,XN,XN,2,CAUX2)                                 
  57 :        COP(3) = -1d0* XD* ( CAUX1 +  XNU*CAUX2)
  58 : *                                                                               
  59 : * 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)                             
  60 : *                                                                               
  60 :        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,3,CAUX1)                                 
  60 :        CALL ELPDDR(CG,R,XR,XN,XT,XT,XN,3,CAUX2)                                 
  61 :        COP(4) = -1d0* XD* ( CAUX1 +  (2d0 - XNU) *CAUX2)
  62 : *****                                                                           
  63 : *      TERMES DERIVES                                                           
  64 : *****                                                                           
  65 : * 7-   CALCUL DE V,N0                                                           
  66 : *                                                                               
  67 : *                                                                               
  68 : *                                                                               
  69 : *                                                                               
  69 :        CALL ELPDDR(CG,R,XR,XN0,XN,XN,XN,1,CAUX1)                                
  70 :        COP(5) = - CAUX1                                                         
  71 : *                                                                               
  72 : * 8-   CALCUL DE V,NN0                                                          
  73 : *                                                                               
  73 :        CALL ELPDDR(CG,R,XR,XN,XN0,XN,XN,2,CAUX1)                                
  74 :        COP(6) =  - CAUX1                                                        
  75 : *                                                                               
  76 : * 9-   CALCUL DE MN,N0                                                          
  77 : *                                                                               
  77 :        CALL ELPDDR(CG,R,XR,XN,XN,XN0,XN,3,CAUX1)                                
  77 :        CALL ELPDDR(CG,R,XR,XT,XT,XN0,XN,3,CAUX2)                                
  78 :        COP(7) = - 1d0* XD* ( - CAUX1 -  XNU*CAUX2)
  79 : *                                                                               
  80 : * 9-   CALCUL DE KN,N0                                                          
  81 : *                                                                               
  81 :        CALL ELPDDR(CG,R,XR,XN,XN,XN,XN0,4,CAUX1)                                
  81 :        CALL ELPDDR(CG,R,XR,XN,XT,XT,XN0,4,CAUX2)                                
  82 :        COP(8) = -1d0* XD* (-  CAUX1 -  (2d0 - XNU) *CAUX2)
  83 :                                                                                 
  84 :        RETURN                                                                   
  85 :       END                                                                       
  86 :