Esope source files

   1 : C TRJARF    SOURCE    CHAT      05/01/13    03:48:35     5004           
   2 :       SUBROUTINE TRJARF(NAR,NF1,NF2,ITP)
   3 : C                                                                               
   4 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   5 : C                                                                               
   6 : C    CHERCHE LES NUMEROS DE FACES QUI CORRESPONDENT A UNE ARETE
   7 : C    DANS L ELEMENT DE REFERENCE
   8 : C                                                                               
   9 : C    E   NAR NUMERO DE L'ARETE
  10 : C    S   NF1 NF2 NUMEROS DES FACES ADJACENTES A CETTE ARETE
  11 : C    E   ITP TYPE DE L'ELEMENT
  12 : C                                                                               
  13 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  14 : C                                                                               
  15 : C                                                                               
  16 :       IMPLICIT INTEGER(I-N)
  17 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  18 : C                                                                               
  19 : C                                                                               
  20 : C                                                                               
  21 :       DIMENSION NFC8(2,12),NFP6(2,9), NFT4(2,6)
  22 :       DATA NFC8/6,1, 1,4, 4,2, 2,6, 3,1, 1,5, 5,2, 2,3,
  23 :      * 6,3, 3,4, 4,5, 5,6 /
  24 :       DATA NFP6/1,5, 5,2, 1,3, 3,2, 1,4, 4,2, 5,3, 3,4, 4,5 /
  25 :       DATA NFT4/1,2, 2,4, 4,1, 2,3, 3,4, 1,3/
  26 : C                                                                               
  27 : C                                                                               
  28 : C***  PRISMES                                                                   
  29 : C                                                                               
  30 :             IF(ITP.EQ.16)THEN
  31 :                  NF1=NFP6(1,NAR)
  32 :                  NF2=NFP6(2,NAR)
  33 : C                                                                               
  34 : C***  CUBES                                                                     
  35 : C                                                                               
  36 :             ELSEIF(ITP.EQ.14)THEN
  37 :                  NF1=NFC8(1,NAR)
  38 :                  NF2=NFC8(2,NAR)
  39 : C                                                                               
  40 : C***  TETRAEDRE
  41 : C                                                                               
  42 :             ELSEIF(ITP.EQ.23)THEN
  43 :                  NF1=NFT4(1,NAR)
  44 :                  NF2=NFT4(2,NAR)
  45 :           ELSE
  45 :               CALL ERREUR(16)
  46 :           ENDIF
  47 : C                                                                               
  48 :       RETURN
  49 :       END                                                                       
  50 :  
  51 :