Esope source files

   1 : C SAUCON    SOURCE    CHAT      05/01/13    03:11:11     5004           
   2 :       SUBROUTINE SAUCON(NDIM,ITY1,ITYG,NOEL1,NOUN
   3 :      $     ,DIAM,UELEM,XDEP,XDEP2,DTREEL,XARI2,IZSH,UREEL,ITEST)
   4 : 
   5 : *************************************************************************
   6 : *** SP 'SAUCON' : fait avancer la particule apres un pas de tps par
   7 : *** convection explicite uniquement dans l'element reel considere
   8 : ***
   9 : *** APPELES 1 = aucun
  10 : *** APPELES 2 = 'VITAPP'
  11 : ***
  12 : *** E=  'NDIM' dimension de l'espace
  13 : ***     'ITY1' entier caracterisant le type de l'element
  14 : ***     'ITYG' entier caracterisant la geometrie de l'element
  15 : ***     'NOEL1' nombre de noeuds de l'element considere
  16 : ***     'NOUN' nombre de flux consideres
  17 : ***     'DIAM' "longueur caracteristique" de l'element
  18 : ***     'UELEM' valeurs des flux aux faces
  19 : ***     'XDEP' coordonnees reference du pt de depart
  20 : ***     'XDEP2' coordonnees reelles du pt de depart
  21 : ***     'DTREEL' pas de tps reel pour avancee de la particule
  22 : ***
  23 : *** S = 'XARI2' coordonnees reelles d'arrivee
  24 : ***     'ITEST' vaut 0 si jacobien nul (ds approximation vitesse), 1 sinon
  25 : ***     'UREEL' vitesse convective reelle
  26 : ***
  27 : *** E/S = 'IZSH' segment content pour l'elemt considere et au pt depart :
  28 : ***             'XYZL' coordonnees reelles des noeuds (E)
  29 : ***             'SHP' fonctions de forme et derivees ds elemt ref (S)
  30 : ***             'SHY' fonctions de base et derivees ds elemt ref (S)
  31 : ***
  32 : ***          auteur Cyril Nou
  33 : *************************************************************************
  34 : 
  35 :       IMPLICIT INTEGER(I-N)
  36 :       IMPLICIT REAL*8 (A-H,O-Z)
  37 : *-INC SMELEME
  38 : *      POINTEUR IZCENT.MELEME,IELTFA.MELEME
  39 :       SEGMENT IZSH
  40 :           REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  41 :       ENDSEGMENT
  42 :       DIMENSION XDEP(3),XDEP2(3),UELEM(NOUN)
  43 :       DIMENSION XARI2(3),UREF(3),UREEL(3)
  44 : 
  45 : ***   calcul de la vitesse reelle au pt de depart
  45 :       CALL VITAPP(NDIM,ITY1,ITYG,NOEL1
  46 :      $     ,NOUN,DIAM,UELEM,XDEP,UREF,UREEL,IZSH,ITEST)
  47 : ***   calcul pt reel d'arrivee apres pas de tps 'DTREEL'
  48 :       DO 10 I=1,NDIM
  49 :          XARI2(I)=XDEP2(I)+UREEL(I)*DTREEL
  50 :  10   CONTINUE
  51 :       RETURN                                                                    
  52 :       END                                                                       
  53 :  
  54 :  
  55 : 
  56 : 
  57 :  
  58 :