1 | SUBROUTINE nflxtr(pdtime,pmfu,pmfd,pen_u,pde_u,pen_d,pde_d, |
---|
2 | . pplay,paprs,x,dx) |
---|
3 | IMPLICIT NONE |
---|
4 | c===================================================================== |
---|
5 | c Objet : Melange convectif de traceurs a partir des flux de masse |
---|
6 | c Date : 13/12/1996 -- 13/01/97 |
---|
7 | c Auteur: O. Boucher (LOA) sur inspiration de Z. X. Li (LMD), |
---|
8 | c Brinkop et Sausen (1996) et Boucher et al. (1996). |
---|
9 | c ATTENTION : meme si cette routine se veut la plus generale possible, |
---|
10 | c elle a herite de certaines notations et conventions du |
---|
11 | c schema de Tiedtke (1993). |
---|
12 | c --En particulier, les couches sont numerotees de haut en bas !!! |
---|
13 | c Ceci est valable pour les flux |
---|
14 | c mais pas pour les entrees x, pplay, paprs !!!! |
---|
15 | c --pmfu est positif, pmfd est negatif |
---|
16 | c --Tous les flux d'entrainements et de detrainements sont positifs |
---|
17 | c contrairement au schema de Tiedtke d'ou les changements de signe!!!! |
---|
18 | c===================================================================== |
---|
19 | c |
---|
20 | #include "dimensions.h" |
---|
21 | #include "dimphy.h" |
---|
22 | #include "YOMCST.h" |
---|
23 | #include "YOECUMF.h" |
---|
24 | c |
---|
25 | REAL pdtime |
---|
26 | c--les flux sont definis au 1/2 niveaux |
---|
27 | c--pmfu(klev+1) et pmfd(klev+1) sont implicitement nuls |
---|
28 | REAL pmfu(klon,klev) ! flux de masse dans le panache montant |
---|
29 | REAL pmfd(klon,klev) ! flux de masse dans le panache descendant |
---|
30 | REAL pen_u(klon,klev) ! flux entraine dans le panache montant |
---|
31 | REAL pde_u(klon,klev) ! flux detraine dans le panache montant |
---|
32 | REAL pen_d(klon,klev) ! flux entraine dans le panache descendant |
---|
33 | REAL pde_d(klon,klev) ! flux detraine dans le panache descendant |
---|
34 | |
---|
35 | REAL pplay(klon,klev) ! pression aux couches (bas en haut) |
---|
36 | REAL paprs(klon,klev+1) ! pression aux 1/2 couches (bas en haut) |
---|
37 | REAL x(klon,klev) ! q de traceur (bas en haut) |
---|
38 | REAL dx(klon,klev) ! tendance de traceur (bas en haut) |
---|
39 | c |
---|
40 | c--flux convectifs mais en variables locales |
---|
41 | REAL zmfu(klon,klev+1) |
---|
42 | REAL zmfd(klon,klev+1) |
---|
43 | REAL zen_u(klon,klev) |
---|
44 | REAL zde_u(klon,klev) |
---|
45 | REAL zen_d(klon,klev) |
---|
46 | REAL zde_d(klon,klev) |
---|
47 | real zmfe |
---|
48 | c |
---|
49 | c--variables locales |
---|
50 | c--les flux de x sont definis aux 1/2 niveaux |
---|
51 | c--xu et xd sont definis aux niveaux complets |
---|
52 | REAL xu(klon,klev) ! q de traceurs dans le panache montant |
---|
53 | REAL xd(klon,klev) ! q de traceurs dans le panache descendant |
---|
54 | REAL zmfux(klon,klev+1) ! flux de x dans le panache montant |
---|
55 | REAL zmfdx(klon,klev+1) ! flux de x dans le panache descendant |
---|
56 | REAL zmfex(klon,klev+1) ! flux de x dans l'environnement |
---|
57 | INTEGER i, k |
---|
58 | REAL zmfmin |
---|
59 | PARAMETER (zmfmin=1.E-10) |
---|
60 | |
---|
61 | c ========================================= |
---|
62 | c |
---|
63 | c |
---|
64 | c Extension des flux UP et DN sur klev+1 niveaux |
---|
65 | c ========================================= |
---|
66 | do k=1,klev |
---|
67 | do i=1,klon |
---|
68 | zmfu(i,k)=pmfu(i,k) |
---|
69 | zmfd(i,k)=pmfd(i,k) |
---|
70 | enddo |
---|
71 | enddo |
---|
72 | do i=1,klon |
---|
73 | zmfu(i,klev+1)=0. |
---|
74 | zmfd(i,klev+1)=0. |
---|
75 | enddo |
---|
76 | |
---|
77 | c--modif pour diagnostiquer les detrainements |
---|
78 | c ========================================= |
---|
79 | c on privilegie l'ajustement de l'entrainement dans l'ascendance. |
---|
80 | |
---|
81 | do k=1, klev |
---|
82 | do i=1, klon |
---|
83 | zen_d(i,k)=pen_d(i,k) |
---|
84 | zde_u(i,k)=pde_u(i,k) |
---|
85 | zde_d(i,k) =-zmfd(i,k+1)+zmfd(i,k)+zen_d(i,k) |
---|
86 | zen_u(i,k) = zmfu(i,k+1)-zmfu(i,k)+zde_u(i,k) |
---|
87 | enddo |
---|
88 | enddo |
---|
89 | c |
---|
90 | c--calcul des flux dans le panache montant |
---|
91 | c ========================================= |
---|
92 | c |
---|
93 | c Dans la premiere couche, on prend q comme valeur de qu |
---|
94 | c |
---|
95 | do i=1, klon |
---|
96 | zmfux(i,1)=0.0 |
---|
97 | enddo |
---|
98 | c |
---|
99 | c Autres couches |
---|
100 | do k=1,klev |
---|
101 | do i=1, klon |
---|
102 | if ((zmfu(i,k+1)+zde_u(i,k)).lt.zmfmin) THEN |
---|
103 | xu(i,k)=x(i,k) |
---|
104 | else |
---|
105 | xu(i,k)=(zmfux(i,k)+zen_u(i,k)*x(i,k)) |
---|
106 | s /(zmfu(i,k+1)+zde_u(i,k)) |
---|
107 | endif |
---|
108 | zmfux(i,k+1)=zmfu(i,k+1)*xu(i,k) |
---|
109 | enddo |
---|
110 | enddo |
---|
111 | c |
---|
112 | c--calcul des flux dans le panache descendant |
---|
113 | c ========================================= |
---|
114 | c |
---|
115 | do i=1, klon |
---|
116 | zmfdx(i,klev+1)=0.0 |
---|
117 | enddo |
---|
118 | c |
---|
119 | do k=klev,1,-1 |
---|
120 | do i=1, klon |
---|
121 | if ((zde_d(i,k)-zmfd(i,k)).lt.zmfmin) THEN |
---|
122 | xd(i,k)=x(i,k) |
---|
123 | else |
---|
124 | xd(i,k)=(zmfdx(i,k+1)-zen_d(i,k)*x(i,k)) / |
---|
125 | . (zmfd(i,k)-zde_d(i,k)) |
---|
126 | endif |
---|
127 | zmfdx(i,k)=zmfd(i,k)*xd(i,k) |
---|
128 | enddo |
---|
129 | enddo |
---|
130 | c |
---|
131 | c--introduction du flux de retour dans l'environnement |
---|
132 | c ========================================= |
---|
133 | c |
---|
134 | do k=2, klev |
---|
135 | do i=1, klon |
---|
136 | zmfe=-zmfu(i,k)-zmfd(i,k) |
---|
137 | if (zmfe.le.0.) then |
---|
138 | zmfex(i,k)= zmfe*x(i,k) |
---|
139 | else |
---|
140 | zmfex(i,k)= zmfe*x(i,k-1) |
---|
141 | endif |
---|
142 | enddo |
---|
143 | enddo |
---|
144 | |
---|
145 | do i=1, klon |
---|
146 | zmfex(i,1)=0. |
---|
147 | zmfex(i,klev+1)=0. |
---|
148 | enddo |
---|
149 | c |
---|
150 | c--calcul final des tendances |
---|
151 | c |
---|
152 | do k=1, klev |
---|
153 | do i=1, klon |
---|
154 | dx(i,k)=RG/(paprs(i,k)-paprs(i,k+1))*pdtime* |
---|
155 | . ( zmfux(i,k) - zmfux(i,k+1) + |
---|
156 | . zmfdx(i,k) - zmfdx(i,k+1) + |
---|
157 | . zmfex(i,k) - zmfex(i,k+1) ) |
---|
158 | enddo |
---|
159 | enddo |
---|
160 | c |
---|
161 | return |
---|
162 | end |
---|