Esope source files

   1 : C TRJCN5    SOURCE    CHAT      07/10/29    21:29:58     5940           
   2 :       SUBROUTINE TRJCN5(ITYP,IZSH,IJK)
   3 : C                                                                               
   4 : C                                                                               
   5 : C                                                                               
   6 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC           
   7 : C
   8 : C    Controle l'orientation des elements utilisés pour les 
   9 : C           trajectoires ISENS=1 sens direct
  10 : C                        ISENS=-1 sens inverse
  11 : C     on arete le calcul si tous les elements ne sont pas de meme sens
  12 : C     ITYP TYPE DE L ELEMENT
  13 : C     SHP FONCTION DE FORME
  14 : C     XYZL COORDONNEES REELLES DES NOEUDS DE L ELEMENT
  15 : C     XR COORDONNEE DU CENTRE DE GRAVITE DANS L ELEMENT DE REFERENCE
  16 : C                                                                               
  17 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC          
  18 : C                                                                               
  19 : C                                                                               
  20 :       IMPLICIT INTEGER(I-N)
  21 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  22 : C
  23 : C                                                                               
  24 :       DIMENSION XR(3)
  25 : -INC CCOPTIO
  26 : -INC CCGEOME
  27 :       SEGMENT IZSH
  28 :             REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  29 :       ENDSEGMENT
  30 : C                                                                               
  31 : C     SAVE ISENS
  32 : C     DATA ISENS/0/
  33 : C                                                                               
  34 :  
  35 :       IF(ITYP.EQ.4)THEN
  36 :            XR(1)=1.D0/3.D0
  37 :            XR(2)=1.D0/3.D0
  38 :            XR(3)=0.D0
  39 :            NBNO=3
  40 :       ELSEIF(ITYP.EQ.8)THEN
  41 :            XR(1)=0.D0
  42 :            XR(2)=0.D0
  43 :            XR(3)=0.D0
  44 :            NBNO=4
  45 :       ELSEIF(ITYP.EQ.14)THEN
  46 :            XR(1)=0.D0
  47 :            XR(2)=0.D0
  48 :            XR(3)=0.D0
  49 :            NBNO=8
  50 :       ELSEIF(ITYP.EQ.16)THEN
  51 :            XR(1)=1.D0/3.D0
  52 :            XR(2)=1.D0/3.D0
  53 :            XR(3)=0.D0
  54 :            NBNO=6
  55 :       ELSEIF(ITYP.EQ.23)THEN
  56 :            XR(1)=1.D0/3.D0
  57 :            XR(2)=1.D0/3.D0
  58 :            XR(3)=1.D0/3.D0
  59 :            NBNO=4
  60 :       ELSE
  61 :            MOTERR(1:4)=NOMS(ITYP)(1:4)
  62 :            MOTERR(5:12)='TRJCN5  '
  62 :            CALL ERREUR(86)
  63 :            RETURN
  64 :       ENDIF
  64 :       CALL SHAPE(XR(1),XR(2),XR(3),ITYP,SHP,IRET)
  65 :        MNO9=SHP(/2)
  65 :       CALL JACOBI(XYZL,SHP,IDIM,NBNO,DJAC)
  66 :       IJK=1
  67 :       IF(DJAC.LT.0.D0)IJK=-1
  68 : C     IF(IJK.EQ.ISENS)RETURN
  69 : C     IF(ISENS.EQ.0)THEN
  70 : C          ISENS=IJK
  71 : C          RETURN
  72 : C     ENDIF
  72 : C     CALL ERREUR(318)
  73 :       RETURN                                                                    
  74 :       END                                                                       
  75 :  
  76 :  
  77 :  
  78 :  
  79 :  
  80 :  
  81 :  
  82 :  
  83 :  
  84 :