Esope source files

   1 : C CHOINT    SOURCE    CHAT      05/01/12    22:00:51     5004           
   2 :       SUBROUTINE CHOINT(NDIM,ITYG,XDEP2,XCON2,XARI2,IZSH
   3 :      $     ,XINTER,XNORMA,ITEST,JFACE)
   4 : *********************************************************************
   5 : *** SP 'CHOINT' : a partir d'une avancee de particule par convection
   6 : *** + diffusion qui traverse une ou plusieurs faces, permet de choisir
   7 : *** la face considérée pour le recalibrage de la particule.
   8 : ***
   9 : *** APPELES 1 = aucun
  10 : *** APPELES 2 = 'FACTRA', 'FACTR1', 'FACTR2', 'INTPRO'
  11 : ***
  12 : *** E = 'NDIM' dimension de l'espace
  13 : ***     'ITYG' entier caracterisant la géometrie de l'element 
  14 : ***     'XDEP2' coordonnees reelles du pt depart
  15 : ***     'XCON2' coordonnees reelles du pt apres convection
  16 : ***     'XARI2' coordonnees reelles du pt apres diffusion (arrivee)
  17 : ***     'IZSH' segmt content coord reelles noeuds elemt considéré
  18 : ***
  19 : *** S = 'XINTER' coord reelles du pt intersection associé à face 'JFACE'
  20 : ***     'XNORMA' vecteur unitaire normal à face de n° 'JFACE'
  21 : ***     'ITEST' vaut 1 si intersection considérée lors de la convection
  22 : ***             , 2 lors de la diffusion, 0 sinon (pas de traversee !)
  23 : ***
  24 : *** E/S = 'JFACE' n° face local traversee considérée pour recalibrage
  25 : ***
  26 : ***    Auteur Cyril Nou
  27 : *************************************************************************
  28 : 
  29 :       IMPLICIT INTEGER(I-N)
  30 :       IMPLICIT REAL*8 (A-H,O-Z)
  31 :       SEGMENT IZSH
  32 :          REAL*8 SHP(6,MNO9),SHY(12,MNO9),XYZL(3,MNO9)
  33 :       ENDSEGMENT
  34 :       DIMENSION XDEP2(3),XCON2(3),XARI2(3),XINTER(3),XNORMA(3)
  35 :       DIMENSION IFACE(6),XINT(3,6),XN(3,6)
  36 : ***   on postule pas de traversee au depart
  37 :       ITEST=0
  38 : 
  38 :       CALL FACTRA(NDIM,ITYG,XCON2,XARI2,IZSH,IFACE,XINT,XN,NBFAC)
  39 :       IF (NBFAC.NE.0) THEN
  40 : ***      cas ou face(s) traversee(s) lors du saut diffusif
  40 :          CALL INTPRO(NDIM,XARI2,NBFAC,IFACE,XINT,XN,JFACE,XINTER,XNORMA)
  41 :          ITEST=2
  42 :       ELSE
  42 :          CALL FACTRA(NDIM,ITYG,XDEP2,XCON2,IZSH,IFACE,XINT,XN,NBFAC)
  43 :          IF (NBFAC.NE.0) THEN
  44 : ***         cas ou face(s) traversee(s) lors du saut convectif
  44 :             CALL INTPRO(NDIM,XCON2,NBFAC
  45 :      $           ,IFACE,XINT,XN,JFACE,XINTER,XNORMA)
  46 :             ITEST=1
  47 :          ENDIF
  48 :       ENDIF
  49 : 
  50 :       RETURN
  51 :       END
  52 :  
  53 :