Télécharger choint.eso

Retour à la liste

Numérotation des lignes :

  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