1 |
guez |
296 |
module climb_hq_up_m |
2 |
|
|
|
3 |
|
|
IMPLICIT none |
4 |
|
|
|
5 |
|
|
contains |
6 |
|
|
|
7 |
guez |
299 |
subroutine climb_hq_up(d_t, d_q, cq, dq, ch, dh, flux_t, flux_q, pkf, t, q) |
8 |
guez |
296 |
|
9 |
guez |
299 |
use comconst, only: dtphys |
10 |
guez |
296 |
USE dimphy, ONLY: klev |
11 |
|
|
USE suphec_m, ONLY: rcpd |
12 |
|
|
|
13 |
|
|
REAL, intent(out):: d_t(:, :) ! (knon, klev) incrementation de "t" |
14 |
|
|
REAL, intent(out):: d_q(:, :) ! (knon, klev) incrementation de "q" |
15 |
|
|
REAL, intent(in), dimension(:, :):: cq, dq, ch, dh ! (knon, klev) |
16 |
|
|
|
17 |
|
|
REAL, intent(in):: flux_t(:) ! (knon) |
18 |
|
|
! (diagnostic) flux de chaleur sensible (Cp T) à la surface, |
19 |
|
|
! positif vers le bas, W / m2 |
20 |
|
|
|
21 |
|
|
REAL, intent(out):: flux_q(:) ! (knon) |
22 |
|
|
! flux de la vapeur d'eau à la surface, en kg / (m**2 s) |
23 |
|
|
|
24 |
|
|
REAL, intent(in):: pkf(:, :) ! (knon, klev) |
25 |
|
|
REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K) |
26 |
|
|
REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg) |
27 |
|
|
|
28 |
|
|
! Local: |
29 |
|
|
REAL h(size(flux_t), klev) ! (knon, klev) enthalpie potentielle |
30 |
|
|
INTEGER k |
31 |
|
|
REAL local_q(size(flux_t), klev) ! (knon, klev) |
32 |
|
|
|
33 |
|
|
!---------------------------------------------------------------------- |
34 |
|
|
|
35 |
guez |
299 |
h(:, 1) = ch(:, 1) + dh(:, 1) * flux_t * dtphys |
36 |
|
|
local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtphys |
37 |
guez |
296 |
|
38 |
|
|
DO k = 2, klev |
39 |
|
|
h(:, k) = ch(:, k) + dh(:, k) * h(:, k - 1) |
40 |
|
|
local_q(:, k) = cq(:, k) + dq(:, k) * local_q(:, k - 1) |
41 |
|
|
ENDDO |
42 |
|
|
|
43 |
|
|
d_t = h / pkf / RCPD - t |
44 |
|
|
d_q = local_q - q |
45 |
|
|
|
46 |
|
|
end subroutine climb_hq_up |
47 |
|
|
|
48 |
|
|
end module climb_hq_up_m |