phylmd/Interface_surf/climb_hq_up.f

module climb_hq_up_m

  IMPLICIT none

contains

  subroutine climb_hq_up(d_t, d_q, cq, dq, ch, dh, flux_t, flux_q, dtime, pkf, &
       t, q)

    USE dimphy, ONLY: klev
    USE suphec_m, ONLY: rcpd

    REAL, intent(out):: d_t(:, :) ! (knon, klev) incrementation de "t"
    REAL, intent(out):: d_q(:, :) ! (knon, klev) incrementation de "q"
    REAL, intent(in), dimension(:, :):: cq, dq, ch, dh ! (knon, klev)

    REAL, intent(in):: flux_t(:) ! (knon)
    ! (diagnostic) flux de chaleur sensible (Cp T) à la surface,
    ! positif vers le bas, W / m2

    REAL, intent(out):: flux_q(:) ! (knon)
    ! flux de la vapeur d'eau à la surface, en kg / (m**2 s)

    REAL, intent(in):: dtime ! intervalle du temps (s)
    REAL, intent(in):: pkf(:, :) ! (knon, klev)
    REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K)
    REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg)

    ! Local:
    REAL h(size(flux_t), klev) ! (knon, klev) enthalpie potentielle
    INTEGER k
    REAL local_q(size(flux_t), klev) ! (knon, klev)

    !----------------------------------------------------------------------

    h(:, 1) = ch(:, 1) + dh(:, 1) * flux_t * dtime
    local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtime

    DO k = 2, klev
       h(:, k) = ch(:, k) + dh(:, k) * h(:, k - 1)
       local_q(:, k) = cq(:, k) + dq(:, k) * local_q(:, k - 1)
    ENDDO

    d_t = h / pkf / RCPD - t
    d_q = local_q - q

  end subroutine climb_hq_up
  
end module climb_hq_up_m
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, dtime, pkf, &
8	t, q)
9
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):: dtime ! intervalle du temps (s)
25	REAL, intent(in):: pkf(:, :) ! (knon, klev)
26	REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K)
27	REAL, intent(in):: q(:, :) ! (knon, klev) humidite specifique (kg / kg)
28
29	! Local:
30	REAL h(size(flux_t), klev) ! (knon, klev) enthalpie potentielle
31	INTEGER k
32	REAL local_q(size(flux_t), klev) ! (knon, klev)
33
34	!----------------------------------------------------------------------
35
36	h(:, 1) = ch(:, 1) + dh(:, 1) * flux_t * dtime
37	local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtime
38
39	DO k = 2, klev
40	h(:, k) = ch(:, k) + dh(:, k) * h(:, k - 1)
41	local_q(:, k) = cq(:, k) + dq(:, k) * local_q(:, k - 1)
42	ENDDO
43
44	d_t = h / pkf / RCPD - t
45	d_q = local_q - q
46
47	end subroutine climb_hq_up
48
49	end module climb_hq_up_m