5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE HBTM(paprs, pplay, t2m, q2m, ustar, flux_t, flux_q, u, v, t, q, & |
SUBROUTINE HBTM(paprs, pplay, t2m, q2m, ustar, flux_t, flux_q, u, v, t, q, & |
8 |
pblh, cape, EauLiq, ctei, pblT, therm, trmb1, trmb2, trmb3, plcl) |
pblh, cape, EauLiq, ctei, pblT, therm, plcl) |
9 |
|
|
10 |
! D'apr\'es Holstag et Boville et Troen et Mahrt |
! D'apr\'es Holstag et Boville et Troen et Mahrt |
11 |
! JAS 47 BLM |
! JAS 47 BLM |
44 |
REAL, intent(in):: q2m(klon) |
REAL, intent(in):: q2m(klon) |
45 |
REAL, intent(in):: ustar(:) ! (knon) |
REAL, intent(in):: ustar(:) ! (knon) |
46 |
REAL, intent(in):: flux_t(:), flux_q(:) ! (knon) flux à la surface |
REAL, intent(in):: flux_t(:), flux_q(:) ! (knon) flux à la surface |
47 |
|
REAL, intent(in):: u(:, :) ! (knon, klev) vitesse U (m/s) |
48 |
REAL, intent(in):: u(klon, klev) ! vitesse U (m/s) |
REAL, intent(in):: v(:, :) ! (knon, klev) vitesse V (m/s) |
49 |
REAL, intent(in):: v(klon, klev) ! vitesse V (m/s) |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
50 |
|
REAL, intent(in):: q(:, :) ! (knon, klev) vapeur d'eau (kg/kg) |
|
! temperature (K) |
|
|
REAL, intent(in):: t(klon, klev) |
|
|
! vapeur d'eau (kg/kg) |
|
|
REAL, intent(in):: q(klon, klev) |
|
51 |
|
|
52 |
REAL, intent(out):: pblh(:) ! (knon) |
REAL, intent(out):: pblh(:) ! (knon) |
53 |
! Cape du thermique |
! Cape du thermique |
59 |
REAL pblT(klon) |
REAL pblT(klon) |
60 |
! thermal virtual temperature excess |
! thermal virtual temperature excess |
61 |
REAL therm(klon) |
REAL therm(klon) |
|
REAL trmb1(klon), trmb2(klon), trmb3(klon) |
|
62 |
REAL plcl(klon) |
REAL plcl(klon) |
63 |
|
|
64 |
! Local: |
! Local: |
201 |
plcl(i) = 6000. |
plcl(i) = 6000. |
202 |
! Lambda = -u*^3 / (alpha.g.kvon.<w'Theta'v> |
! Lambda = -u*^3 / (alpha.g.kvon.<w'Theta'v> |
203 |
obklen(i) = -t(i, 1)*ustar(i)**3/(RG*vk*heatv(i)) |
obklen(i) = -t(i, 1)*ustar(i)**3/(RG*vk*heatv(i)) |
|
trmb1(i) = 0. |
|
|
trmb2(i) = 0. |
|
|
trmb3(i) = 0. |
|
204 |
ENDDO |
ENDDO |
205 |
|
|
206 |
! PBL height calculation: Search for level of pbl. Scan upward |
! PBL height calculation: Search for level of pbl. Scan upward |