C T2IBT SOURCE CHAT 05/01/13 03:30:04 5004 SUBROUTINE T2IBT(BOITE,ITRI,NBE,COORD,NBN) C ********************************************************************** C OBJET : TRIANGULATION D'UN CARRE CONTENANT LA BOITE C EN ENTREE : C BOITE : LA BOITE QUE DOIT CONTENIR LE MAILLAGE C EN SORTIE : C ITRI : LA DEFINITION DES TRIANGLES (LEURS POINTS) C NBE : LE NOMBRE D'ELEMENTS C COORD : TABLEAU DES COORDONNEES C NBN : NOMBRE DE POINTS DE LA TRIANGULATION INITIALE C ********************************************************************** IMPLICIT INTEGER(I-N) REAL*8 BOITE(4),COORD(*) INTEGER NBN,ITRI(*),NBE C REAL*8 R, CTR(2), COEF INTEGER I C NBN = -1 NBE = -1 C R = 0.0D0 DO 10 I=1,2 CTR(I) = (BOITE(2+I) + BOITE(I)) / 2.D0 R = MAX(R,(BOITE(2+I) - BOITE(I))) 10 CONTINUE C COEF = SQRT(2.) * 50. COEF = 2.D0 C COORD(1) = CTR(1) - R * COEF COORD(2) = CTR(2) - R * COEF C COORD(3) = CTR(1) + R * COEF COORD(4) = CTR(2) - R * COEF C COORD(5) = CTR(1) + R * COEF COORD(6) = CTR(2) + R * COEF C COORD(7) = CTR(1) - R * COEF COORD(8) = CTR(2) + R * COEF C ITRI(1) = 1 ITRI(2) = 2 ITRI(3) = 3 ITRI(4) = 3 ITRI(5) = 4 ITRI(6) = 1 C NBN = 4 NBE = 2 C END