Télécharger kfpa.eso

Retour à la liste

Numérotation des lignes :

  1. C KFPA SOURCE CHAT 05/01/13 00:55:58 5004
  2. SUBROUTINE KFPA
  3. C*****************************************************************
  4. C
  5. C OBJET : calcul de la vitesse de depot d'aerosols AK dans chaque
  6. C element de la ligne paroi .
  7. C
  8. C SYNTAXE : AK = KFPA NU YP UET NORM ROG RAP
  9. C
  10. C AK : CHPOINT SCAL CENTRE (vitesse de depot)
  11. C NU : FLOTTANT (viscosite)
  12. C YP : FLOTTANT (epaisseur de la couche limite)
  13. C UET : CHPOINT SCAL CENTRE (vitesse de frottement)
  14. C NORM : CHPOINT VECT FACE (normale a la paroi)
  15. C ROG : POINT (masse volumique * g)
  16. C RAP : FLOTTANT (rayon des particules)
  17. C
  18. C*****************************************************************
  19. C
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8 (A-H,O-Z)
  22. CHARACTER*8 TYPE
  23. C
  24. -INC CCOPTIO
  25. -INC SMCHPOI
  26. POINTEUR MUET.MCHPOI, MNORM.MCHPOI, MAK.MCHPOI
  27. POINTEUR NUET.MPOVAL, NNORM.MPOVAL, NAK.MPOVAL
  28. -INC SMCOORD
  29. -INC SMELEME
  30. POINTEUR IGEO.MELEME,IGEO0.MELEME
  31. C
  32. DIMENSION ROG(3)
  33. C
  34. CALL LIRREE(ANU,1,IRET)
  35. IF (IRET.EQ.0) RETURN
  36. CALL LIRREE(YP,1,IRET)
  37. IF (IRET.EQ.0) RETURN
  38. CALL LIROBJ('CHPOINT',MUET,1,IRET)
  39. IF (IRET.EQ.0) RETURN
  40. CALL LIROBJ('CHPOINT',MNORM,1,IRET)
  41. IF (IRET.EQ.0) RETURN
  42. CALL LIROBJ('POINT',IP,1,IRET)
  43. IF (IRET.EQ.0) RETURN
  44. CALL LIRREE(RAP,1,IRET)
  45. IF (IRET.EQ.0) RETURN
  46. C
  47. C
  48. CALL LICHT(MUET,NUET,TYPE,IGEO)
  49. CALL LICHT(MNORM,NNORM,TYPE,IGEO0)
  50. IF (IGEO.NE.IGEO0)THEN
  51. WRITE(6,100)
  52. RETURN
  53. ENDIF
  54. C
  55. TYPE='CENTRE'
  56. NC = 1
  57. CALL CRCHPT(TYPE,IGEO,NC,MAK)
  58. CALL LICHT(MAK,NAK,TYPE,IGEO)
  59. N = NAK.VPOCHA(/1)
  60. C
  61. DO 10 I=1,IDIM
  62. ROG(I) = XCOOR((IP-1)*(IDIM+1)+I)
  63. 10 CONTINUE
  64. CALL XFPA(ANU,YP,ROG,RAP,N,
  65. 1 NUET.VPOCHA,NNORM.VPOCHA,NAK.VPOCHA)
  66. C
  67. CALL ECROBJ('CHPOINT',MAK)
  68. C
  69. SEGDES MAK,NAK,MUET,NUET,MNORM,NNORM
  70. C
  71. RETURN
  72. 100 FORMAT(2X,'Les chpoints n''ont pas le meme SPG')
  73. END
  74.  
  75.  
  76.  
  77.  
  78.  

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