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	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	guez	341	REAL, intent(out):: d_t(:, :) ! (knon, klev) variation of air temperature t
14	guez	296	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