4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
subroutine climb_hq_up(d_t, d_q, cq, dq, ch, dh, flux_t, flux_q, dtime, pkf, & |
subroutine climb_hq_up(d_t, d_q, cq, dq, ch, dh, flux_t, flux_q, pkf, t, q) |
|
t, q) |
|
8 |
|
|
9 |
|
use comconst, only: dtphys |
10 |
USE dimphy, ONLY: klev |
USE dimphy, ONLY: klev |
11 |
USE suphec_m, ONLY: rcpd |
USE suphec_m, ONLY: rcpd |
12 |
|
|
21 |
REAL, intent(out):: flux_q(:) ! (knon) |
REAL, intent(out):: flux_q(:) ! (knon) |
22 |
! flux de la vapeur d'eau à la surface, en kg / (m**2 s) |
! flux de la vapeur d'eau à la surface, en kg / (m**2 s) |
23 |
|
|
|
REAL, intent(in):: dtime ! intervalle du temps (s) |
|
24 |
REAL, intent(in):: pkf(:, :) ! (knon, klev) |
REAL, intent(in):: pkf(:, :) ! (knon, klev) |
25 |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
26 |
REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg) |
REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg) |
32 |
|
|
33 |
!---------------------------------------------------------------------- |
!---------------------------------------------------------------------- |
34 |
|
|
35 |
h(:, 1) = ch(:, 1) + dh(:, 1) * flux_t * dtime |
h(:, 1) = ch(:, 1) + dh(:, 1) * flux_t * dtphys |
36 |
local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtime |
local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtphys |
37 |
|
|
38 |
DO k = 2, klev |
DO k = 2, klev |
39 |
h(:, k) = ch(:, k) + dh(:, k) * h(:, k - 1) |
h(:, k) = ch(:, k) + dh(:, k) * h(:, k - 1) |