C ELPOP1    SOURCE    CHAT      05/01/12    23:37:51     5004
      SUBROUTINE ELPOP1(P0,XN0,QI,XN,XT,XD,XNU,XOP)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
************************************************************************
*
*    PLAQUES EN STATIQUE PAR EQUATION INTEGRALE:
*
*       CALCUL DES OPERATEURS PLAQUE ET DE LEURS DERIVEES
*       POUR UN COUPLE DE PTS P0 , QI
*
************************************************************************
*
*
          DIMENSION P0(2)
          DIMENSION QI(2)
          DIMENSION XOP(8)
*
          DIMENSION XN (2)
          DIMENSION XN0(2)
*
          DIMENSION XT (2)
          DIMENSION XT0(2)
          DIMENSION XG (5)
          DIMENSION XR (2)
*
* 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R
*
         R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5
         XR(1) = (QI(1) - P0(1)) / R
         XR(2) = (QI(2) - P0(2)) / R
*
*
*
*
* 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R
*
       CALL ELPGFS(R,XG)
*
* 3-   CALCUL DE V
*
       XOP(1) = XG(1)
*
* 4-   CALCUL DE V,N
*
       CALL ELPDER(XG,R,XR,XN,XN,XN,XN,1,AUX1)
       XOP(2) = AUX1
*
* 5-   CALCUL DE MN= -D * (V,NN + NU* V,TT)
*
       CALL ELPDER(XG,R,XR,XN,XN,XN,XN,2,AUX1)
       CALL ELPDER(XG,R,XR,XT,XT,XN,XN,2,AUX2)
       XOP(3) = -1* XD* ( AUX1 +  XNU*AUX2)
*
* 6-   CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)
*
       CALL ELPDER(XG,R,XR,XN,XN,XN,XN,3,AUX1)
       CALL ELPDER(XG,R,XR,XN,XT,XT,XN,3,AUX2)
       XOP(4) = -1* XD* ( AUX1 +  (2. - XNU) *AUX2)
*****
*      TERMES DERIVES
*****
* 7-   CALCUL DE V,N0
*
*
*
*
       CALL ELPDER(XG,R,XR,XN0,XN,XN,XN,1,AUX1)
       XOP(5) = - AUX1
*
* 8-   CALCUL DE V,NN0
*
       CALL ELPDER(XG,R,XR,XN,XN0,XN,XN,2,AUX1)
       XOP(6) =  - AUX1
*
* 9-   CALCUL DE MN,N0
*
       CALL ELPDER(XG,R,XR,XN,XN,XN0,XN,3,AUX1)
       CALL ELPDER(XG,R,XR,XT,XT,XN0,XN,3,AUX2)
       XOP(7) = - 1* XD* ( - AUX1 -  XNU*AUX2)
*
* 9-   CALCUL DE KN,N0
*
       CALL ELPDER(XG,R,XR,XN,XN,XN,XN0,4,AUX1)
       CALL ELPDER(XG,R,XR,XN,XT,XT,XN0,4,AUX2)
       XOP(8) = -1* XD* (-  AUX1 -  (2. - XNU) *AUX2)

       RETURN
      END