Télécharger la plaquette Cast3M


fr en

Télécharger choint.eso

Retour à la liste

Numérotation des lignes : 

choint

  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.  
  39.       CALL FACTRA(NDIM,ITYG,XCON2,XARI2,IZSH,IFACE,XINT,XN,NBFAC)
  40.       IF (NBFAC.NE.0) THEN
  41. ***      cas ou face(s) traversee(s) lors du saut diffusif
  42.          CALL INTPRO(NDIM,XARI2,NBFAC,IFACE,XINT,XN,JFACE,XINTER,XNORMA)
  43.          ITEST=2
  44.       ELSE
  45.          CALL FACTRA(NDIM,ITYG,XDEP2,XCON2,IZSH,IFACE,XINT,XN,NBFAC)
  46.          IF (NBFAC.NE.0) THEN
  47. ***         cas ou face(s) traversee(s) lors du saut convectif
  48.             CALL INTPRO(NDIM,XCON2,NBFAC
  49.      $           ,IFACE,XINT,XN,JFACE,XINTER,XNORMA)
  50.             ITEST=1
  51.          ENDIF
  52.       ENDIF
  53.  
  54.       RETURN
  55.       END
  56.  
  57.  
  58.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales