Sources Esope | Cast3M

trjcn5
C TRJCN5    SOURCE    BP208322  16/11/18    21:21:48     9177           
      SUBROUTINE TRJCN5(ITYP,IZSH,IJK)
C
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C    Controle l'orientation des elements utilisés pour les
C           trajectoires ISENS=1 sens direct
C                        ISENS=-1 sens inverse
C     on arete le calcul si tous les elements ne sont pas de meme sens
C     ITYP TYPE DE L ELEMENT
C     SHP FONCTION DE FORME
C     XYZL COORDONNEES REELLES DES NOEUDS DE L ELEMENT
C     XR COORDONNEE DU CENTRE DE GRAVITE DANS L ELEMENT DE REFERENCE
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C
      DIMENSION XR(3)
 
-INC PPARAM
-INC CCOPTIO
-INC CCGEOME
      SEGMENT IZSH
            REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
      ENDSEGMENT
C
C     SAVE ISENS
C     DATA ISENS/0/
C
 
      IF(ITYP.EQ.4)THEN
           XR(1)=1.D0/3.D0
           XR(2)=1.D0/3.D0
           XR(3)=0.D0
           NBNO=3
      ELSEIF(ITYP.EQ.8)THEN
           XR(1)=0.D0
           XR(2)=0.D0
           XR(3)=0.D0
           NBNO=4
      ELSEIF(ITYP.EQ.14)THEN
           XR(1)=0.D0
           XR(2)=0.D0
           XR(3)=0.D0
           NBNO=8
      ELSEIF(ITYP.EQ.16)THEN
           XR(1)=1.D0/3.D0
           XR(2)=1.D0/3.D0
           XR(3)=0.D0
           NBNO=6
      ELSEIF(ITYP.EQ.23)THEN
           XR(1)=1.D0/3.D0
           XR(2)=1.D0/3.D0
           XR(3)=1.D0/3.D0
           NBNO=4
      ELSE
           MOTERR(1:4)=NOMS(ITYP)(1:4)
           MOTERR(5:12)='TRJCN5  '
           CALL ERREUR(86)
           RETURN
      ENDIF
      CALL SHAPE(XR(1),XR(2),XR(3),ITYP,SHP,IRET)
       MNO9=SHP(/2)
      CALL JACOBI(XYZL,SHP,IDIM,NBNO,DJAC)
      IJK=1
      IF(DJAC.LT.0.D0)IJK=-1
C     IF(IJK.EQ.ISENS)RETURN
C     IF(ISENS.EQ.0)THEN
C          ISENS=IJK
C          RETURN
C     ENDIF
C     CALL ERREUR(318)
      RETURN
      END