Numérotation des lignes :

C NOEREF    SOURCE    CHAT      05/01/13    01:58:11     5004      SUBROUTINE NOEREF(XYZ,ITYP,INO)CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC     RECHERCHE SI LE POINT XYZ EST UN NOEUD DANS LES COORDONNEESC       DE REFERENCEC     INO=0  CE N EST PAS UN NOEUDC     INO=1  C'EST UN NOEUDCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8 (A-H,O-Z)C      DIMENSION XYZ(3)C      EPS=1.D-14C      NDIM=2C        ON CHERCHE SI ON A AFFAIRE    A UN NOEUD DU MAILLAGE       ITEST=0       INO=0      IF(ITYP.EQ.4.OR.ITYP.EQ.6.OR.ITYP.EQ.7)THENC              L ELEMENT DE REFERENCE EST UN TRIANGLE             DO 30 I=1,NDIM             IF(ABS(XYZ(I)).LT.EPS)ITEST=ITEST+1             IF(ABS(1.D0-XYZ(I)).LT.EPS)ITEST=ITEST+1   30        CONTINUE      ELSEIF(ITYP.EQ.8.OR.ITYP.EQ.11)THENC              L ELEMENT DE REFERENCE EST UN CARRE             DO 35 I=1,NDIM             IF(ABS(1.D0+XYZ(I)).LT.EPS)ITEST=ITEST+1             IF(ABS(1.D0-XYZ(I)).LT.EPS)ITEST=ITEST+1   35        CONTINUE      ELSEIF(ITYP.EQ.16)THENC              L ELEMENT DE REFERENCE EST UN PRISME             NDIM=3             IF((ABS(XYZ(3)-1.D0).LT.EPS).OR.     *                 (ABS(XYZ(3)+1.D0).LT.EPS))ITEST=ITEST+1             DO 40 I=1,2             IF(ABS(XYZ(I)).LT.EPS)ITEST=ITEST+1             IF(ABS(1.D0-XYZ(I)).LT.EPS)ITEST=ITEST+1   40        CONTINUE      ELSEIF(ITYP.EQ.14)THENC              L ELEMENT DE REFERENCE EST UN CUBE             NDIM=3             DO 45 I=1,NDIM             IF(ABS(1.D0+XYZ(I)).LT.EPS)ITEST=ITEST+1             IF(ABS(1.D0-XYZ(I)).LT.EPS)ITEST=ITEST+1   45        CONTINUE      ELSEIF(ITYP.EQ.23)THENC              L ELEMENT DE REFERENCE EST UN TETRAEDRE             NDIM=3             DO 50 I=1,NDIM             IF(ABS(XYZ(I)).LT.EPS)ITEST=ITEST+1             IF(ABS(1.D0-XYZ(I)).LT.EPS)ITEST=ITEST+1   50        CONTINUE      ENDIF      IF(ITEST.EQ.NDIM)INO=1C      RETURN      END

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