Sources Esope | Cast3M

caluhh
C CALUHH    SOURCE    CHAT      05/01/12    21:48:09     5004
      SUBROUTINE CALUHH(X,Y,PG,NG)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
C   CALCUL DES COORDONNEES ET DES POIDS DES POINTS DE HAMMER SUR LE
C   TRIANGLE DE REFERENCE
C
C   NG     NOMBRE DE POINTS D'INTEGRATION
C
C   X,Y   COORDONNEES DES POINTS DE HAMMER SUR LE TRIANGLE DE REFERENCE
C         DE SOMMETS (0.,0.) , (SQRT(2),0) , (0,SQRT(2))
C   PG     POIDS DES POINTS DE HAMMER
C
C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
      REAL*8 X(NG),Y(NG)
      DIMENSION  PG(NG)
      REAL*8 A1,A2,B1,B2,C1,P1,P2,P3,R2
C***
      R2=SQRT(2.D0)
 
      IF(NG.NE.7) GO TO 10
      A1=(9.D0-2.D0*SQRT(15.D0))/21.D0
      A2=(9.D0+2.D0*SQRT(15.D0))/21.D0
      B1=(6.D0+SQRT(15.D0))/21.D0
      B2=(6.D0-SQRT(15.D0))/21.D0
      C1=1.D0/3.D0
C
      P1=9.D0/40.D0
      P2=(155.D0+SQRT(15.D0))/1200.D0
      P3=(155.D0-SQRT(15.D0))/1200.D0
C
      X(1)=C1
      Y(1)=C1
      X(2)=B1
      Y(2)=B1
      X(3)=A1
      Y(3)=B1
      X(4)=B1
      Y(4)=A1
      X(5)=B2
      Y(5)=B2
      X(6)=A2
      Y(6)=B2
      X(7)=B2
      Y(7)=A2
      DO 71 L=1,NG
      X(L)=X(L)*R2
      Y(L)=Y(L)*R2
   71 CONTINUE
C
      PG(1)=P1
      PG(2)=P2
      PG(3)=P2
      PG(4)=P2
      PG(5)=P3
      PG(6)=P3
      PG(7)=P3
C
      RETURN
   10 CONTINUE
      IF(NG.NE.1) GO TO 30
      X(1)=1.D0/3.D0*R2
      Y(1)=1.D0/3.D0*R2
      PG(1)=1.D0
      RETURN
   30 CONTINUE
      IF(NG.NE.3) CALL ARRET(0)
      X(1)=0.D0
      Y(1)=0.5D0
      X(2)=0.5D0
      Y(2)=0.5D0
      X(3)=0.5D0
      Y(3)=0.D0
      DO 31 L=1,NG
      X(L)=X(L)*R2
      Y(L)=Y(L)*R2
   31 CONTINUE
      PG(1)=1.D0/3.D0
      PG(2)=1.D0/3.D0
      PG(3)=1.D0/3.D0
      RETURN
C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
      END