Esope source files

   1 : C TRJVPO    SOURCE    CHAT      05/01/13    03:52:04     5004           
   2 :       SUBROUTINE TRJVPO(UNOEU,UE,XY,ITYP,IDIM,NOEL1,SHP)
   3 : C*******************************************************************
   4 : C      CALCUL DE LA VITESSE AU POINT XY
   5 : C                                                                               
   6 : C********************************************************************           
   7 : C   ENTREES
   8 : C     UNOEU VITESSES AUX NOEUDS DE L ELEMENT                                    
   9 : C     XY COORDONNEES DE REFERENCES DU POINT DE CALCUL                           
  10 : C     ITYP TYPE DE L ELEMENT                                                    
  11 : C     NOEL1  NOMBRE DE NOEUDS DE L ELEMENT POSITION                             
  12 : C   SORTIES
  13 : C     UE VITESSE DANS L ELEMENT DE REFERENCE POUR LE POINT
  14 : C        DE COORDONNEES XY
  15 : C
  16 : C     SHP(1,I) FONCTIONS DE FORMES AU POINT DE CALCUL
  17 : C     SHP(2,I) DERIVEE DES FONCTIONS DE FORME / XSI AU POINT DE CALCUL
  18 : C     SHP(3,I)    "     "      "     "    "   / ETA "    "   "    "
  19 : C     SHP(4,I)    "     "      "      "   "   /DZETA "   "    "   "
  20 : C ********************************************************************          
  21 : C                                                                               
  22 : C                                                                               
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  25 : C
  26 : C                                                                               
  27 : C                                                                               
  28 :       DIMENSION UNOEU(*),XY(*),UE(*)
  29 :       DIMENSION SHP(6,*)
  30 : C                                                                               
  31 : C                                                                               
  32 : C***  CALCUL DES XI ET DES DXI                                                  
  33 : C                                                                               
  33 :       CALL SHAPE(XY(1),XY(2),XY(3),ITYP,SHP,IRET)
  34 : C                                                                               
  35 : C***  VITESSE                                                                   
  36 : C                                                                               
  37 :       DO 51 I=1,IDIM                                                            
  38 :       UE(I)=0.D0
  39 :       DO 52 K=1,NOEL1                                                           
  40 :       J=(K-1)*IDIM+I
  41 :       UE(I)=UE(I)+SHP(1,K)*UNOEU(J)
  42 :   52  CONTINUE                                                                  
  43 :   51  CONTINUE                                                                  
  44 : C     write(6,*) 'vitesse',(ue(i),i=1,idim)
  45 : C                                                                               
  46 :   50  CONTINUE                                                                  
  47 : C                                                                               
  48 :       RETURN                                                                    
  49 :       END                                                                       
  50 :  
  51 :