Télécharger la plaquette Cast3M


fr en

Télécharger kfpa.eso

Retour à la liste

Numérotation des lignes : 

kfpa

  1. C KFPA      SOURCE    CB215821  20/11/25    13:31:22     10792          
  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.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. -INC SMCHPOI
  28.       POINTEUR  MUET.MCHPOI, MNORM.MCHPOI, MAK.MCHPOI
  29.       POINTEUR  NUET.MPOVAL, NNORM.MPOVAL, NAK.MPOVAL
  30. -INC SMCOORD
  31. -INC SMELEME
  32.       POINTEUR  IGEO.MELEME,IGEO0.MELEME
  33. C
  34.       DIMENSION ROG(3)
  35. C
  36.       CALL LIRREE(ANU,1,IRET)
  37.       IF (IRET.EQ.0) RETURN
  38.       CALL LIRREE(YP,1,IRET)
  39.       IF (IRET.EQ.0) RETURN
  40.       CALL LIROBJ('CHPOINT',MUET,1,IRET)
  41.       IF (IRET.EQ.0) RETURN
  42.       CALL LIROBJ('CHPOINT',MNORM,1,IRET)
  43.       IF (IRET.EQ.0) RETURN
  44.       CALL LIROBJ('POINT',IP,1,IRET)
  45.       IF (IRET.EQ.0) RETURN
  46.       CALL LIRREE(RAP,1,IRET)
  47.       IF (IRET.EQ.0) RETURN
  48. C
  49. C
  50.       CALL LICHT(MUET,NUET,TYPE,IGEO)
  51.       CALL LICHT(MNORM,NNORM,TYPE,IGEO0)
  52.       IF (IGEO.NE.IGEO0)THEN
  53.         WRITE(6,100)
  54.         RETURN
  55.       ENDIF
  56. C
  57.       TYPE='CENTRE'
  58.       NC = 1
  59.       CALL CRCHPT(TYPE,IGEO,NC,MAK)
  60.       CALL LICHT(MAK,NAK,TYPE,IGEO)
  61.       N = NAK.VPOCHA(/1)
  62. C
  63.       segact mcoord
  64.       DO 10 I=1,IDIM
  65.             ROG(I) = XCOOR((IP-1)*(IDIM+1)+I)
  66.  10   CONTINUE
  67.       CALL XFPA(ANU,YP,ROG,RAP,N,
  68.      1          NUET.VPOCHA,NNORM.VPOCHA,NAK.VPOCHA)
  69. C
  70.       CALL ECROBJ('CHPOINT',MAK)
  71. C
  72.       SEGDES MAK,NAK,MUET,NUET,MNORM,NNORM
  73. C
  74.       RETURN
  75.  100  FORMAT(2X,'Les chpoints n''ont pas le meme SPG')
  76.       END
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  

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