Esope source files

   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
  33 :       CALL LIRREE(ANU,1,IRET)
  34 :       IF (IRET.EQ.0) RETURN 
  34 :       CALL LIRREE(YP,1,IRET)
  35 :       IF (IRET.EQ.0) RETURN
  35 :       CALL LIROBJ('CHPOINT',MUET,1,IRET)
  36 :       IF (IRET.EQ.0) RETURN
  36 :       CALL LIROBJ('CHPOINT',MNORM,1,IRET)
  37 :       IF (IRET.EQ.0) RETURN 
  37 :       CALL LIROBJ('POINT',IP,1,IRET)
  38 :       IF (IRET.EQ.0) RETURN
  38 :       CALL LIRREE(RAP,1,IRET)
  39 :       IF (IRET.EQ.0) RETURN
  40 : C
  41 : C
  41 :       CALL LICHT(MUET,NUET,TYPE,IGEO)
  41 :       CALL LICHT(MNORM,NNORM,TYPE,IGEO0) 
  42 :       IF (IGEO.NE.IGEO0)THEN
  43 :         WRITE(6,100)
  44 :         RETURN
  45 :       ENDIF
  46 : C
  47 :       TYPE='CENTRE'
  48 :       NC = 1
  48 :       CALL CRCHPT(TYPE,IGEO,NC,MAK)
  48 :       CALL LICHT(MAK,NAK,TYPE,IGEO)
  49 :       N = NAK.VPOCHA(/1)
  50 : C
  51 :       DO 10 I=1,IDIM
  52 :             ROG(I) = XCOOR((IP-1)*(IDIM+1)+I) 
  53 :  10   CONTINUE
  53 :       CALL XFPA(ANU,YP,ROG,RAP,N,
  54 :      1          NUET.VPOCHA,NNORM.VPOCHA,NAK.VPOCHA)
  55 : C
  55 :       CALL ECROBJ('CHPOINT',MAK)
  56 : C
  57 :       SEGDES MAK,NAK,MUET,NUET,MNORM,NNORM
  58 : C
  59 :       RETURN
  60 :  100  FORMAT(2X,'Les chpoints n''ont pas le meme SPG')
  61 :       END
  62 :  
  63 :  
  64 :  
  65 :