Directory: | ./ |
---|---|
File: | phys/cv3p_mixing.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 238 | 258 | 92.2% |
Branches: | 276 | 316 | 87.3% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | 81010721 | SUBROUTINE cv3p_mixing(nloc, ncum, nd, na, ntra, icb, nk, inb, & | |
2 | 240 | ph, t, rr, rs, u, v, tra, h, lv, lf, frac, qta, & | |
3 | 240 | unk, vnk, hp, tv, tvp, ep, clw, sig, & | |
4 | 240 | Ment, Qent, hent, uent, vent, nent, & | |
5 | Sigij, elij, supmax, Ments, Qents, traent) | ||
6 | ! ************************************************************** | ||
7 | ! * | ||
8 | ! CV3P_MIXING : compute mixed draught properties and, * | ||
9 | ! within a scaling factor, mixed draught * | ||
10 | ! mass fluxes. * | ||
11 | ! written by : VTJ Philips,JY Grandpeix, 21/05/2003, 09.14.15* | ||
12 | ! modified by : * | ||
13 | ! ************************************************************** | ||
14 | |||
15 | USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level | ||
16 | USE ioipsl_getin_p_mod, ONLY: getin_p | ||
17 | USE add_phys_tend_mod, ONLY: fl_cor_ebil | ||
18 | |||
19 | IMPLICIT NONE | ||
20 | |||
21 | include "cvthermo.h" | ||
22 | include "cv3param.h" | ||
23 | include "YOMCST2.h" | ||
24 | include "cvflag.h" | ||
25 | |||
26 | !inputs: | ||
27 | INTEGER, INTENT (IN) :: ncum, nd, na | ||
28 | INTEGER, INTENT (IN) :: ntra, nloc | ||
29 | INTEGER, DIMENSION (nloc), INTENT (IN) :: icb, inb, nk | ||
30 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: sig | ||
31 | REAL, DIMENSION (nloc), INTENT (IN) :: unk, vnk | ||
32 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: qta | ||
33 | REAL, DIMENSION (nloc, nd+1), INTENT (IN) :: ph | ||
34 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: t, rr, rs | ||
35 | REAL, DIMENSION (nloc, nd), INTENT (IN) :: u, v | ||
36 | REAL, DIMENSION (nloc, nd, ntra), INTENT (IN) :: tra ! input of convect3 | ||
37 | REAL, DIMENSION (nloc, na), INTENT (IN) :: lv | ||
38 | REAL, DIMENSION (nloc, na), INTENT (IN) :: lf | ||
39 | REAL, DIMENSION (nloc, na), INTENT (IN) :: frac !ice fraction in condensate | ||
40 | REAL, DIMENSION (nloc, na), INTENT (IN) :: h !liquid water static energy of environment | ||
41 | REAL, DIMENSION (nloc, na), INTENT (IN) :: hp !liquid water static energy of air shed from adiab. asc. | ||
42 | REAL, DIMENSION (nloc, na), INTENT (IN) :: tv, tvp | ||
43 | REAL, DIMENSION (nloc, na), INTENT (IN) :: ep, clw | ||
44 | |||
45 | !outputs: | ||
46 | REAL, DIMENSION (nloc, na, na), INTENT (OUT) :: Ment, Qent | ||
47 | REAL, DIMENSION (nloc, na, na), INTENT (OUT) :: uent, vent | ||
48 | REAL, DIMENSION (nloc, na, na), INTENT (OUT) :: Sigij, elij | ||
49 | REAL, DIMENSION (nloc, na), INTENT (OUT) :: supmax ! Highest mixing fraction of mixed | ||
50 | ! updraughts with the sign of (h-hp) | ||
51 | REAL, DIMENSION (nloc, nd, nd, ntra), INTENT (OUT) :: traent | ||
52 | REAL, DIMENSION (nloc, nd, nd), INTENT (OUT) :: Ments, Qents | ||
53 | REAL, DIMENSION (nloc, nd, nd), INTENT (OUT) :: hent | ||
54 | INTEGER, DIMENSION (nloc, nd), INTENT (OUT) :: nent | ||
55 | |||
56 | !local variables: | ||
57 | INTEGER i, j, k, il, im, jm | ||
58 | INTEGER num1, num2 | ||
59 | REAL :: rti, bf2, anum, denom, dei, altem, cwat, stemp | ||
60 | REAL :: alt, delp, delm | ||
61 | 480 | REAL, DIMENSION (nloc) :: Qmixmax, Rmixmax, sqmrmax | |
62 | 480 | REAL, DIMENSION (nloc) :: Qmixmin, Rmixmin, sqmrmin | |
63 | 480 | REAL, DIMENSION (nloc) :: signhpmh | |
64 | 480 | REAL, DIMENSION (nloc) :: Sx | |
65 | REAL :: Scrit2 | ||
66 | 480 | REAL, DIMENSION (nloc) :: Smid, Sjmin, Sjmax | |
67 | 480 | REAL, DIMENSION (nloc) :: Sbef, sup, smin | |
68 | 480 | REAL, DIMENSION (nloc) :: ASij, ASij_inv, smax, Scrit | |
69 | 480 | REAL, DIMENSION (nloc, nd, nd) :: Sij | |
70 | 480 | REAL, DIMENSION (nloc, nd) :: csum | |
71 | REAL :: awat | ||
72 | REAL :: cpm !Mixed draught heat capacity | ||
73 | REAL :: Tm !Mixed draught temperature | ||
74 | 480 | LOGICAL, DIMENSION (nloc) :: lwork | |
75 | |||
76 | REAL amxupcrit, df, ff | ||
77 | INTEGER nstep | ||
78 | |||
79 | INTEGER,SAVE :: igout=1 | ||
80 | !$OMP THREADPRIVATE(igout) | ||
81 | |||
82 | ! -- Mixing probability distribution functions | ||
83 | |||
84 | REAL Qcoef1, Qcoef2, QFF, QFFF, Qmix, Rmix, Qmix1, Rmix1, Qmix2, Rmix2, F | ||
85 | |||
86 | Qcoef1(F) = tanh(F/gammas) | ||
87 | Qcoef2(F) = (tanh(F/gammas)+gammas*log(cosh((1.-F)/gammas)/cosh(F/gammas))) | ||
88 | QFF(F) = max(min(F,1.), 0.) | ||
89 | QFFf(F) = min(QFF(F), scut) | ||
90 | Qmix1(F) = (tanh((QFF(F)-Fmax)/gammas)+Qcoef1max)/Qcoef2max | ||
91 | Rmix1(F) = (gammas*log(cosh((QFF(F)-Fmax)/gammas))+QFF(F)*Qcoef1max)/Qcoef2max | ||
92 | Qmix2(F) = -log(1.-QFFf(F))/scut | ||
93 | Rmix2(F) = (QFFf(F)+(1.-QFF(F))*log(1.-QFFf(F)))/scut | ||
94 | Qmix(F) = qqa1*Qmix1(F) + qqa2*Qmix2(F) | ||
95 | Rmix(F) = qqa1*Rmix1(F) + qqa2*Rmix2(F) | ||
96 | |||
97 | INTEGER, SAVE :: ifrst | ||
98 | DATA ifrst/0/ | ||
99 | !$OMP THREADPRIVATE(ifrst) | ||
100 | |||
101 | |||
102 | ! ===================================================================== | ||
103 | ! --- INITIALIZE VARIOUS ARRAYS USED IN THE COMPUTATIONS | ||
104 | ! ===================================================================== | ||
105 | |||
106 | ! -- Initialize mixing PDF coefficients | ||
107 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 239 times.
|
240 | IF (ifrst==0) THEN |
108 | 1 | ifrst = 1 | |
109 | 1 | Qcoef1max = Qcoef1(Fmax) | |
110 | 1 | Qcoef2max = Qcoef2(Fmax) | |
111 | !<jyg | ||
112 | 1 | print*, 'fmax, gammas, qqa1, qqa2, Qcoef1max, Qcoef2max ', & | |
113 | 2 | fmax, gammas, qqa1, qqa2, Qcoef1max, Qcoef2max | |
114 | !>jyg | ||
115 | ! | ||
116 | END IF | ||
117 | |||
118 | |||
119 | ! ori do 360 i=1,ncum*nlp | ||
120 |
2/2✓ Branch 0 taken 6480 times.
✓ Branch 1 taken 240 times.
|
6720 | DO j = 1, nl |
121 |
2/2✓ Branch 0 taken 3431619 times.
✓ Branch 1 taken 6480 times.
|
3438339 | DO i = 1, ncum |
122 | 3438099 | nent(i, j) = 0 | |
123 | ! in convect3, m is computed in cv3_closure | ||
124 | ! ori m(i,1)=0.0 | ||
125 | END DO | ||
126 | END DO | ||
127 | |||
128 | ! ori do 400 k=1,nlp | ||
129 | ! ori do 390 j=1,nlp | ||
130 |
2/2✓ Branch 0 taken 6480 times.
✓ Branch 1 taken 240 times.
|
6720 | DO j = 1, nl |
131 |
2/2✓ Branch 0 taken 174960 times.
✓ Branch 1 taken 6480 times.
|
181680 | DO k = 1, nl |
132 |
2/2✓ Branch 0 taken 92653713 times.
✓ Branch 1 taken 174960 times.
|
92835153 | DO i = 1, ncum |
133 | 92653713 | Qent(i, k, j) = rr(i, j) | |
134 | 92653713 | uent(i, k, j) = u(i, j) | |
135 | 92653713 | vent(i, k, j) = v(i, j) | |
136 | 92653713 | elij(i, k, j) = 0.0 | |
137 | 92828673 | hent(i, k, j) = 0.0 | |
138 | !AC! Ment(i,k,j)=0.0 | ||
139 | !AC! Sij(i,k,j)=0.0 | ||
140 | END DO | ||
141 | END DO | ||
142 | END DO | ||
143 | |||
144 | !AC! | ||
145 |
6/6✓ Branch 0 taken 9360 times.
✓ Branch 1 taken 240 times.
✓ Branch 2 taken 365040 times.
✓ Branch 3 taken 9360 times.
✓ Branch 4 taken 193314537 times.
✓ Branch 5 taken 365040 times.
|
193689177 | Ment(1:ncum, 1:nd, 1:nd) = 0.0 |
146 |
6/6✓ Branch 0 taken 9360 times.
✓ Branch 1 taken 240 times.
✓ Branch 2 taken 365040 times.
✓ Branch 3 taken 9360 times.
✓ Branch 4 taken 193314537 times.
✓ Branch 5 taken 365040 times.
|
193689177 | Sij(1:ncum, 1:nd, 1:nd) = 0.0 |
147 | !AC! | ||
148 | !ym | ||
149 |
6/6✓ Branch 0 taken 9360 times.
✓ Branch 1 taken 240 times.
✓ Branch 2 taken 365040 times.
✓ Branch 3 taken 9360 times.
✓ Branch 4 taken 193314537 times.
✓ Branch 5 taken 365040 times.
|
193689177 | Sigij(1:ncum, 1:nd, 1:nd) = 0.0 |
150 | !ym | ||
151 | |||
152 | !jyg! DO k = 1, ntra | ||
153 | !jyg! DO j = 1, nd ! instead nlp | ||
154 | !jyg! DO i = 1, nd ! instead nlp | ||
155 | !jyg! DO il = 1, ncum | ||
156 | !jyg! traent(il, i, j, k) = tra(il, j, k) | ||
157 | !jyg! END DO | ||
158 | !jyg! END DO | ||
159 | !jyg! END DO | ||
160 | !jyg! END DO | ||
161 | |||
162 | ! ===================================================================== | ||
163 | ! --- CALCULATE ENTRAINED AIR MASS FLUX (Ment), TOTAL WATER MIXING | ||
164 | ! --- RATIO (QENT), TOTAL CONDENSED WATER (elij), AND MIXING | ||
165 | ! --- FRACTION (Sij) | ||
166 | ! ===================================================================== | ||
167 | |||
168 |
2/2✓ Branch 0 taken 6240 times.
✓ Branch 1 taken 240 times.
|
6480 | DO i = minorig + 1, nl |
169 | |||
170 |
1/2✓ Branch 0 taken 6240 times.
✗ Branch 1 not taken.
|
6240 | IF (ok_entrain) THEN |
171 |
2/2✓ Branch 0 taken 168480 times.
✓ Branch 1 taken 6240 times.
|
174720 | DO j = minorig, nl |
172 |
2/2✓ Branch 0 taken 89222094 times.
✓ Branch 1 taken 168480 times.
|
89396814 | DO il = 1, ncum |
173 | IF ((i>=icb(il)) .AND. (i<=inb(il)) .AND. (j>=(icb(il)-1)) & | ||
174 |
8/8✓ Branch 0 taken 78601536 times.
✓ Branch 1 taken 10620558 times.
✓ Branch 2 taken 30907008 times.
✓ Branch 3 taken 47694528 times.
✓ Branch 4 taken 27574473 times.
✓ Branch 5 taken 3332535 times.
✓ Branch 6 taken 14357294 times.
✓ Branch 7 taken 13217179 times.
|
89390574 | .AND. (j<=inb(il))) THEN |
175 | |||
176 | !! rti = qnk(il) - ep(il, i)*clw(il, i) | ||
177 | 14357294 | rti = qta(il,i-1) - ep(il, i)*clw(il, i) | |
178 | 14357294 | bf2 = 1. + lv(il, j)*lv(il, j)*rs(il, j)/(rrv*t(il,j)*t(il,j)*cpd) | |
179 | !jyg(from aj)< | ||
180 |
1/2✓ Branch 0 taken 14357294 times.
✗ Branch 1 not taken.
|
14357294 | IF (cvflag_ice) THEN |
181 | ! print*,cvflag_ice,'cvflag_ice dans do 700' | ||
182 |
2/2✓ Branch 0 taken 5723139 times.
✓ Branch 1 taken 8634155 times.
|
14357294 | IF (t(il,j)<=263.15) THEN |
183 | bf2 = 1. + (lf(il,j)+lv(il,j))*(lv(il,j)+frac(il,j)* & | ||
184 | 5723139 | lf(il,j))*rs(il, j)/(rrv*t(il,j)*t(il,j)*cpd) | |
185 | END IF | ||
186 | END IF | ||
187 | !>jyg | ||
188 | 14357294 | anum = h(il, j) - hp(il, i) + (cpv-cpd)*t(il, j)*(rti-rr(il,j)) | |
189 | 14357294 | denom = h(il, i) - hp(il, i) + (cpd-cpv)*(rr(il,i)-rti)*t(il, j) | |
190 | dei = denom | ||
191 |
2/2✓ Branch 0 taken 114 times.
✓ Branch 1 taken 14357180 times.
|
14357294 | IF (abs(dei)<0.01) dei = 0.01 |
192 | 14357294 | Sij(il, i, j) = anum/dei | |
193 | 14357294 | Sij(il, i, i) = 1.0 | |
194 | 14357294 | altem = Sij(il, i, j)*rr(il, i) + (1.-Sij(il,i,j))*rti - rs(il, j) | |
195 | 14357294 | altem = altem/bf2 | |
196 | 14357294 | cwat = clw(il, j)*(1.-ep(il,j)) | |
197 | stemp = Sij(il, i, j) | ||
198 |
8/8✓ Branch 0 taken 9791411 times.
✓ Branch 1 taken 4565883 times.
✓ Branch 2 taken 1793584 times.
✓ Branch 3 taken 7997827 times.
✓ Branch 4 taken 565250 times.
✓ Branch 5 taken 1228334 times.
✓ Branch 6 taken 6031093 times.
✓ Branch 7 taken 7097867 times.
|
14357294 | IF ((stemp<0.0 .OR. stemp>1.0 .OR. altem>cwat) .AND. j>i) THEN |
199 | !jyg(from aj)< | ||
200 |
1/2✓ Branch 0 taken 6031093 times.
✗ Branch 1 not taken.
|
6031093 | IF (cvflag_ice) THEN |
201 | 6031093 | anum = anum - (lv(il,j)+frac(il,j)*lf(il,j))*(rti-rs(il,j)-cwat*bf2) | |
202 | 6031093 | denom = denom + (lv(il,j)+frac(il,j)*lf(il,j))*(rr(il,i)-rti) | |
203 | ELSE | ||
204 | ✗ | anum = anum - lv(il, j)*(rti-rs(il,j)-cwat*bf2) | |
205 | ✗ | denom = denom + lv(il, j)*(rr(il,i)-rti) | |
206 | END IF | ||
207 | !>jyg | ||
208 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 6031092 times.
|
6031093 | IF (abs(denom)<0.01) denom = 0.01 |
209 | 6031093 | Sij(il, i, j) = anum/denom | |
210 | 6031093 | altem = Sij(il, i, j)*rr(il, i) + (1.-Sij(il,i,j))*rti - rs(il, j) | |
211 | 6031093 | altem = altem - (bf2-1.)*cwat | |
212 | END IF | ||
213 |
2/2✓ Branch 0 taken 12499941 times.
✓ Branch 1 taken 1857353 times.
|
14357294 | IF (Sij(il,i,j)>0.0) THEN |
214 | !!! Ment(il,i,j)=m(il,i) | ||
215 | 12499941 | Ment(il, i, j) = 1. | |
216 | elij(il, i, j) = altem | ||
217 | 12499941 | elij(il, i, j) = amax1(0.0, elij(il,i,j)) | |
218 | 12499941 | nent(il, i) = nent(il, i) + 1 | |
219 | END IF | ||
220 | |||
221 | 14357294 | Sij(il, i, j) = amax1(0.0, Sij(il,i,j)) | |
222 | 14357294 | Sij(il, i, j) = amin1(1.0, Sij(il,i,j)) | |
223 |
2/2✓ Branch 0 taken 35272582 times.
✓ Branch 1 taken 39592218 times.
|
74864800 | ELSE IF (j > i) THEN |
224 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 35272582 times.
|
35272582 | IF (prt_level >= 10) THEN |
225 | ✗ | print *,'cv3p_mixing i, j, Sij given by the no-precip eq. ', i, j, Sij(il,i,j) | |
226 | ENDIF | ||
227 | END IF ! new | ||
228 | END DO | ||
229 | END DO | ||
230 | ELSE ! (ok_entrain) | ||
231 | ✗ | DO il = 1,ncum | |
232 | ✗ | nent(il,i) = 0 | |
233 | ENDDO | ||
234 | ENDIF ! (ok_entrain) | ||
235 | |||
236 | !jygdebug< | ||
237 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 6240 times.
|
6240 | IF (prt_level >= 10) THEN |
238 | ✗ | print *,'cv3p_mixing i, nent(i), icb, inb ',i, nent(igout,i), icb(igout), inb(igout) | |
239 | ✗ | IF (nent(igout,i) .gt. 0) THEN | |
240 | ✗ | print *,'i,(j,Sij(i,j),j=icb-1,inb) ',i,(j,Sij(igout,i,j),j=icb(igout)-1,inb(igout)) | |
241 | ENDIF | ||
242 | ENDIF | ||
243 | !>jygdebug | ||
244 | |||
245 | ! *** if no air can entrain at level i assume that updraft detrains *** | ||
246 | ! *** at that level and calculate detrained air flux and properties *** | ||
247 | |||
248 | |||
249 | ! @ do 170 i=icb(il),inb(il) | ||
250 | |||
251 |
2/2✓ Branch 0 taken 3304522 times.
✓ Branch 1 taken 6240 times.
|
3311002 | DO il = 1, ncum |
252 |
5/6✓ Branch 0 taken 2911168 times.
✓ Branch 1 taken 393354 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1766464 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1144704 times.
|
3310762 | IF ((i>=icb(il)) .AND. (i<=inb(il)) .AND. (nent(il,i)==0)) THEN |
253 | ! @ if(nent(il,i).eq.0)then | ||
254 | !!! Ment(il,i,i)=m(il,i) | ||
255 | ✗ | Ment(il, i, i) = 1. | |
256 | !! Qent(il, i, i) = qnk(il) - ep(il, i)*clw(il, i) | ||
257 | ✗ | Qent(il, i, i) = qta(il,i-1) - ep(il, i)*clw(il, i) | |
258 | ✗ | uent(il, i, i) = unk(il) | |
259 | ✗ | vent(il, i, i) = vnk(il) | |
260 | ✗ | IF (fl_cor_ebil .GE. 2) THEN | |
261 | ✗ | hent(il, i, i) = hp(il,i) | |
262 | ENDIF | ||
263 | ✗ | elij(il, i, i) = clw(il, i)*(1.-ep(il,i)) | |
264 | ✗ | Sij(il, i, i) = 0.0 | |
265 | END IF | ||
266 | END DO | ||
267 | END DO ! i = minorig + 1, nl | ||
268 | |||
269 | !jyg! DO j = 1, ntra | ||
270 | !jyg! DO i = minorig + 1, nl | ||
271 | !jyg! DO il = 1, ncum | ||
272 | !jyg! IF (i>=icb(il) .AND. i<=inb(il) .AND. nent(il,i)==0) THEN | ||
273 | !jyg! traent(il, i, i, j) = tra(il, nk(il), j) | ||
274 | !jyg! END IF | ||
275 | !jyg! END DO | ||
276 | !jyg! END DO | ||
277 | !jyg! END DO | ||
278 | |||
279 |
2/2✓ Branch 0 taken 6480 times.
✓ Branch 1 taken 240 times.
|
6720 | DO j = minorig, nl |
280 |
2/2✓ Branch 0 taken 174960 times.
✓ Branch 1 taken 6480 times.
|
181680 | DO i = minorig, nl |
281 |
2/2✓ Branch 0 taken 92653713 times.
✓ Branch 1 taken 174960 times.
|
92835153 | DO il = 1, ncum |
282 | IF ((j>=(icb(il)-1)) .AND. (j<=inb(il)) .AND. & | ||
283 |
8/8✓ Branch 0 taken 82033155 times.
✓ Branch 1 taken 10620558 times.
✓ Branch 2 taken 34338627 times.
✓ Branch 3 taken 47694528 times.
✓ Branch 4 taken 29340937 times.
✓ Branch 5 taken 4997690 times.
✓ Branch 6 taken 14357294 times.
✓ Branch 7 taken 14983643 times.
|
92828673 | (i>=icb(il)) .AND. (i<=inb(il))) THEN |
284 | 14357294 | Sigij(il, i, j) = Sij(il, i, j) | |
285 | END IF | ||
286 | END DO | ||
287 | END DO | ||
288 | END DO | ||
289 | ! @ enddo | ||
290 | |||
291 | ! @170 continue | ||
292 | |||
293 | ! ===================================================================== | ||
294 | ! --- NORMALIZE ENTRAINED AIR MASS FLUXES | ||
295 | ! --- TO REPRESENT EQUAL PROBABILITIES OF MIXING | ||
296 | ! ===================================================================== | ||
297 | |||
298 | 240 | CALL zilch(csum, nloc*nd) | |
299 | |||
300 |
2/2✓ Branch 0 taken 127097 times.
✓ Branch 1 taken 240 times.
|
127337 | DO il = 1, ncum |
301 | 127337 | lwork(il) = .FALSE. | |
302 | END DO | ||
303 | |||
304 | ! --------------------------------------------------------------- | ||
305 |
2/2✓ Branch 0 taken 6240 times.
✓ Branch 1 taken 240 times.
|
6480 | DO i = minorig + 1, nl !Loop on origin level "i" |
306 | ! --------------------------------------------------------------- | ||
307 | |||
308 | num1 = 0 | ||
309 |
2/2✓ Branch 0 taken 3304522 times.
✓ Branch 1 taken 6240 times.
|
3310762 | DO il = 1, ncum |
310 |
4/4✓ Branch 0 taken 2911168 times.
✓ Branch 1 taken 393354 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1766464 times.
|
3310762 | IF (i>=icb(il) .AND. i<=inb(il)) num1 = num1 + 1 |
311 | END DO | ||
312 |
2/2✓ Branch 0 taken 4633 times.
✓ Branch 1 taken 1607 times.
|
6240 | IF (num1<=0) GO TO 789 |
313 | |||
314 | |||
315 | !JYG1 Find maximum of SIJ for J>I, if any. | ||
316 | |||
317 |
2/2✓ Branch 0 taken 4605202 times.
✓ Branch 1 taken 4633 times.
|
4609835 | Sx(:) = 0. |
318 | |||
319 |
2/2✓ Branch 0 taken 2453448 times.
✓ Branch 1 taken 4633 times.
|
2458081 | DO il = 1, ncum |
320 |
4/4✓ Branch 0 taken 2187191 times.
✓ Branch 1 taken 266257 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1042487 times.
|
2458081 | IF (i>=icb(il) .AND. i<=inb(il)) THEN |
321 | 1144704 | signhpmh(il) = sign(1., hp(il,i)-h(il,i)) | |
322 | 1144704 | Sbef(il) = max(0., signhpmh(il)) | |
323 | END IF | ||
324 | END DO | ||
325 | |||
326 |
2/2✓ Branch 0 taken 68765 times.
✓ Branch 1 taken 4633 times.
|
73398 | DO j = i + 1, nl |
327 |
2/2✓ Branch 0 taken 36415670 times.
✓ Branch 1 taken 68765 times.
|
36489068 | DO il = 1, ncum |
328 |
6/6✓ Branch 0 taken 30299975 times.
✓ Branch 1 taken 6115695 times.
✓ Branch 2 taken 19251122 times.
✓ Branch 3 taken 11048853 times.
✓ Branch 4 taken 6033943 times.
✓ Branch 5 taken 13217179 times.
|
36484435 | IF (i>=icb(il) .AND. i<=inb(il) .AND. j<=inb(il)) THEN |
329 |
2/2✓ Branch 0 taken 984135 times.
✓ Branch 1 taken 5049808 times.
|
6033943 | IF (Sbef(il)<Sij(il,i,j)) THEN |
330 | 984135 | Sx(il) = max(Sij(il,i,j), Sx(il)) | |
331 | END IF | ||
332 | 6033943 | Sbef(il) = Sij(il, i, j) | |
333 | END IF | ||
334 | END DO | ||
335 | END DO | ||
336 | |||
337 | |||
338 |
2/2✓ Branch 0 taken 2453448 times.
✓ Branch 1 taken 4633 times.
|
2458081 | DO il = 1, ncum |
339 |
4/4✓ Branch 0 taken 2187191 times.
✓ Branch 1 taken 266257 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1042487 times.
|
2458081 | IF (i>=icb(il) .AND. i<=inb(il)) THEN |
340 | 1144704 | lwork(il) = (nent(il,i)/=0) | |
341 | !! rti = qnk(il) - ep(il, i)*clw(il, i) | ||
342 | 1144704 | rti = qta(il,i-1) - ep(il, i)*clw(il, i) | |
343 | !jyg< | ||
344 |
1/2✓ Branch 0 taken 1144704 times.
✗ Branch 1 not taken.
|
1144704 | IF (cvflag_ice) THEN |
345 | |||
346 | anum = h(il, i) - hp(il, i) - (lv(il,i)+frac(il,i)*lf(il,i))* & | ||
347 | 1144704 | (rti-rs(il,i)) + (cpv-cpd)*t(il, i)*(rti-rr(il,i)) | |
348 | denom = h(il, i) - hp(il, i) + (lv(il,i)+frac(il,i)*lf(il,i))* & | ||
349 | 1144704 | (rr(il,i)-rti) + (cpd-cpv)*t(il, i)*(rr(il,i)-rti) | |
350 | ELSE | ||
351 | |||
352 | anum = h(il, i) - hp(il, i) - lv(il, i)*(rti-rs(il,i)) + & | ||
353 | ✗ | (cpv-cpd)*t(il, i)*(rti-rr(il,i)) | |
354 | denom = h(il, i) - hp(il, i) + lv(il, i)*(rr(il,i)-rti) + & | ||
355 | ✗ | (cpd-cpv)*t(il, i)*(rr(il,i)-rti) | |
356 | END IF | ||
357 | !>jyg | ||
358 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1144704 times.
|
1144704 | IF (abs(denom)<0.01) denom = 0.01 |
359 | 1144704 | Scrit(il) = min(anum/denom, 1.) | |
360 | 1144704 | alt = rti - rs(il, i) + Scrit(il)*(rr(il,i)-rti) | |
361 | |||
362 | !JYG1 Find new critical value Scrit2 | ||
363 | ! such that : Sij > Scrit2 => mixed draught will detrain at J<I | ||
364 | ! Sij < Scrit2 => mixed draught will detrain at J>I | ||
365 | |||
366 | Scrit2 = min(Scrit(il), Sx(il))*max(0., -signhpmh(il)) + & | ||
367 | 1144704 | Scrit(il)*max(0., signhpmh(il)) | |
368 | |||
369 | 1144704 | Scrit(il) = Scrit2 | |
370 | |||
371 | !JYG Correction pour la nouvelle logique; la correction pour ALT | ||
372 | ! est un peu au hazard | ||
373 |
2/2✓ Branch 0 taken 193829 times.
✓ Branch 1 taken 950875 times.
|
1144704 | IF (Scrit(il)<=0.0) Scrit(il) = 0.0 |
374 |
2/2✓ Branch 0 taken 750019 times.
✓ Branch 1 taken 394685 times.
|
1144704 | IF (alt<=0.0) Scrit(il) = 1.0 |
375 | |||
376 | 1144704 | smax(il) = 0.0 | |
377 | 1144704 | ASij(il) = 0.0 | |
378 | 1144704 | sup(il) = 0. ! upper S-value reached by descending draughts | |
379 | END IF | ||
380 | END DO | ||
381 | |||
382 | ! --------------------------------------------------------------- | ||
383 |
2/2✓ Branch 0 taken 125091 times.
✓ Branch 1 taken 4633 times.
|
129724 | DO j = minorig, nl !Loop on destination level "j" |
384 | ! --------------------------------------------------------------- | ||
385 | |||
386 | num2 = 0 | ||
387 |
2/2✓ Branch 0 taken 66243096 times.
✓ Branch 1 taken 125091 times.
|
66368187 | DO il = 1, ncum |
388 | IF (i>=icb(il) .AND. i<=inb(il) .AND. & | ||
389 |
9/10✓ Branch 0 taken 59054157 times.
✓ Branch 1 taken 7188939 times.
✓ Branch 2 taken 30907008 times.
✓ Branch 3 taken 28147149 times.
✓ Branch 4 taken 27574473 times.
✓ Branch 5 taken 3332535 times.
✓ Branch 6 taken 14357294 times.
✓ Branch 7 taken 13217179 times.
✓ Branch 8 taken 14357294 times.
✗ Branch 9 not taken.
|
66243096 | j>=(icb(il)-1) .AND. j<=inb(il) .AND. & |
390 | 14482385 | lwork(il)) num2 = num2 + 1 | |
391 | END DO | ||
392 |
2/2✓ Branch 0 taken 94002 times.
✓ Branch 1 taken 31089 times.
|
125091 | IF (num2<=0) GO TO 175 |
393 | |||
394 | ! ----------------------------------------------- | ||
395 |
2/2✓ Branch 0 taken 42368 times.
✓ Branch 1 taken 51634 times.
|
94002 | IF (j>i) THEN |
396 | ! ----------------------------------------------- | ||
397 |
2/2✓ Branch 0 taken 42368 times.
✓ Branch 1 taken 22436073 times.
|
22478441 | DO il = 1, ncum |
398 | IF (i>=icb(il) .AND. i<=inb(il) .AND. & | ||
399 |
8/10✓ Branch 0 taken 17864651 times.
✓ Branch 1 taken 4571422 times.
✓ Branch 2 taken 12722753 times.
✓ Branch 3 taken 5141898 times.
✓ Branch 4 taken 12722753 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 6033943 times.
✓ Branch 7 taken 6688810 times.
✓ Branch 8 taken 6033943 times.
✗ Branch 9 not taken.
|
22436073 | j>=(icb(il)-1) .AND. j<=inb(il) .AND. & |
400 | 42368 | lwork(il)) THEN | |
401 |
2/2✓ Branch 0 taken 5153725 times.
✓ Branch 1 taken 880218 times.
|
6033943 | IF (Sij(il,i,j)>0.0) THEN |
402 | 5153725 | Smid(il) = min(Sij(il,i,j), Scrit(il)) | |
403 | 5153725 | Sjmax(il) = Smid(il) | |
404 | 5153725 | Sjmin(il) = Smid(il) | |
405 |
4/4✓ Branch 0 taken 4633253 times.
✓ Branch 1 taken 520472 times.
✓ Branch 2 taken 2400980 times.
✓ Branch 3 taken 2232273 times.
|
5153725 | IF (Smid(il)<smin(il) .AND. Sij(il,i,j+1)<Smid(il)) THEN |
406 | 2400980 | smin(il) = Smid(il) | |
407 | 2400980 | Sjmax(il) = min((Sij(il,i,j+1)+Sij(il,i,j))/2., Sij(il,i,j), Scrit(il)) | |
408 | 2400980 | Sjmin(il) = max((Sbef(il)+Sij(il,i,j))/2., Sij(il,i,j)) | |
409 | 2400980 | Sjmin(il) = min(Sjmin(il), Scrit(il)) | |
410 | 2400980 | Sbef(il) = Sij(il, i, j) | |
411 | END IF | ||
412 | END IF | ||
413 | END IF | ||
414 | END DO | ||
415 | ! ----------------------------------------------- | ||
416 |
2/2✓ Branch 0 taken 4633 times.
✓ Branch 1 taken 47001 times.
|
51634 | ELSE IF (j==i) THEN |
417 | ! ----------------------------------------------- | ||
418 |
2/2✓ Branch 0 taken 4633 times.
✓ Branch 1 taken 2453448 times.
|
2458081 | DO il = 1, ncum |
419 | IF (i>=icb(il) .AND. i<=inb(il) .AND. & | ||
420 |
7/10✓ Branch 0 taken 2187191 times.
✓ Branch 1 taken 266257 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1042487 times.
✓ Branch 4 taken 1144704 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1144704 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 1144704 times.
✗ Branch 9 not taken.
|
2453448 | j>=(icb(il)-1) .AND. j<=inb(il) .AND. & |
421 | 4633 | lwork(il)) THEN | |
422 |
1/2✓ Branch 0 taken 1144704 times.
✗ Branch 1 not taken.
|
1144704 | IF (Sij(il,i,j)>0.0) THEN |
423 | 1144704 | Smid(il) = 1. | |
424 | Sjmin(il) = max((Sij(il,i,j-1)+Smid(il))/2., Scrit(il))*max(0., -signhpmh(il)) + & | ||
425 | 1144704 | min((Sij(il,i,j+1)+Smid(il))/2., Scrit(il))*max(0., signhpmh(il)) | |
426 | 1144704 | Sjmin(il) = max(Sjmin(il), sup(il)) | |
427 | 1144704 | Sjmax(il) = 1. | |
428 | |||
429 | ! - preparation des variables Scrit, Smin et Sbef pour la partie j>i | ||
430 | 1144704 | Scrit(il) = min(Sjmin(il), Sjmax(il), Scrit(il)) | |
431 | |||
432 | 1144704 | smin(il) = 1. | |
433 | 1144704 | Sbef(il) = max(0., signhpmh(il)) | |
434 | 1144704 | supmax(il, i) = sign(Scrit(il), -signhpmh(il)) | |
435 | END IF | ||
436 | END IF | ||
437 | END DO | ||
438 | ! ----------------------------------------------- | ||
439 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 47001 times.
|
47001 | ELSE IF (j<i) THEN |
440 | ! ----------------------------------------------- | ||
441 |
2/2✓ Branch 0 taken 47001 times.
✓ Branch 1 taken 24889521 times.
|
24936522 | DO il = 1, ncum |
442 | IF (i>=icb(il) .AND. i<=inb(il) .AND. & | ||
443 |
8/10✓ Branch 0 taken 24348791 times.
✓ Branch 1 taken 540730 times.
✓ Branch 2 taken 9364655 times.
✓ Branch 3 taken 14984136 times.
✓ Branch 4 taken 7178647 times.
✓ Branch 5 taken 2186008 times.
✓ Branch 6 taken 7178647 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 7178647 times.
✗ Branch 9 not taken.
|
24889521 | j>=(icb(il)-1) .AND. j<=inb(il) .AND. & |
444 | 47001 | lwork(il)) THEN | |
445 |
2/2✓ Branch 0 taken 6201512 times.
✓ Branch 1 taken 977135 times.
|
7178647 | IF (Sij(il,i,j)>0.0) THEN |
446 | 6201512 | Smid(il) = max(Sij(il,i,j), Scrit(il)) | |
447 | 6201512 | Sjmax(il) = Smid(il) | |
448 | 6201512 | Sjmin(il) = Smid(il) | |
449 |
4/4✓ Branch 0 taken 6201509 times.
✓ Branch 1 taken 3 times.
✓ Branch 2 taken 58118 times.
✓ Branch 3 taken 6143391 times.
|
6201512 | IF (Smid(il)>smax(il) .AND. Sij(il,i,j+1)>Smid(il)) THEN |
450 | 58118 | smax(il) = Smid(il) | |
451 | 58118 | Sjmax(il) = max((Sij(il,i,j+1)+Sij(il,i,j))/2., Sij(il,i,j)) | |
452 | 58118 | Sjmax(il) = max(Sjmax(il), Scrit(il)) | |
453 | 58118 | Sjmin(il) = min((Sbef(il)+Sij(il,i,j))/2., Sij(il,i,j)) | |
454 | 58118 | Sjmin(il) = max(Sjmin(il), Scrit(il)) | |
455 | 58118 | Sbef(il) = Sij(il, i, j) | |
456 | END IF | ||
457 |
2/2✓ Branch 0 taken 47816 times.
✓ Branch 1 taken 6153696 times.
|
6201512 | IF (abs(Sjmin(il)-Sjmax(il))>1.E-10) & |
458 | 47816 | sup(il) = max(Sjmin(il), Sjmax(il), sup(il)) | |
459 | END IF | ||
460 | END IF | ||
461 | END DO | ||
462 | ! ----------------------------------------------- | ||
463 | END IF | ||
464 | ! ----------------------------------------------- | ||
465 | |||
466 | |||
467 |
2/2✓ Branch 0 taken 49779042 times.
✓ Branch 1 taken 94002 times.
|
49873044 | DO il = 1, ncum |
468 | IF (i>=icb(il) .AND. i<=inb(il) .AND. & | ||
469 |
9/10✓ Branch 0 taken 44400633 times.
✓ Branch 1 taken 5378409 times.
✓ Branch 2 taken 23232112 times.
✓ Branch 3 taken 21168521 times.
✓ Branch 4 taken 21046104 times.
✓ Branch 5 taken 2186008 times.
✓ Branch 6 taken 14357294 times.
✓ Branch 7 taken 6688810 times.
✓ Branch 8 taken 14357294 times.
✗ Branch 9 not taken.
|
49779042 | j>=(icb(il)-1) .AND. j<=inb(il) .AND. & |
470 | 94002 | lwork(il)) THEN | |
471 |
2/2✓ Branch 0 taken 12499941 times.
✓ Branch 1 taken 1857353 times.
|
14357294 | IF (Sij(il,i,j)>0.0) THEN |
472 | !! rti = qnk(il) - ep(il, i)*clw(il, i) | ||
473 | 12499941 | rti = qta(il,i-1) - ep(il, i)*clw(il, i) | |
474 | 12499941 | Qmixmax(il) = Qmix(Sjmax(il)) | |
475 | 12499941 | Qmixmin(il) = Qmix(Sjmin(il)) | |
476 | 12499941 | Rmixmax(il) = Rmix(Sjmax(il)) | |
477 | 12499941 | Rmixmin(il) = Rmix(Sjmin(il)) | |
478 | 12499941 | sqmrmax(il) = Sjmax(il)*Qmix(Sjmax(il)) - Rmix(Sjmax(il)) | |
479 | 12499941 | sqmrmin(il) = Sjmin(il)*Qmix(Sjmin(il)) - Rmix(Sjmin(il)) | |
480 | |||
481 | 12499941 | Ment(il, i, j) = abs(Qmixmax(il)-Qmixmin(il))*Ment(il, i, j) | |
482 | |||
483 | ! Sigij(i,j) is the 'true' mixing fraction of mixture Ment(i,j) | ||
484 |
2/2✓ Branch 0 taken 3358784 times.
✓ Branch 1 taken 9141157 times.
|
12499941 | IF (abs(Qmixmax(il)-Qmixmin(il))>1.E-10) THEN |
485 | 3358784 | Sigij(il, i, j) = (sqmrmax(il)-sqmrmin(il))/(Qmixmax(il)-Qmixmin(il)) | |
486 | ELSE | ||
487 | 9141157 | Sigij(il, i, j) = 0. | |
488 | END IF | ||
489 | |||
490 | ! -- Compute Qent, uent, vent according to the true mixing fraction | ||
491 | 12499941 | Qent(il, i, j) = (1.-Sigij(il,i,j))*rti + Sigij(il, i, j)*rr(il, i) | |
492 | 12499941 | uent(il, i, j) = (1.-Sigij(il,i,j))*unk(il) + Sigij(il, i, j)*u(il, i) | |
493 | 12499941 | vent(il, i, j) = (1.-Sigij(il,i,j))*vnk(il) + Sigij(il, i, j)*v(il, i) | |
494 | |||
495 | ! -- Compute liquid water static energy of mixed draughts | ||
496 | ! IF (j .GT. i) THEN | ||
497 | ! awat=elij(il,i,j)-(1.-ep(il,j))*clw(il,j) | ||
498 | ! awat=amax1(awat,0.0) | ||
499 | ! ELSE | ||
500 | ! awat = 0. | ||
501 | ! ENDIF | ||
502 | ! Hent(il,i,j) = (1.-Sigij(il,i,j))*HP(il,i) | ||
503 | ! : + Sigij(il,i,j)*H(il,i) | ||
504 | ! : + (LV(il,j)+(cpd-cpv)*t(il,j))*awat | ||
505 | !IM 301008 beg | ||
506 | 12499941 | hent(il, i, j) = (1.-Sigij(il,i,j))*hp(il, i) + Sigij(il, i, j)*h(il, i) | |
507 | |||
508 | !jyg< | ||
509 | ! elij(il, i, j) = Qent(il, i, j) - rs(il, j) | ||
510 | ! elij(il, i, j) = elij(il, i, j) + & | ||
511 | ! ((h(il,j)-hent(il,i,j))*rs(il,j)*lv(il,j) / & | ||
512 | ! ((cpd*(1.-Qent(il,i,j))+Qent(il,i,j)*cpv)*rrv*t(il,j)*t(il,j))) | ||
513 | ! elij(il, i, j) = elij(il, i, j) / & | ||
514 | ! (1.+lv(il,j)*lv(il,j)*rs(il,j) / & | ||
515 | ! ((cpd*(1.-Qent(il,i,j))+Qent(il,i,j)*cpv)*rrv*t(il,j)*t(il,j))) | ||
516 | ! | ||
517 | ! Computation of condensate amount Elij, taking into account the ice fraction frac | ||
518 | ! Warning : the same saturation humidity rs is used over both liquid water and ice; this | ||
519 | ! should be corrected. | ||
520 | ! | ||
521 | ! Heat capacity of mixed draught | ||
522 | 12499941 | cpm = cpd+Qent(il,i,j)*(cpv-cpd) | |
523 | ! | ||
524 |
3/4✓ Branch 0 taken 12499941 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5006923 times.
✓ Branch 3 taken 7493018 times.
|
12499941 | IF (cvflag_ice .and. frac(il,j) .gt. 0.) THEN |
525 | 5006923 | elij(il, i, j) = Qent(il, i, j) - rs(il, j) | |
526 | elij(il, i, j) = elij(il, i, j) + & | ||
527 | (h(il,j)-hent(il,i,j)+(cpv-cpd)*(Qent(il,i,j)-rr(il,j))*t(il,j))* & | ||
528 | 5006923 | rs(il,j)*lv(il,j) / (cpm*rrv*t(il,j)*t(il,j)) | |
529 | elij(il, i, j) = elij(il, i, j) / & | ||
530 | (1.+(lv(il,j)+frac(il,j)*lf(il,j))*lv(il,j)*rs(il,j) / & | ||
531 | 5006923 | (cpm*rrv*t(il,j)*t(il,j))) | |
532 | ELSE | ||
533 | 7493018 | elij(il, i, j) = Qent(il, i, j) - rs(il, j) | |
534 | elij(il, i, j) = elij(il, i, j) + & | ||
535 | (h(il,j)-hent(il,i,j)+(cpv-cpd)*(Qent(il,i,j)-rr(il,j))*t(il,j))* & | ||
536 | 7493018 | rs(il,j)*lv(il,j) / (cpm*rrv*t(il,j)*t(il,j)) | |
537 | elij(il, i, j) = elij(il, i, j) / & | ||
538 | (1.+lv(il,j)*lv(il,j)*rs(il,j) / & | ||
539 | 7493018 | (cpm*rrv*t(il,j)*t(il,j))) | |
540 | ENDIF | ||
541 | !>jyg | ||
542 | 12499941 | elij(il, i, j) = max(elij(il,i,j), 0.) | |
543 | |||
544 | 12499941 | elij(il, i, j) = min(elij(il,i,j), Qent(il,i,j)) | |
545 | |||
546 |
2/2✓ Branch 0 taken 5153725 times.
✓ Branch 1 taken 7346216 times.
|
12499941 | IF (j>i) THEN |
547 | 5153725 | awat = elij(il, i, j) - (1.-ep(il,j))*clw(il, j) | |
548 | 5153725 | awat = amax1(awat, 0.0) | |
549 | ELSE | ||
550 | awat = 0. | ||
551 | END IF | ||
552 | |||
553 | ! print *,h(il,j)-hent(il,i,j),LV(il,j)*rs(il,j)/(cpd*rrv*t(il,j)* | ||
554 | ! : t(il,j)) | ||
555 | |||
556 | !jyg< | ||
557 | ! hent(il, i, j) = hent(il, i, j) + (lv(il,j)+(cpd-cpv)*t(il,j))*awat | ||
558 | ! Mixed draught temperature at level j | ||
559 |
3/4✓ Branch 0 taken 12499941 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 5006923 times.
✓ Branch 3 taken 7493018 times.
|
12499941 | IF (cvflag_ice .and. frac(il,j) .gt. 0.) THEN |
560 | 5006923 | Tm = t(il,j) + (Qent(il,i,j)-elij(il,i,j)-rs(il,j))*rrv*t(il,j)*t(il,j)/(lv(il,j)*rs(il,j)) | |
561 | 5006923 | hent(il, i, j) = hent(il, i, j) + (lv(il,j)+frac(il,j)*lf(il,j)+(cpd-cpv)*Tm)*awat | |
562 | ELSE | ||
563 | 7493018 | Tm = t(il,j) + (Qent(il,i,j)-elij(il,i,j)-rs(il,j))*rrv*t(il,j)*t(il,j)/(lv(il,j)*rs(il,j)) | |
564 | 7493018 | hent(il, i, j) = hent(il, i, j) + (lv(il,j)+(cpd-cpv)*Tm)*awat | |
565 | ENDIF | ||
566 | !>jyg | ||
567 | |||
568 | !IM 301008 end | ||
569 | |||
570 | ! print *,'mix : i,j,hent(il,i,j),Sigij(il,i,j) ', | ||
571 | ! : i,j,hent(il,i,j),Sigij(il,i,j) | ||
572 | |||
573 | ! -- ASij is the integral of P(F) over the relevant F interval | ||
574 | ASij(il) = ASij(il) + abs(Qmixmax(il)*(1.-Sjmax(il))+Rmixmax(il) - & | ||
575 | 12499941 | Qmixmin(il)*(1.-Sjmin(il))-Rmixmin(il)) | |
576 | |||
577 | END IF | ||
578 | END IF | ||
579 | END DO | ||
580 | !jyg! DO k = 1, ntra | ||
581 | !jyg! DO il = 1, ncum | ||
582 | !jyg! IF ((i>=icb(il)) .AND. (i<=inb(il)) .AND. & | ||
583 | !jyg! (j>=(icb(il)-1)) .AND. (j<=inb(il)) .AND. & | ||
584 | !jyg! lwork(il)) THEN | ||
585 | !jyg! IF (Sij(il,i,j)>0.0) THEN | ||
586 | !jyg! traent(il, i, j, k) = Sigij(il, i, j)*tra(il, i, k) + & | ||
587 | !jyg! (1.-Sigij(il,i,j))*tra(il, nk(il), k) | ||
588 | !jyg! END IF | ||
589 | !jyg! END IF | ||
590 | !jyg! END DO | ||
591 | !jyg! END DO | ||
592 | |||
593 | ! -- If I=J (detrainement and entrainement at the same level), then only the | ||
594 | ! -- adiabatic ascent part of the mixture is considered | ||
595 |
2/2✓ Branch 0 taken 4633 times.
✓ Branch 1 taken 89369 times.
|
94002 | IF (i==j) THEN |
596 |
2/2✓ Branch 0 taken 2453448 times.
✓ Branch 1 taken 4633 times.
|
2458081 | DO il = 1, ncum |
597 | IF (i>=icb(il) .AND. i<=inb(il) .AND. & | ||
598 |
7/10✓ Branch 0 taken 2187191 times.
✓ Branch 1 taken 266257 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1042487 times.
✓ Branch 4 taken 1144704 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1144704 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 1144704 times.
✗ Branch 9 not taken.
|
2453448 | j>=(icb(il)-1) .AND. j<=inb(il) .AND. & |
599 | 4633 | lwork(il)) THEN | |
600 |
1/2✓ Branch 0 taken 1144704 times.
✗ Branch 1 not taken.
|
1144704 | IF (Sij(il,i,j)>0.0) THEN |
601 | !! rti = qnk(il) - ep(il, i)*clw(il, i) | ||
602 | 1144704 | rti = qta(il,i-1) - ep(il, i)*clw(il, i) | |
603 | !!! Ment(il,i,i) = m(il,i)*abs(Qmixmax(il)*(1.-Sjmax(il)) | ||
604 | Ment(il, i, i) = abs(Qmixmax(il)*(1.-Sjmax(il))+Rmixmax(il) - & | ||
605 | 1144704 | Qmixmin(il)*(1.-Sjmin(il))-Rmixmin(il)) | |
606 | 1144704 | Qent(il, i, i) = rti | |
607 | 1144704 | uent(il, i, i) = unk(il) | |
608 | 1144704 | vent(il, i, i) = vnk(il) | |
609 | 1144704 | hent(il, i, i) = hp(il, i) | |
610 | 1144704 | elij(il, i, i) = clw(il, i)*(1.-ep(il,i)) | |
611 | 1144704 | Sigij(il, i, i) = 0. | |
612 | END IF | ||
613 | END IF | ||
614 | END DO | ||
615 | !jyg! DO k = 1, ntra | ||
616 | !jyg! DO il = 1, ncum | ||
617 | !jyg! IF ((i>=icb(il)) .AND. (i<=inb(il)) .AND. & | ||
618 | !jyg! (j>=(icb(il)-1)) .AND. (j<=inb(il)) .AND. & | ||
619 | !jyg! lwork(il)) THEN | ||
620 | !jyg! IF (Sij(il,i,j)>0.0) THEN | ||
621 | !jyg! traent(il, i, i, k) = tra(il, nk(il), k) | ||
622 | !jyg! END IF | ||
623 | !jyg! END IF | ||
624 | !jyg! END DO | ||
625 | !jyg! END DO | ||
626 | |||
627 | END IF | ||
628 | |||
629 | ! --------------------------------------------------------------- | ||
630 | 4633 | 175 END DO ! End loop on destination level "j" | |
631 | ! --------------------------------------------------------------- | ||
632 | |||
633 |
2/2✓ Branch 0 taken 2453448 times.
✓ Branch 1 taken 4633 times.
|
2458081 | DO il = 1, ncum |
634 |
5/6✓ Branch 0 taken 2187191 times.
✓ Branch 1 taken 266257 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1042487 times.
✓ Branch 4 taken 1144704 times.
✗ Branch 5 not taken.
|
2458081 | IF (i>=icb(il) .AND. i<=inb(il) .AND. lwork(il)) THEN |
635 | 1144704 | ASij(il) = amax1(1.0E-16, ASij(il)) | |
636 | !jyg+lluis< | ||
637 | !! ASij(il) = 1.0/ASij(il) | ||
638 | 1144704 | ASij_inv(il) = 1.0/ASij(il) | |
639 | ! IF the F-interval spanned by possible mixtures is less than 0.01, no mixing occurs | ||
640 |
2/2✓ Branch 0 taken 70188 times.
✓ Branch 1 taken 1074516 times.
|
1144704 | IF (ASij_inv(il) > 100.) ASij_inv(il) = 0. |
641 | !>jyg+lluis | ||
642 | 1144704 | csum(il, i) = 0.0 | |
643 | END IF | ||
644 | END DO | ||
645 | |||
646 |
2/2✓ Branch 0 taken 125091 times.
✓ Branch 1 taken 4633 times.
|
129724 | DO j = minorig, nl |
647 |
2/2✓ Branch 0 taken 66243096 times.
✓ Branch 1 taken 125091 times.
|
66372820 | DO il = 1, ncum |
648 | IF (i>=icb(il) .AND. i<=inb(il) .AND. lwork(il) .AND. & | ||
649 |
9/10✓ Branch 0 taken 59054157 times.
✓ Branch 1 taken 7188939 times.
✓ Branch 2 taken 30907008 times.
✓ Branch 3 taken 28147149 times.
✓ Branch 4 taken 30907008 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 27574473 times.
✓ Branch 7 taken 3332535 times.
✓ Branch 8 taken 14357294 times.
✓ Branch 9 taken 13217179 times.
|
66368187 | j>=(icb(il)-1) .AND. j<=inb(il)) THEN |
650 | !jyg Ment(il, i, j) = Ment(il, i, j)*ASij(il) | ||
651 | 14357294 | Ment(il, i, j) = Ment(il, i, j)*ASij_inv(il) | |
652 | END IF | ||
653 | END DO | ||
654 | END DO | ||
655 | |||
656 |
2/2✓ Branch 0 taken 125091 times.
✓ Branch 1 taken 4633 times.
|
129724 | DO j = minorig, nl |
657 |
2/2✓ Branch 0 taken 66243096 times.
✓ Branch 1 taken 125091 times.
|
66372820 | DO il = 1, ncum |
658 | IF (i>=icb(il) .AND. i<=inb(il) .AND. lwork(il) .AND. & | ||
659 |
9/10✓ Branch 0 taken 59054157 times.
✓ Branch 1 taken 7188939 times.
✓ Branch 2 taken 30907008 times.
✓ Branch 3 taken 28147149 times.
✓ Branch 4 taken 30907008 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 27574473 times.
✓ Branch 7 taken 3332535 times.
✓ Branch 8 taken 14357294 times.
✓ Branch 9 taken 13217179 times.
|
66368187 | j>=(icb(il)-1) .AND. j<=inb(il)) THEN |
660 | 14357294 | csum(il, i) = csum(il, i) + Ment(il, i, j) | |
661 | END IF | ||
662 | END DO | ||
663 | END DO | ||
664 | |||
665 |
2/2✓ Branch 0 taken 4633 times.
✓ Branch 1 taken 2453448 times.
|
2458081 | DO il = 1, ncum |
666 |
7/8✓ Branch 0 taken 2187191 times.
✓ Branch 1 taken 266257 times.
✓ Branch 2 taken 1144704 times.
✓ Branch 3 taken 1042487 times.
✓ Branch 4 taken 1144704 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 73348 times.
✓ Branch 7 taken 1071356 times.
|
2459688 | IF (i>=icb(il) .AND. i<=inb(il) .AND. lwork(il) .AND. csum(il,i)<1.) THEN |
667 | ! cc : .and. csum(il,i).lt.m(il,i) ) then | ||
668 | 73348 | nent(il, i) = 0 | |
669 | ! cc Ment(il,i,i)=m(il,i) | ||
670 | 73348 | Ment(il, i, i) = 1. | |
671 | !! Qent(il, i, i) = qnk(il) - ep(il, i)*clw(il, i) | ||
672 | 73348 | Qent(il, i, i) = qta(il,i-1) - ep(il, i)*clw(il, i) | |
673 | 73348 | uent(il, i, i) = unk(il) | |
674 | 73348 | vent(il, i, i) = vnk(il) | |
675 | 73348 | elij(il, i, i) = clw(il, i)*(1.-ep(il,i)) | |
676 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 73348 times.
|
73348 | IF (fl_cor_ebil .GE. 2) THEN |
677 | ✗ | hent(il, i, i) = hp(il,i) | |
678 | ✗ | Sigij(il, i, i) = 0.0 | |
679 | ELSE | ||
680 | 73348 | Sij(il, i, i) = 0.0 | |
681 | ENDIF | ||
682 | END IF | ||
683 | END DO ! il | ||
684 | |||
685 | !jyg! DO j = 1, ntra | ||
686 | !jyg! DO il = 1, ncum | ||
687 | !jyg! IF (i>=icb(il) .AND. i<=inb(il) .AND. lwork(il) .AND. csum(il,i)<1.) THEN | ||
688 | !jyg!! cc : .and. csum(il,i).lt.m(il,i) ) then | ||
689 | !jyg! traent(il, i, i, j) = tra(il, nk(il), j) | ||
690 | !jyg! END IF | ||
691 | !jyg! END DO | ||
692 | !jyg! END DO | ||
693 | |||
694 | ! --------------------------------------------------------------- | ||
695 | 240 | 789 END DO ! End loop on origin level "i" | |
696 | ! --------------------------------------------------------------- | ||
697 | |||
698 | |||
699 | 240 | RETURN | |
700 | END SUBROUTINE cv3p_mixing | ||
701 | |||
702 |