Esope source files

   1 : C VHYNOE    SOURCE    CHAT      05/01/13    04:07:09     5004           
   2 :       SUBROUTINE VHYNOE(UELEM,IZNOEU,IZUN,NDIM,XY,ITP,IEL)
   3 : C                                                                               
   4 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   5 : C                                                                               
   6 : C     CALCUL DE LA VITESSE REELE POUR UN NOEUD DU MAILLAGE EN
   7 : C     FORMULATION HYBRIDE
   8 : C                                                                               
   9 : C     ENTREES
  10 : C          IZNOEU POINTEUR 
  11 : C              XELE(NDIM,NOEL) COORDONNES DES NOEUDS DE L ELEMENT CONSIDERE
  12 : C          IZUN POINTEUR FLUX AUX CENTRES DES FACES DE L ELEMENT
  13 : C          ITP             TYPE DE L'ELEMENT CONSIDERE
  14 : C          XY          coordonnees DU NOEUD POUR LEQUEL ON CALCULERA
  15 : C                         LA VITESSE
  16 : C     SORTIE
  17 : C          UELEM VITESSE REELLE AU NOEUD INOELO
  18 : C                                                                               
  19 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  20 : C                                                                               
  21 : C                                                                               
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  24 : C                                                                               
  25 : C                                                                               
  26 : C                                                                               
  27 :       SEGMENT  IZNOEU
  28 :             REAL*8 XELE(IDIM,NOEL)  
  29 :             INTEGER NOEGLO(NOEL)
  30 :       ENDSEGMENT
  31 :       SEGMENT IZUN
  32 :            REAL*8 UN(I1,I2,I3)
  33 :       ENDSEGMENT
  34 :       SEGMENT HYTRAV
  35 :             REAL*8 XEL(3,NOEL),SHP(6,NOEL),SHY(12,NOEL)
  36 :             REAL*8 ZJAC(NDIM,NDIM)
  37 :       ENDSEGMENT
  38 : C                                                                               
  39 : C                                                                               
  40 :       DIMENSION XY(3),UELEM(*),VECT(3)
  41 : C                                                                               
  42 : C                                                                               
  43 : C                     FORMULATION HYBRIDE ON VA CALCULER LA VITESSE
  44 : C                   REELLE AU NOEUD DU MAILLAGE
  44 :            CALL INITD(UELEM,NDIM,0.D0)
  45 :            NOEL=XELE(/2)
  46 :            SEGINI HYTRAV
  46 :            CALL SHAPE(XY(1),XY(2),XY(3),ITP,SHP,IRET)
  47 : C          write(6,*)'shp ',(shp(1,j),j=1,noel)
  48 :            DO 40 J=1,NOEL
  49 :            DO 35 I=1,NDIM
  50 :            XEL(I,J)=XELE(I,J)
  51 :    35      CONTINUE
  52 :    40      CONTINUE
  52 :            CALL MATJAC(XEL,SHP,NDIM,NOEL,ZJAC,DJAC)
  53 : C          write(6,*)' zjac',((zjac(i,j),i=1,ndim),j=1,ndim)
  54 :            NTYP=NUMHYB(ITP)
  54 :            CALL HYSHP(XY(1),XY(2),XY(3),NTYP,NOEL,12,SHY,IRET)
  55 : C           write(6,*)' shy ',((shy(i,j),i=1,ndim),j=1,noel)
  56 :            I2=UN(/2)
  57 :            DO 50 IDDL=1,I2
  57 :                 CALL INITD(VECT,NDIM,0.D0)
  58 :                 DO 60 J=1,NDIM
  59 :                      DO 55 I=1,NDIM
  60 :                         VECT(I)=VECT(I)+ZJAC(I,J)*SHY(J,IDDL)
  61 :    55                CONTINUE
  62 :    60           CONTINUE
  63 :                 COEF=UN(1,IDDL,IEL)/DJAC
  64 : C       write(6,*)' flux ',iddl,iel,UN(1,IDDL,IEL),coef,djac
  65 :                      DO 65 I=1,NDIM
  66 :                           UELEM(I)=UELEM(I)+COEF*VECT(I)
  67 :    65                CONTINUE
  68 :    50           CONTINUE
  69 :            SEGSUP HYTRAV
  70 :       RETURN
  71 :       END                                                                       
  72 :