Esope source files

   1 : C T2IBT     SOURCE    CHAT      05/01/13    03:30:04     5004           
   2 :       SUBROUTINE T2IBT(BOITE,ITRI,NBE,COORD,NBN)
   3 : C     **********************************************************************
   4 : C     OBJET : TRIANGULATION D'UN CARRE CONTENANT LA BOITE
   5 : C     EN ENTREE  :
   6 : C        BOITE   :  LA BOITE QUE DOIT CONTENIR LE MAILLAGE
   7 : C     EN SORTIE  :
   8 : C        ITRI    :  LA DEFINITION DES TRIANGLES (LEURS POINTS)
   9 : C        NBE     :  LE NOMBRE D'ELEMENTS
  10 : C        COORD   :  TABLEAU DES COORDONNEES
  11 : C        NBN     :  NOMBRE DE POINTS DE LA TRIANGULATION INITIALE
  12 : C     **********************************************************************
  13 :       IMPLICIT INTEGER(I-N)
  14 :       REAL*8       BOITE(4),COORD(*)
  15 :       INTEGER    NBN,ITRI(*),NBE
  16 : C
  17 :       REAL*8     R, CTR(2), COEF
  18 :       INTEGER I
  19 : C
  20 :       NBN = -1
  21 :       NBE = -1
  22 : C
  23 :       R = 0.0D0
  24 :       DO 10 I=1,2
  25 :         CTR(I) = (BOITE(2+I) + BOITE(I)) / 2.D0
  26 :         R = MAX(R,(BOITE(2+I) - BOITE(I)))
  27 :    10 CONTINUE
  28 : C      COEF = SQRT(2.) * 50.
  29 :       COEF = 2.D0
  30 : C
  31 :       COORD(1) = CTR(1) - R * COEF
  32 :       COORD(2) = CTR(2) - R * COEF
  33 : C
  34 :       COORD(3) = CTR(1) + R * COEF
  35 :       COORD(4) = CTR(2) - R * COEF
  36 : C
  37 :       COORD(5) = CTR(1) + R * COEF
  38 :       COORD(6) = CTR(2) + R * COEF
  39 : C
  40 :       COORD(7) = CTR(1) - R * COEF
  41 :       COORD(8) = CTR(2) + R * COEF
  42 : C
  43 :       ITRI(1) = 1
  44 :       ITRI(2) = 2
  45 :       ITRI(3) = 3
  46 :       ITRI(4) = 3
  47 :       ITRI(5) = 4
  48 :       ITRI(6) = 1
  49 : C
  50 :       NBN = 4
  51 :       NBE = 2
  52 : C
  53 :       END
  54 :  
  55 :  
  56 :