1 : C CALG2 SOURCE CHAT 05/01/12 21:46:02 5004 2 : SUBROUTINE CALG2(A,B,C,D,NG,H,U,X,Y,PG) 3 : C*********************************************************************** 4 : C CALCUL DES COORDONNEES ET DES POIDS DES POINTS DE GAUSS SUR LE PAVE 5 : C (A,B) X (C,D) LES CONNAISSANT SUR LE SEGMENT (-1,1) 6 : C 7 : C NG NOMBRE DE POINTS DE GAUSS SUR LE SEGMENT (-1,1) 8 : C U(NG) ABCISSES DES POINTS DE GAUSS 9 : C H(NG) POIDS DES POINTS DE GAUSS 10 : C 11 : C X,Y COORDONNEES DES POINTS DE GAUSS SUR LE PAVE (A,B) X (C,D) 12 : C PG POIDS DES POINTS DE GAUSS SUR LE PAVE (A,B) X (C,D) 13 : C 14 : C*********************************************************************** 15 : IMPLICIT INTEGER(I-N) 16 : IMPLICIT REAL*8 (A-H,O-Z) 17 : REAL*8 A,B,C,D,H(NG),U(NG),X(NG,NG),Y(NG,NG) 18 : REAL*8 PG(NG,NG) 19 : 20 : C*** 21 : 22 : BMA=(B-A)/2.D0 23 : BPA=(B+A)/2.D0 24 : DMC=(D-C)/2.D0 25 : DPC=(D+C)/2.D0 26 : C 27 : DO 1 LY=1,NG 28 : DO 1 LX=1,NG 29 : X(LX,LY)=BMA*U(LX)+BPA 30 : Y(LX,LY)=DMC*U(LY)+DPC 31 : PG(LX,LY)=BMA*DMC*H(LX)*H(LY) 32 : 1 CONTINUE 33 : 34 : C WRITE(6,*)' SUB CALG2 : POIDS ET POSITIONS DES PTS DE GAUSS ' 35 : C WRITE(6,*)' SUB CALG2 : H,U,X,Y ET PG ' 36 : C WRITE(6,1002)H 37 : C WRITE(6,1002)U 38 : C WRITE(6,1002)X 39 : C WRITE(6,1002)Y 40 : C WRITE(6,1002)PG 41 : 1002 FORMAT(10(1X,1PD11.4)) 42 : RETURN 43 : END 44 : 45 :
© Cast3M 2003 - All rights reserved.
Disclaimer