--- trunk/Sources/phylmd/hbtm.f 2016/08/30 12:52:46 206 +++ trunk/Sources/phylmd/hbtm.f 2018/01/22 15:02:56 252 @@ -5,7 +5,7 @@ contains SUBROUTINE HBTM(paprs, pplay, t2m, q2m, ustar, flux_t, flux_q, u, v, t, q, & - pblh, cape, EauLiq, ctei, pblT, therm, trmb1, trmb2, trmb3, plcl) + pblh, cape, EauLiq, ctei, pblT, therm, plcl) ! D'apr\'es Holstag et Boville et Troen et Mahrt ! JAS 47 BLM @@ -42,13 +42,12 @@ REAL, intent(in):: t2m(klon) ! temperature a 2 m ! q a 2 et 10m REAL, intent(in):: q2m(klon) - REAL, intent(in):: ustar(klon) + REAL, intent(in):: ustar(:) ! (knon) REAL, intent(in):: flux_t(:), flux_q(:) ! (knon) flux à la surface - ! vitesse U (m/s) - REAL, intent(in):: u(klon, klev) - ! vitesse V (m/s) - REAL, intent(in):: v(klon, klev) + REAL, intent(in):: u(klon, klev) ! vitesse U (m/s) + REAL, intent(in):: v(klon, klev) ! vitesse V (m/s) + ! temperature (K) REAL, intent(in):: t(klon, klev) ! vapeur d'eau (kg/kg) @@ -64,7 +63,6 @@ REAL pblT(klon) ! thermal virtual temperature excess REAL therm(klon) - REAL trmb1(klon), trmb2(klon), trmb3(klon) REAL plcl(klon) ! Local: @@ -207,9 +205,6 @@ plcl(i) = 6000. ! Lambda = -u*^3 / (alpha.g.kvon. obklen(i) = -t(i, 1)*ustar(i)**3/(RG*vk*heatv(i)) - trmb1(i) = 0. - trmb2(i) = 0. - trmb3(i) = 0. ENDDO ! PBL height calculation: Search for level of pbl. Scan upward