1 : C GEOCST SOURCE CHAT 05/01/13 00:17:47 5004 2 : SUBROUTINE GEOCST(XEL,GEOM) 3 : C======================================================================= 4 : C 5 : C CALCULE DES CONSTANTES D UN TRIANGLE POUR LE DKT 6 : C 7 : C ENTREES 8 : C XEL(3,3)= COORDONNEES DE L ELEMENT 9 : C SORTIES 10 : C GEOM(17)=DES CONSTANTES POUR LE DKT 11 : C 12 : C======================================================================= 13 : IMPLICIT INTEGER(I-N) 14 : IMPLICIT REAL*8(A-H,O-Z) 15 : DIMENSION XEL(3,*),GEOM(*) 16 : DATA TROIS/3.D0/,SIX/6.D0/ 17 : X12=XEL(1,1)-XEL(1,2) 18 : X23=XEL(1,2)-XEL(1,3) 19 : X31=XEL(1,3)-XEL(1,1) 20 : Y12=XEL(2,1)-XEL(2,2) 21 : Y23=XEL(2,2)-XEL(2,3) 22 : Y31=XEL(2,3)-XEL(2,1) 23 : AIREDO=-Y31*X12+Y12*X31 24 : VJ11=Y31/AIREDO 25 : VJ12=Y12/AIREDO 26 : VJ21=-X31/AIREDO 27 : VJ22=-X12/AIREDO 28 : XSQ12=X12*X12 29 : XSQ23=X23*X23 30 : XSQ31=X31*X31 31 : YSQ12=Y12*Y12 32 : YSQ23=Y23*Y23 33 : YSQ31=Y31*Y31 34 : DIST12=XSQ12+YSQ12 35 : DIST23=XSQ23+YSQ23 36 : DIST31=XSQ31+YSQ31 37 : CL4=-SIX*X23/DIST23 38 : CL5=-SIX*X31/DIST31 39 : CL6=-SIX*X12/DIST12 40 : SL4=-SIX*Y23/DIST23 41 : SL5=-SIX*Y31/DIST31 42 : SL6=-SIX*Y12/DIST12 43 : D4=TROIS*X23*Y23/DIST23 44 : D5=TROIS*X31*Y31/DIST31 45 : D6=TROIS*X12*Y12/DIST12 46 : SU4=TROIS*YSQ23/DIST23 47 : SU5=TROIS*YSQ31/DIST31 48 : SU6=TROIS*YSQ12/DIST12 49 : GEOM(1)=CL4 50 : GEOM(2)=CL5 51 : GEOM(3)=CL6 52 : GEOM(4)=SL4 53 : GEOM(5)=SL5 54 : GEOM(6)=SL6 55 : GEOM(7)=D4 56 : GEOM(8)=D5 57 : GEOM(9)=D6 58 : GEOM(10)=SU4 59 : GEOM(11)=SU5 60 : GEOM(12)=SU6 61 : GEOM(13)=VJ11 62 : GEOM(14)=VJ12 63 : GEOM(15)=VJ21 64 : GEOM(16)=VJ22 65 : GEOM(17)=AIREDO 66 : RETURN 67 : END 68 :
© Cast3M 2003 - All rights reserved.
Disclaimer