1 : C CALUHH SOURCE CHAT 05/01/12 21:48:09 5004 2 : SUBROUTINE CALUHH(X,Y,PG,NG) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 6 : C CALCUL DES COORDONNEES ET DES POIDS DES POINTS DE HAMMER SUR LE 7 : C TRIANGLE DE REFERENCE 8 : C 9 : C NG NOMBRE DE POINTS D'INTEGRATION 10 : C 11 : C X,Y COORDONNEES DES POINTS DE HAMMER SUR LE TRIANGLE DE REFERENCE 12 : C DE SOMMETS (0.,0.) , (SQRT(2),0) , (0,SQRT(2)) 13 : C PG POIDS DES POINTS DE HAMMER 14 : C 15 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 16 : REAL*8 X(NG),Y(NG) 17 : DIMENSION PG(NG) 18 : REAL*8 A1,A2,B1,B2,C1,P1,P2,P3,R2 19 : C*** 20 : R2=SQRT(2.D0) 21 : 22 : IF(NG.NE.7) GO TO 10 23 : A1=(9.D0-2.D0*SQRT(15.D0))/21.D0 24 : A2=(9.D0+2.D0*SQRT(15.D0))/21.D0 25 : B1=(6.D0+SQRT(15.D0))/21.D0 26 : B2=(6.D0-SQRT(15.D0))/21.D0 27 : C1=1.D0/3.D0 28 : C 29 : P1=9.D0/40.D0 30 : P2=(155.D0+SQRT(15.D0))/1200.D0 31 : P3=(155.D0-SQRT(15.D0))/1200.D0 32 : C 33 : X(1)=C1 34 : Y(1)=C1 35 : X(2)=B1 36 : Y(2)=B1 37 : X(3)=A1 38 : Y(3)=B1 39 : X(4)=B1 40 : Y(4)=A1 41 : X(5)=B2 42 : Y(5)=B2 43 : X(6)=A2 44 : Y(6)=B2 45 : X(7)=B2 46 : Y(7)=A2 47 : DO 71 L=1,NG 48 : X(L)=X(L)*R2 49 : Y(L)=Y(L)*R2 50 : 71 CONTINUE 51 : C 52 : PG(1)=P1 53 : PG(2)=P2 54 : PG(3)=P2 55 : PG(4)=P2 56 : PG(5)=P3 57 : PG(6)=P3 58 : PG(7)=P3 59 : C 60 : RETURN 61 : 10 CONTINUE 62 : IF(NG.NE.1) GO TO 30 63 : X(1)=1.D0/3.D0*R2 64 : Y(1)=1.D0/3.D0*R2 65 : PG(1)=1.D0 66 : RETURN 67 : 30 CONTINUE 67 : IF(NG.NE.3) CALL ARRET(0) 68 : X(1)=0.D0 69 : Y(1)=0.5D0 70 : X(2)=0.5D0 71 : Y(2)=0.5D0 72 : X(3)=0.5D0 73 : Y(3)=0.D0 74 : DO 31 L=1,NG 75 : X(L)=X(L)*R2 76 : Y(L)=Y(L)*R2 77 : 31 CONTINUE 78 : PG(1)=1.D0/3.D0 79 : PG(2)=1.D0/3.D0 80 : PG(3)=1.D0/3.D0 81 : RETURN 82 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 83 : END 84 : 85 : 86 :
© Cast3M 2003 - All rights reserved.
Disclaimer