18 |
use climb_hq_up_m, only: climb_hq_up |
use climb_hq_up_m, only: climb_hq_up |
19 |
USE dimphy, ONLY: klev, klon |
USE dimphy, ONLY: klev, klon |
20 |
USE interfsurf_hq_m, ONLY: interfsurf_hq |
USE interfsurf_hq_m, ONLY: interfsurf_hq |
21 |
|
USE suphec_m, ONLY: rkappa |
22 |
|
|
23 |
REAL, intent(in):: dtime ! intervalle du temps (s) |
REAL, intent(in):: dtime ! intervalle du temps (s) |
24 |
integer, intent(in):: julien ! jour de l'annee en cours |
integer, intent(in):: julien ! jour de l'annee en cours |
100 |
|
|
101 |
! Local: |
! Local: |
102 |
|
|
103 |
|
INTEGER k |
104 |
REAL evap(size(knindex)) ! (knon) evaporation au sol |
REAL evap(size(knindex)) ! (knon) evaporation au sol |
105 |
REAL, dimension(size(knindex), klev):: cq, dq, ch, dh ! (knon, klev) |
REAL, dimension(size(knindex), klev):: cq, dq, ch, dh ! (knon, klev) |
106 |
REAL pkf(size(knindex), klev) ! (knon, klev) |
REAL pkf(size(knindex), klev) ! (knon, klev) |
108 |
|
|
109 |
!---------------------------------------------------------------- |
!---------------------------------------------------------------- |
110 |
|
|
111 |
|
forall (k = 1:klev) pkf(:, k) = (paprs(:, 1) / pplay(:, k))**RKAPPA |
112 |
|
! (La pression de r\'ef\'erence est celle au sol.) |
113 |
|
|
114 |
call climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, dtime, & |
call climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, dtime, & |
115 |
delp, q) |
delp, q) |
116 |
CALL interfsurf_hq(dtime, julien, rmu0, nisurf, knindex, debut, tsoil, & |
CALL interfsurf_hq(dtime, julien, rmu0, nisurf, knindex, debut, tsoil, & |
117 |
qsol, u1lay, v1lay, t(:, 1), q(:, 1), tq_cdrag, ch(:, 1), cq(:, 1), & |
qsol, u1lay, v1lay, t(:, 1), q(:, 1), tq_cdrag, ch(:, 1), cq(:, 1), & |
118 |
dh(:, 1), dq(:, 1), precip_rain, precip_snow, rugos, rugoro, snow, & |
dh(:, 1), dq(:, 1), precip_rain, precip_snow, rugos, rugoro, snow, & |