Télécharger la plaquette Cast3M


fr en

Télécharger s3faso.eso

Retour à la liste

Numérotation des lignes : 

s3faso

  1. C S3FASO    SOURCE    CHAT      06/03/29    21:32:51     5360
  2. C
  3.       FUNCTION S3FASO(I,N,ISENS,ISOM)
  4. C     ************************************************************
  5. C     OBJET : K SOMMETS DE LA FACE (INDICE RELATIF)
  6. C     EN ENTREE:
  7. C      I     :  LES INDICES DES FACES DE L'ELEMENT
  8. C      N     :  (4) NOMBRE DE NOEUDS DE L'ELEMENT
  9. C                 TETRA(4),PYRAM(5),PRISME(6),HEXA(8)
  10. C     EN SORTIE:
  11. C      ISOM  :  INDICE DES NOEUDS DE LA FACE (SENS DIRECT)
  12. C     CONDITION D'APPLICATION : TETRAEDRE ET HEXAEDRE SEULEMENT
  13. C     ************************************************************
  14.       IMPLICIT INTEGER(I-N)
  15.       INTEGER S3FASO
  16.       INTEGER I,N,ISENS,ISOM(*)
  17. C
  18.       COMMON /STRTET/ ITA2F(4,4),IT2FA(4,4),IT3SF(3,4)
  19.       INTEGER ITA2F, IT2FA, IT3SF
  20.       COMMON /STRHEX/ IQ4SF(4,6),IQ3FS(3,8)
  21.       INTEGER IQ4SF,IQ3FS
  22.       INTEGER J,iarr
  23. C
  24.       GOTO (1,1,1,100,1,1,1,200) N
  25. C        =====================
  26. C     --- ELEMENT NON RECONNU ---
  27. C        =====================
  28.     1 S3FASO = 0
  29.       iarr = -1
  30.       CALL DSERRE(1,iarr,'S3FASO',' TYPE D ELEMENT INCONNU')
  31.       GOTO 9999
  32. C        ==================
  33. C     --- CAS DU TETRAEDRE ---
  34. C        ==================
  35.   100 S3FASO = 3
  36.       IF( ISENS .EQ. 1 )THEN
  37.           DO 110 J=1,S3FASO
  38.             ISOM(J) = IT3SF(J,I)
  39.   110     CONTINUE
  40.         ELSE
  41.           DO 120 J=1,S3FASO
  42.             ISOM(J) = IT3SF(S3FASO+1-J,I)
  43.   120     CONTINUE
  44.         ENDIF
  45.         GOTO 9999
  46. C        ==================
  47. C     --- CAS DE L'HEXAEDRE ---
  48. C        ==================
  49.   200 S3FASO = 4
  50.       IF( ISENS .EQ. 1 )THEN
  51.           DO 210 J=1,S3FASO
  52.             ISOM(J) = IQ4SF(J,I)
  53.   210     CONTINUE
  54.         ELSE
  55.           DO 220 J=1,S3FASO
  56.             ISOM(J) = IQ4SF(S3FASO+1-J,I)
  57.   220     CONTINUE
  58.         ENDIF
  59.         GOTO 9999
  60. C
  61.  9999 END
  62.  
  63.  
  64.  
  65.  

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