Esope source files

   1 : C TRJVLO    SOURCE    CHAT      05/01/13    03:51:59     5004           
   2 :       SUBROUTINE TRJVLO(IZNOEU,IZUN,IZTRAV,ITYP,UELEM,NUPOIN,IEL,IZSH)
   3 : C                                                                               
   4 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC           
   5 : C                                                                               
   6 : C     CALCUL DE LA VITESSE (REELLE) UELEM A L'INTERIEUR D'UN ELEMENT
   7 : C     EN FONCTION DES VITESSES NODALES
   8 : C                                                                               
   9 : C     IZNOEU COORDONNEES ET NUMEROS DES NOEUDS DE L'ELEMENT
  10 : C     IZUN VITESSES AUX NOEUDS DE L'ELEMENT
  11 : C     IZTRAV COORDONNEES DES NOEUDS (ON S'INTERESSE AU NUPOIN-IEME )
  12 : C     ITYP TYPE DE L'ELEMENT
  13 : C     NUPOIN NUMERO DE LA PARTICULE A LAQUELLE ON S'INTERESSE
  14 : C     IEL NUMERO DE L ELEMENT
  15 : C     ITYP TYPE DE LA VITESSE
  16 : C     IZSH SEGMENT DE TRAVAIL POUR LES FONCTIONS DE FORME
  17 : C                                                                               
  18 : C                                                                               
  19 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC       
  20 : C                                                                               
  21 : C                                                                               
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  24 : C                                                                               
  25 : 
  26 :       SEGMENT  IZNOEU 
  27 :             REAL*8 XELE(IDIM,NOEL)  
  28 :             INTEGER NOEGLO(NOEL)
  29 :       ENDSEGMENT
  30 :       SEGMENT IZTRAV
  31 :            REAL*8  COOR(NDIM,NPART)
  32 :       ENDSEGMENT
  33 : C
  34 :       SEGMENT IZSH
  35 :             REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  36 :       ENDSEGMENT
  37 : C                                                                               
  38 :       SEGMENT IZUN
  39 :            REAL*8 UN(I1,I2,I3)
  40 :       ENDSEGMENT
  41 : C                                                                               
  42 :       DIMENSION UELEM(3),XYREF(3),XYZP(3)
  43 : C                                                                               
  44 :       IDIM=XELE(/1)                                                             
  45 :       NOEL=XELE(/2)                                                             
  46 : C                                                                               
  47 :       DO 1 ID=1,IDIM
  48 :         XYZP(ID)=COOR(ID,NUPOIN)
  49 :    1  CONTINUE
  50 : C                                                                               
  51 : C                                                                               
  52 :        IER=0
  53 : C                     TRI3
  53 :         IF(ITYP.EQ.4)CALL TRJ302(XELE,XYZP,XYREF)
  54 : C                     QUA4
  54 :         IF(ITYP.EQ.8)CALL TRJ402(XELE,XYZP,XYREF,IER)
  55 : C                     TRI6
  55 :         IF(ITYP.EQ.6)CALL TRJ602(XELE,XYZP,XYREF,IER)
  56 : C                     TRI7
  56 :         IF(ITYP.EQ.7)CALL TRJ702(XELE,XYZP,XYREF,IER)
  57 : C                     QUA9
  57 :         IF(ITYP.EQ.11)CALL TRJ902(XELE,XYZP,XYREF,IER)
  58 : C                     PRI6
  58 :         IF(ITYP.EQ.16)CALL TRJ603(XELE,XYZP,XYREF,IER)
  59 : C                     CUB8
  59 :         IF(ITYP.EQ.14)CALL TRJ803(XELE,XYZP,XYREF,IER)
  60 : C                     TET4
  60 :         IF(ITYP.EQ.23)CALL TRJ403(XELE,XYZP,XYREF)
  61 :       IF(IER.NE.0)THEN
  61 :                    CALL ERREUR(IER)
  62 :       ENDIF
  63 : C                                                                               
  64 : C                                                                               
  64 :         CALL TRJVPO(UN(1,1,IEL),UELEM,XYREF,ITYP,IDIM,NOEL,SHP)
  65 :       RETURN                                                                    
  66 :       END                                                                       
  67 :  
  68 :  
  69 :