1 : C CALG3 SOURCE CHAT 05/01/12 21:46:06 5004 2 : SUBROUTINE CALG3(A,B,C,D,E,F,NG,H,U,X,Y,Z,PG) 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 GAUSS SUR LE PAVE 7 : C (A,B) X (C,D) X (E,F) LES CONNAISSANT SUR LE SEGMENT (-1,1) 8 : C 9 : C NG NOMBRE DE POINTS DE GAUSS SUR LE SEGMENT (-1,1) 10 : C U(NG) ABCISSES DES POINTS DE GAUSS 11 : C H(NG) POIDS DES POINTS DE GAUSS 12 : C 13 : C X,Y,Z COORDONNEES DES POINTS DE GAUSS SUR LE PAVE (A,B)X(C,D)X(E,F) 14 : C PG POIDS DES POINTS DE GAUSS SUR LE PAVE (A,B)X(C,D)X(E,F) 15 : C 16 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 17 : 18 : REAL*8 A,B,C,D,E,F 19 : REAL*8 H(NG),U(NG) 20 : REAL*8 X(NG,NG,NG),Y(NG,NG,NG),Z(NG,NG,NG) 21 : DIMENSION PG(NG,NG,NG) 22 : 23 : 24 : C*** 25 : BMA=(B-A)/2.D0 26 : BPA=(B+A)/2.D0 27 : DMC=(D-C)/2.D0 28 : DPC=(D+C)/2.D0 29 : FME=(F-E)/2.D0 30 : FPE=(F+E)/2.D0 31 : C 32 : DO 1 LZ=1,NG 33 : DO 1 LY=1,NG 34 : DO 1 LX=1,NG 35 : X(LX,LY,LZ)=BMA*U(LX)+BPA 36 : Y(LX,LY,LZ)=DMC*U(LY)+DPC 37 : Z(LX,LY,LZ)=FME*U(LZ)+FPE 38 : PG(LX,LY,LZ)=BMA*DMC*FME*H(LX)*H(LY)*H(LZ) 39 : 1 CONTINUE 40 : C WRITE(6,*)' SUB CALG3 : POIDS ET POSITIONS DES PTS DE GAUSS ' 41 : C WRITE(6,*)' SUB CALG3 : H,U,X,Y ET PG ' 42 : C WRITE(6,1002)H 43 : C WRITE(6,1002)U 44 : C WRITE(6,1002)X 45 : C WRITE(6,1002)Y 46 : C WRITE(6,1002)PG 47 : 1002 FORMAT(10(1X,1PD11.4)) 48 : RETURN 49 : END 50 : 51 :
© Cast3M 2003 - All rights reserved.
Disclaimer