C KFPA      SOURCE    CHAT      05/01/13    00:55:58     5004           
      SUBROUTINE KFPA
C*****************************************************************
C
C OBJET : calcul de la vitesse de depot d'aerosols AK dans chaque
C         element de la ligne paroi .
C
C SYNTAXE : AK = KFPA NU YP UET NORM ROG RAP 
C
C AK   : CHPOINT SCAL CENTRE (vitesse de depot)
C NU   : FLOTTANT (viscosite)
C YP   : FLOTTANT (epaisseur de la couche limite)
C UET  : CHPOINT SCAL CENTRE (vitesse de frottement)
C NORM : CHPOINT VECT FACE (normale a la paroi)
C ROG  : POINT (masse volumique * g)
C RAP  : FLOTTANT (rayon des particules)
C
C*****************************************************************
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      CHARACTER*8 TYPE
C
-INC CCOPTIO
-INC SMCHPOI
      POINTEUR  MUET.MCHPOI, MNORM.MCHPOI, MAK.MCHPOI
      POINTEUR  NUET.MPOVAL, NNORM.MPOVAL, NAK.MPOVAL
-INC SMCOORD
-INC SMELEME
      POINTEUR  IGEO.MELEME,IGEO0.MELEME
C
      DIMENSION ROG(3)
C
      CALL LIRREE(ANU,1,IRET)
      IF (IRET.EQ.0) RETURN 
      CALL LIRREE(YP,1,IRET)
      IF (IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT',MUET,1,IRET)
      IF (IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT',MNORM,1,IRET)
      IF (IRET.EQ.0) RETURN 
      CALL LIROBJ('POINT',IP,1,IRET)
      IF (IRET.EQ.0) RETURN
      CALL LIRREE(RAP,1,IRET)
      IF (IRET.EQ.0) RETURN
C
C
      CALL LICHT(MUET,NUET,TYPE,IGEO)
      CALL LICHT(MNORM,NNORM,TYPE,IGEO0) 
      IF (IGEO.NE.IGEO0)THEN
        WRITE(6,100)
        RETURN
      ENDIF
C
      TYPE='CENTRE'
      NC = 1
      CALL CRCHPT(TYPE,IGEO,NC,MAK)
      CALL LICHT(MAK,NAK,TYPE,IGEO)
      N = NAK.VPOCHA(/1)
C
      DO 10 I=1,IDIM
            ROG(I) = XCOOR((IP-1)*(IDIM+1)+I) 
 10   CONTINUE
      CALL XFPA(ANU,YP,ROG,RAP,N,
     1          NUET.VPOCHA,NNORM.VPOCHA,NAK.VPOCHA)
C
      CALL ECROBJ('CHPOINT',MAK)
C
      SEGDES MAK,NAK,MUET,NUET,MNORM,NNORM
C
      RETURN
 100  FORMAT(2X,'Les chpoints n''ont pas le meme SPG')
      END