Esope source files

   1 : C SESFR1    SOURCE    CHAT      05/01/13    03:14:17     5004           
   2 : C
   3 :       SUBROUTINE SESFR1(IT,IAR,ITRTRI,NBCMAX,IT2,IAR2)
   4 : C     *******************************************************
   5 : C     OBJET : ELEMENT SUIVANT SUR FRONTIERE IDE-1
   6 : C        TRIANGLE SUIVANT SUR ARETE / TETRA SUIVANT SUR FACE
   7 : C     EN ENTREE:
   8 : C      IT     :  LE TRIANGLE INITIAL
   9 : C      IAR    :  L'ARETE INITIALE DE IT
  10 : C      NBCMAX :  NOMBRE DE COTES MAXIMUM DES ELEMENTS DU MAILLAGE
  11 : C      ITRTRI:  TABLEAU DES VOISINS
  12 : C     EN SORTIE: 
  13 : C      IT1 :  LE TRIANGLE SUIVANT
  14 : C      IAR1:  INDICE DE L'ARETE IAR POUR LE TRIANGLE IT1
  15 : C     CONDITION D'APPLICATION : TRIANGLE, QUADRANGLE, TETRA   
  16 : C     *******************************************************
  17 :       IMPLICIT INTEGER(I-N)
  18 :       INTEGER   IT,IAR,NBCMAX,ITRTRI(*),IT2,IAR2
  19 : C
  20 :       INTEGER J,IT1,IAR1,IT0
  21 : C      
  22 :       IAR1 = 0
  23 :       IT1 = ITRTRI((IT-1)*NBCMAX+IAR)
  24 :       IF( IT1 .EQ. 0 )GO TO 20
  25 :       IF( IT1 .LT. 0 )IT1 = -IT1
  26 : C     --- RECHERCHE DE L'ARETE IAR1 DE IT1 ---
  27 :       DO 10 J=1,NBCMAX
  28 :         IT0 = ITRTRI((IT1-1)*NBCMAX + J)
  29 :         IF( IT0.LT. 0 )IT0 = -IT0
  30 :         IF( IT0 .EQ. IT )THEN
  31 :           IAR1 = J
  32 :           GO TO 20
  33 :         ENDIF
  34 :    10 CONTINUE
  35 :    20 IT2  = IT1
  36 :       IAR2 = IAR1
  37 :   999 END        
  38 :  
  39 :