Numérotation des lignes :

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 LEC   TRIANGLE DE REFERENCECC   NG     NOMBRE DE POINTS D'INTEGRATIONCC   X,Y   COORDONNEES DES POINTS DE HAMMER SUR LE TRIANGLE DE REFERENCEC         DE SOMMETS (0.,0.) , (SQRT(2),0) , (0,SQRT(2))C   PG     POIDS DES POINTS DE HAMMERCC>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>      REAL*8 X(NG),Y(NG)      DIMENSION  PG(NG)      REAL*8 A1,A2,B1,B2,C1,P1,P2,P3,R2C***      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.D0C      P1=9.D0/40.D0      P2=(155.D0+SQRT(15.D0))/1200.D0      P3=(155.D0-SQRT(15.D0))/1200.D0C      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 CONTINUEC      PG(1)=P1      PG(2)=P2      PG(3)=P2      PG(4)=P2      PG(5)=P3      PG(6)=P3      PG(7)=P3C      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      RETURNC>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>      END    

© Cast3M 2003 - Tous droits réservés.
Mentions légales