11 |
USE dimens_m |
USE dimens_m |
12 |
USE dimphy |
USE dimphy |
13 |
USE tracstoke |
USE tracstoke |
14 |
USE yomcst |
USE suphec_m |
15 |
USE yoethf |
USE yoethf_m |
16 |
USE fcttre |
USE fcttre |
17 |
USE comfisrtilp |
USE comfisrtilp |
18 |
use numer_rec, only: erf |
use numer_rec, only: nr_erf |
19 |
|
|
20 |
IMPLICIT NONE |
IMPLICIT NONE |
21 |
|
|
24 |
REAL, INTENT (IN) :: & ! intervalle du temps (s) |
REAL, INTENT (IN) :: & ! intervalle du temps (s) |
25 |
dtime |
dtime |
26 |
REAL, INTENT (IN) :: paprs(klon,klev+1) ! pression a inter-couche |
REAL, INTENT (IN) :: paprs(klon,klev+1) ! pression a inter-couche |
27 |
REAL, INTENT (IN) :: pplay(klon,klev) ! pression au milieu de couche |
REAL, INTENT (IN):: pplay(klon,klev) ! pression au milieu de couche |
28 |
REAL t(klon,klev) ! temperature (K) |
REAL, INTENT (IN):: t(klon,klev) ! temperature (K) |
29 |
REAL q(klon,klev) ! humidite specifique (kg/kg) |
REAL q(klon,klev) ! humidite specifique (kg/kg) |
30 |
REAL d_t(klon,klev) ! incrementation de la temperature (K) |
REAL d_t(klon,klev) ! incrementation de la temperature (K) |
31 |
REAL d_q(klon,klev) ! incrementation de la vapeur d'eau |
REAL d_q(klon,klev) ! incrementation de la vapeur d'eau |
338 |
zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.)) |
zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.)) |
339 |
zpdf_e1(i) = zpdf_a(i) - zpdf_b(i) |
zpdf_e1(i) = zpdf_a(i) - zpdf_b(i) |
340 |
zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i)) |
zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i)) |
341 |
zpdf_e1(i) = 1. - erf(zpdf_e1(i)) |
zpdf_e1(i) = 1. - nr_erf(zpdf_e1(i)) |
342 |
zpdf_e2(i) = zpdf_a(i) + zpdf_b(i) |
zpdf_e2(i) = zpdf_a(i) + zpdf_b(i) |
343 |
zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i)) |
zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i)) |
344 |
zpdf_e2(i) = 1. - erf(zpdf_e2(i)) |
zpdf_e2(i) = 1. - nr_erf(zpdf_e2(i)) |
345 |
IF (zpdf_e1(i)<1.E-10) THEN |
IF (zpdf_e1(i)<1.E-10) THEN |
346 |
rneb(i,k) = 0. |
rneb(i,k) = 0. |
347 |
zqn(i) = zqs(i) |
zqn(i) = zqs(i) |