11 |
! From phylmd/fisrtilp.F, version 1.2 2004/11/09 16:55:40 |
! From phylmd/fisrtilp.F, version 1.2 2004/11/09 16:55:40 |
12 |
! Author: Z. X. Li (LMD/CNRS), 20 mars 1995 |
! Author: Z. X. Li (LMD/CNRS), 20 mars 1995 |
13 |
|
|
14 |
! Objet: condensation et précipitation stratiforme, schéma de |
! Objet : condensation et précipitation stratiforme, schéma de |
15 |
! nuage, schéma de condensation à grande échelle (pluie). |
! nuage, schéma de condensation à grande échelle (pluie). |
16 |
|
|
17 |
USE dimphy, ONLY: klev, klon |
USE dimphy, ONLY: klev, klon |
28 |
REAL, INTENT (IN):: paprs(klon, klev+1) ! pression a inter-couche |
REAL, INTENT (IN):: paprs(klon, klev+1) ! pression a inter-couche |
29 |
REAL, INTENT (IN):: pplay(klon, klev) ! pression au milieu de couche |
REAL, INTENT (IN):: pplay(klon, klev) ! pression au milieu de couche |
30 |
REAL, INTENT (IN):: t(klon, klev) ! temperature (K) |
REAL, INTENT (IN):: t(klon, klev) ! temperature (K) |
31 |
REAL q(klon, klev) ! humidite specifique (kg/kg) |
REAL, INTENT (IN):: q(klon, klev) ! humidite specifique (kg/kg) |
32 |
|
LOGICAL ptconv(klon, klev) ! determine la largeur de distribution de vapeur |
33 |
|
REAL ratqs(klon, klev) ! determine la largeur de distribution de vapeur |
34 |
REAL d_t(klon, klev) ! incrementation de la temperature (K) |
REAL d_t(klon, klev) ! incrementation de la temperature (K) |
35 |
REAL d_q(klon, klev) ! incrementation de la vapeur d'eau |
REAL d_q(klon, klev) ! incrementation de la vapeur d'eau |
36 |
REAL d_ql(klon, klev) ! incrementation de l'eau liquide |
REAL d_ql(klon, klev) ! incrementation de l'eau liquide |
37 |
REAL rneb(klon, klev) ! fraction nuageuse |
REAL rneb(klon, klev) ! fraction nuageuse |
38 |
REAL radliq(klon, klev) ! eau liquide utilisee dans rayonnements |
REAL radliq(klon, klev) ! eau liquide utilisee dans rayonnements |
|
REAL rhcl(klon, klev) ! humidite relative en ciel clair |
|
39 |
REAL rain(klon) ! pluies (mm/s) |
REAL rain(klon) ! pluies (mm/s) |
40 |
REAL snow(klon) ! neige (mm/s) |
REAL snow(klon) ! neige (mm/s) |
|
REAL prfl(klon, klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s) |
|
|
REAL psfl(klon, klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s) |
|
|
! Coeffients de fraction lessivee : pour OFF-LINE |
|
41 |
|
|
42 |
|
! Coeffients de fraction lessivee : pour OFF-LINE |
43 |
|
REAL pfrac_impa(klon, klev) |
44 |
REAL pfrac_nucl(klon, klev) |
REAL pfrac_nucl(klon, klev) |
45 |
REAL pfrac_1nucl(klon, klev) |
REAL pfrac_1nucl(klon, klev) |
|
REAL pfrac_impa(klon, klev) |
|
46 |
|
|
47 |
! Fraction d'aerosols lessivee par impaction et par nucleation |
! Fraction d'aerosols lessivee par impaction et par nucleation |
48 |
! POur ON-LINE |
! POur ON-LINE |
49 |
|
REAL frac_nucl(klon, klev) |
50 |
|
|
51 |
|
REAL prfl(klon, klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s) |
52 |
|
REAL psfl(klon, klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s) |
53 |
|
REAL rhcl(klon, klev) ! humidite relative en ciel clair |
54 |
|
|
55 |
|
! Local: |
56 |
|
|
57 |
|
! Fraction d'aerosols lessivee par impaction et par nucleation |
58 |
|
! POur ON-LINE |
59 |
REAL frac_impa(klon, klev) |
REAL frac_impa(klon, klev) |
|
REAL frac_nucl(klon, klev) |
|
60 |
REAL zct(klon), zcl(klon) |
REAL zct(klon), zcl(klon) |
61 |
!AA |
!AA |
62 |
|
|
69 |
PARAMETER (ninter=5) |
PARAMETER (ninter=5) |
70 |
LOGICAL evap_prec ! evaporation de la pluie |
LOGICAL evap_prec ! evaporation de la pluie |
71 |
PARAMETER (evap_prec=.TRUE.) |
PARAMETER (evap_prec=.TRUE.) |
|
REAL ratqs(klon, klev) ! determine la largeur de distribution de vapeur |
|
|
LOGICAL ptconv(klon, klev) ! determine la largeur de distribution de vapeur |
|
72 |
REAL zpdf_sig(klon), zpdf_k(klon), zpdf_delta(klon) |
REAL zpdf_sig(klon), zpdf_k(klon), zpdf_delta(klon) |
73 |
REAL zpdf_a(klon), zpdf_b(klon), zpdf_e1(klon), zpdf_e2(klon) |
REAL zpdf_a(klon), zpdf_b(klon), zpdf_e1(klon), zpdf_e2(klon) |
74 |
|
|