C TRJVPO    SOURCE    CHAT      05/01/13    03:52:04     5004
      SUBROUTINE TRJVPO(UNOEU,UE,XY,ITYP,IDIM,NOEL1,SHP)
C*******************************************************************
C      CALCUL DE LA VITESSE AU POINT XY
C
C********************************************************************
C   ENTREES
C     UNOEU VITESSES AUX NOEUDS DE L ELEMENT
C     XY COORDONNEES DE REFERENCES DU POINT DE CALCUL
C     ITYP TYPE DE L ELEMENT
C     NOEL1  NOMBRE DE NOEUDS DE L ELEMENT POSITION
C   SORTIES
C     UE VITESSE DANS L ELEMENT DE REFERENCE POUR LE POINT
C        DE COORDONNEES XY
C
C     SHP(1,I) FONCTIONS DE FORMES AU POINT DE CALCUL
C     SHP(2,I) DERIVEE DES FONCTIONS DE FORME / XSI AU POINT DE CALCUL
C     SHP(3,I)    "     "      "     "    "   / ETA "    "   "    "
C     SHP(4,I)    "     "      "      "   "   /DZETA "   "    "   "
C ********************************************************************
C
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C
C
      DIMENSION UNOEU(*),XY(*),UE(*)
      DIMENSION SHP(6,*)
C
C
C***  CALCUL DES XI ET DES DXI
C
      CALL SHAPE(XY(1),XY(2),XY(3),ITYP,SHP,IRET)
C
C***  VITESSE
C
      DO 51 I=1,IDIM
      UE(I)=0.D0
      DO 52 K=1,NOEL1
      J=(K-1)*IDIM+I
      UE(I)=UE(I)+SHP(1,K)*UNOEU(J)
  52  CONTINUE
  51  CONTINUE
C     write(6,*) 'vitesse',(ue(i),i=1,idim)
C
  50  CONTINUE
C
      RETURN
      END