Sources/phylmd/cvltr.f

module cvltr_m

  IMPLICIT NONE

contains

  SUBROUTINE cvltr(pdtime, da, phi, mp, paprs, x, upd, dnd, dx)

    ! From LMDZ4/libf/phylmd/cvltr.F, version 1.1 2005/04/15 12:36:17

    USE dimphy, ONLY: klev, klon
    USE suphec_m, ONLY: rg

    ! Objet : convection des traceurs / Kerry Emanuel
    ! Authors: M.-A. Filiberti and J.-Y. Grandpeix

    REAL, intent(in):: pdtime
    real, intent(in):: da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
    REAL, intent(in):: paprs(klon, klev + 1) ! pression aux 1/2 couches
    REAL, intent(in):: x(klon, klev) ! q de traceur (bas en haut)
    REAL, intent(in):: upd(klon, klev) ! saturated updraft mass flux
    REAL, intent(in):: dnd(klon, klev) ! saturated downdraft mass flux
    REAL, intent(out):: dx(klon, klev) ! tendance de traceur (bas en haut)

    ! Local:
    real zed(klon, klev), zmd(klon, klev, klev)
    real za(klon, klev, klev)
    real zmfd(klon, klev), zmfa(klon, klev)
    real zmfp(klon, klev), zmfu(klon, klev)
    integer i, k, j

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

    ! calcul des tendances liees au downdraft

    zed = 0.
    zmfd = 0.
    zmfa = 0.
    zmfu = 0.
    zmfp = 0.
    zmd = 0.
    za = 0.
    ! entrainement
    do k = 1, klev - 1
       do i = 1, klon
          zed(i, k) = max(0., mp(i, k) - mp(i, k + 1))
       end do
    end do

    ! calcul de la matrice d echange
    ! matrice de distribution de la masse entrainee en k

    do k = 1, klev
       do i = 1, klon
          zmd(i, k, k) = zed(i, k)
       end do
    end do
    do k = 2, klev
       do j = k - 1, 1, - 1
          do i = 1, klon
             if(mp(i, j + 1) /= 0) then
                zmd(i, j, k) = zmd(i, j + 1, k) * min(1., mp(i, j)/mp(i, j + 1))
             endif
          end do
       end do
    end do
    do k = 1, klev
       do j = 1, klev - 1
          do i = 1, klon
             za(i, j, k) = max(0., zmd(i, j + 1, k) - zmd(i, j, k))
          end do
       end do
    end do

    ! rajout du terme lie a l'ascendance induite

    do j = 2, klev
       do i = 1, klon
          za(i, j, j - 1) = za(i, j, j - 1) + mp(i, j)
       end do
    end do

    ! tendances

    do k = 1, klev
       do j = 1, klev
          do i = 1, klon
             zmfd(i, j) = zmfd(i, j) + za(i, j, k) * (x(i, k) - x(i, j))
          end do
       end do
    end do

    ! calcul des tendances liees aux flux satures

    do j = 1, klev
       do i = 1, klon
          zmfa(i, j) = da(i, j) * (x(i, 1) - x(i, j))
       end do
    end do
    do k = 1, klev
       do j = 1, klev
          do i = 1, klon
             zmfp(i, j) = zmfp(i, j) + phi(i, j, k) * (x(i, k) - x(i, j))
          end do
       end do
    end do
    do j = 1, klev - 1
       do i = 1, klon
          zmfu(i, j) = max(0., upd(i, j + 1) + dnd(i, j + 1)) &
               * (x(i, j + 1) - x(i, j))
       end do
    end do
    do j = 2, klev
       do i = 1, klon
          zmfu(i, j) = zmfu(i, j) &
               + min(0., upd(i, j) + dnd(i, j)) * (x(i, j) - x(i, j - 1))
       end do
    end do

    ! calcul final des tendances

    do k = 1, klev
       do i = 1, klon
          dx(i, k) = (zmfd(i, k) + zmfu(i, k) &
               + zmfa(i, k) + zmfp(i, k)) * pdtime &
               * RG/(paprs(i, k) - paprs(i, k + 1))
       enddo
    enddo

  end SUBROUTINE cvltr

end module cvltr_m
1	guez	120	module cvltr_m
2	guez	63
3	guez	212	IMPLICIT NONE
4	guez	63
5	guez	120	contains
6	guez	63
7	guez	212	SUBROUTINE cvltr(pdtime, da, phi, mp, paprs, x, upd, dnd, dx)
8	guez	3
9	guez	212	! From LMDZ4/libf/phylmd/cvltr.F, version 1.1 2005/04/15 12:36:17
10	guez	63
11	guez	120	USE dimphy, ONLY: klev, klon
12			USE suphec_m, ONLY: rg
13
14	guez	212	! Objet : convection des traceurs / Kerry Emanuel
15			! Authors: M.-A. Filiberti and J.-Y. Grandpeix
16
17	guez	120	REAL, intent(in):: pdtime
18	guez	212	real, intent(in):: da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
19			REAL, intent(in):: paprs(klon, klev + 1) ! pression aux 1/2 couches
20			REAL, intent(in):: x(klon, klev) ! q de traceur (bas en haut)
21			REAL, intent(in):: upd(klon, klev) ! saturated updraft mass flux
22			REAL, intent(in):: dnd(klon, klev) ! saturated downdraft mass flux
23			REAL, intent(out):: dx(klon, klev) ! tendance de traceur (bas en haut)
24
25			! Local:
26			real zed(klon, klev), zmd(klon, klev, klev)
27			real za(klon, klev, klev)
28			real zmfd(klon, klev), zmfa(klon, klev)
29			real zmfp(klon, klev), zmfu(klon, klev)
30			integer i, k, j
31
32			!------------------------------------------------------------
33
34	guez	120	! calcul des tendances liees au downdraft
35	guez	212
36			zed = 0.
37			zmfd = 0.
38			zmfa = 0.
39			zmfu = 0.
40			zmfp = 0.
41			zmd = 0.
42			za = 0.
43	guez	120	! entrainement
44	guez	212	do k = 1, klev - 1
45			do i = 1, klon
46			zed(i, k) = max(0., mp(i, k) - mp(i, k + 1))
47	guez	120	end do
48			end do
49	guez	212
50	guez	120	! calcul de la matrice d echange
51			! matrice de distribution de la masse entrainee en k
52	guez	212
53			do k = 1, klev
54			do i = 1, klon
55			zmd(i, k, k) = zed(i, k)
56	guez	120	end do
57			end do
58	guez	212	do k = 2, klev
59			do j = k - 1, 1, - 1
60			do i = 1, klon
61			if(mp(i, j + 1) /= 0) then
62			zmd(i, j, k) = zmd(i, j + 1, k) * min(1., mp(i, j)/mp(i, j + 1))
63	guez	120	endif
64			end do
65			end do
66			end do
67	guez	212	do k = 1, klev
68			do j = 1, klev - 1
69			do i = 1, klon
70			za(i, j, k) = max(0., zmd(i, j + 1, k) - zmd(i, j, k))
71	guez	120	end do
72			end do
73			end do
74	guez	212
75			! rajout du terme lie a l'ascendance induite
76
77			do j = 2, klev
78			do i = 1, klon
79			za(i, j, j - 1) = za(i, j, j - 1) + mp(i, j)
80	guez	120	end do
81			end do
82	guez	212
83	guez	120	! tendances
84	guez	212
85			do k = 1, klev
86			do j = 1, klev
87			do i = 1, klon
88			zmfd(i, j) = zmfd(i, j) + za(i, j, k) * (x(i, k) - x(i, j))
89	guez	120	end do
90			end do
91			end do
92	guez	212
93	guez	120	! calcul des tendances liees aux flux satures
94	guez	212
95			do j = 1, klev
96			do i = 1, klon
97			zmfa(i, j) = da(i, j) * (x(i, 1) - x(i, j))
98	guez	120	end do
99			end do
100	guez	212	do k = 1, klev
101			do j = 1, klev
102			do i = 1, klon
103			zmfp(i, j) = zmfp(i, j) + phi(i, j, k) * (x(i, k) - x(i, j))
104	guez	120	end do
105			end do
106			end do
107	guez	212	do j = 1, klev - 1
108			do i = 1, klon
109			zmfu(i, j) = max(0., upd(i, j + 1) + dnd(i, j + 1)) &
110			* (x(i, j + 1) - x(i, j))
111	guez	120	end do
112			end do
113	guez	212	do j = 2, klev
114			do i = 1, klon
115			zmfu(i, j) = zmfu(i, j) &
116			+ min(0., upd(i, j) + dnd(i, j)) * (x(i, j) - x(i, j - 1))
117	guez	120	end do
118			end do
119
120	guez	212	! calcul final des tendances
121
122			do k = 1, klev
123			do i = 1, klon
124			dx(i, k) = (zmfd(i, k) + zmfu(i, k) &
125			+ zmfa(i, k) + zmfp(i, k)) * pdtime &
126			* RG/(paprs(i, k) - paprs(i, k + 1))
127	guez	120	enddo
128			enddo
129
130			end SUBROUTINE cvltr
131
132			end module cvltr_m