phylmd/Interface_surf/climb_hq_up.f90

module climb_hq_up_m

  IMPLICIT none

contains

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

    use comconst, only: dtphys
    USE dimphy, ONLY: klev
    USE suphec_m, ONLY: rcpd

    REAL, intent(out):: d_t(:, :) ! (knon, klev) variation of air temperature 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):: 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 * dtphys
    local_q(:, 1) = cq(:, 1) + dq(:, 1) * flux_q * dtphys

    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, 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) variation of air temperature 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