Esope source files

   1 : C COBAR3    SOURCE    CHAT      05/01/12    22:12:14     5004           
   2 :       SUBROUTINE COBAR3(X,XELE,TLI)
   3 : C                                                                               
   4 : C**********************************************************************         
   5 : C                                                                               
   6 : C   OBJET  : CE SOUS-PROGRAMME CALCULE LES COORDONNEES BARYCENTRIQUES
   7 : C   -----    DANS UN TRI3 AU POINT DONT LES COORDONNEES SONT DANS X
   8 : C                                                                               
   9 : C   ARGUMENTS:                                                                  
  10 : C   ---------                                                                   
  11 : C                                                                               
  12 : C   ENTREE  :
  13 : C              X      = TABLEAU DES COORDONNEES DU POINT                        
  14 : C              XELE   = TABLEAU DES COORDONNEES DU TRI3                         
  15 : C                                                                               
  16 : C   SORTIE  :
  17 : C              TLI    = TABLEAU DES COORDONNEES BARYCENTRIQUES
  18 : C                       DU POINT X
  19 : C                                                                               
  20 : C      ISSU DE FON302 DE TRIOEF
  21 : C         UTILISE PAR L OPERATEUR TRAJ
  22 : C***********************************************************************        
  23 : C                                                                               
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  26 :       DIMENSION X(2),XELE(2,*),NUDROI(3),NOEDRO(2,3)
  27 : C                                                                               
  28 : C                                                                               
  29 :       DIMENSION TLI(*)
  30 : C                                                                               
  31 :       DATA NOEDRO/2,3 ,3,1 ,1,2/
  32 : C                                                                               
  33 : C                                                                               
  34 :       DO 10 I=1,3                                                               
  35 :       NOE1=NOEDRO(1,I)
  36 :       NOE2=NOEDRO(2,I)
  37 :       COEFX =XELE(2,NOE1)-XELE(2,NOE2)
  38 :       COEFY=XELE(1,NOE2)-XELE(1,NOE1)
  39 :       CONS=XELE(1,NOE1)*XELE(2,NOE2)-XELE(1,NOE2)*XELE(2,NOE1)
  40 :       CONSNI=XELE(1,I)*COEFX +XELE(2,I)*COEFY+CONS
  41 :       DROITE=X(1)*COEFX +X(2)*COEFY+CONS
  42 :       TLI(I)=DROITE/CONSNI
  43 :    10 CONTINUE
  44 : C                                                                               
  45 : C                                                                               
  46 : C                                                                               
  47 : C                                                                               
  48 :       RETURN
  49 :       END                                                                       
  50 :