8 |
flux_q, u, v, t, q, pblh, cape, EauLiq, ctei, pblT, therm, trmb1, & |
flux_q, u, v, t, q, pblh, cape, EauLiq, ctei, pblT, therm, trmb1, & |
9 |
trmb2, trmb3, plcl) |
trmb2, trmb3, plcl) |
10 |
|
|
11 |
use dimens_m |
! D'après Holstag et Boville et Troen et Mahrt |
|
use dimphy |
|
|
use SUPHEC_M |
|
|
use yoethf_m |
|
|
use fcttre |
|
|
|
|
|
! D'apres Holstag & Boville et Troen & Mahrt |
|
12 |
! JAS 47 BLM |
! JAS 47 BLM |
13 |
! Algorithme thèse Anne Mathieu |
! Algorithme thèse Anne Mathieu |
14 |
! Critère d'entraînement Peter Duynkerke (JAS 50) |
! Critère d'entraînement Peter Duynkerke (JAS 50) |
29 |
! mais on garde la possibilité de changer si besoin est (jusqu'à présent |
! mais on garde la possibilité de changer si besoin est (jusqu'à présent |
30 |
! la forme de HB avec le 1er niveau modele etait conservee) |
! la forme de HB avec le 1er niveau modele etait conservee) |
31 |
|
|
32 |
|
USE dimphy, ONLY: klev, klon, max |
33 |
|
USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlvtt, rtt, rv |
34 |
|
USE yoethf_m, ONLY: r2es, rvtmp2 |
35 |
|
USE fcttre, ONLY: foeew |
36 |
|
|
37 |
REAL RLvCp, REPS |
REAL RLvCp, REPS |
38 |
! Arguments: |
! Arguments: |
39 |
|
|
446 |
|
|
447 |
! Main level loop to compute the diffusivities and |
! Main level loop to compute the diffusivities and |
448 |
! counter-gradient terms: |
! counter-gradient terms: |
449 |
DO k = 2, isommet |
loop_level: DO k = 2, isommet |
450 |
! Find levels within boundary layer: |
! Find levels within boundary layer: |
451 |
DO i = 1, knon |
DO i = 1, knon |
452 |
unslev(i) = .FALSE. |
unslev(i) = .FALSE. |
516 |
|
|
517 |
! For all layers, compute integral info and CTEI |
! For all layers, compute integral info and CTEI |
518 |
DO i = 1, knon |
DO i = 1, knon |
519 |
if (check(i).or.omegafl(i)) then |
if (check(i) .or. omegafl(i)) then |
520 |
if (.not.Zsat(i)) then |
if (.not. Zsat(i)) then |
521 |
T2 = T2m(i) * s(i, k) |
T2 = T2m(i) * s(i, k) |
522 |
! thermodyn functions |
! thermodyn functions |
523 |
zdelta=MAX(0., SIGN(1., RTT - T2)) |
zdelta=MAX(0., SIGN(1., RTT - T2)) |
542 |
! cette ligne a deja ete faite normalement ? |
! cette ligne a deja ete faite normalement ? |
543 |
endif |
endif |
544 |
ENDDO |
ENDDO |
545 |
end DO |
end DO loop_level |
546 |
|
|
547 |
END SUBROUTINE HBTM |
END SUBROUTINE HBTM |
548 |
|
|