| 1 |
module climb_hq_up_m |
| 2 |
|
| 3 |
IMPLICIT none |
| 4 |
|
| 5 |
contains |
| 6 |
|
| 7 |
subroutine climb_hq_up(d_t, d_q, cq, dq, ch, dh, flux_t, flux_q, pkf, t, q) |
| 8 |
|
| 9 |
use comconst, only: dtphys |
| 10 |
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 |
h(:, 1) = ch(:, 1) + dh(:, 1) * flux_t * dtphys |
| 36 |
local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtphys |
| 37 |
|
| 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 |
Parent Directory
| 
Revision Log