Directory: | ./ |
---|---|
File: | phys/climb_wind_mod.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 79 | 101 | 78.2% |
Branches: | 136 | 212 | 64.2% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | ! | ||
2 | MODULE climb_wind_mod | ||
3 | ! | ||
4 | ! Module to solve the verctical diffusion of the wind components "u" and "v". | ||
5 | ! | ||
6 | USE dimphy | ||
7 | |||
8 | IMPLICIT NONE | ||
9 | |||
10 | SAVE | ||
11 | PRIVATE | ||
12 | |||
13 | REAL, DIMENSION(:), ALLOCATABLE :: alf1, alf2 | ||
14 | !$OMP THREADPRIVATE(alf1,alf2) | ||
15 | REAL, DIMENSION(:,:), ALLOCATABLE :: Kcoefm | ||
16 | !$OMP THREADPRIVATE(Kcoefm) | ||
17 | REAL, DIMENSION(:,:), ALLOCATABLE :: Ccoef_U, Dcoef_U | ||
18 | !$OMP THREADPRIVATE(Ccoef_U, Dcoef_U) | ||
19 | REAL, DIMENSION(:,:), ALLOCATABLE :: Ccoef_V, Dcoef_V | ||
20 | !$OMP THREADPRIVATE(Ccoef_V, Dcoef_V) | ||
21 | REAL, DIMENSION(:), ALLOCATABLE :: Acoef_U, Bcoef_U | ||
22 | !$OMP THREADPRIVATE(Acoef_U, Bcoef_U) | ||
23 | REAL, DIMENSION(:), ALLOCATABLE :: Acoef_V, Bcoef_V | ||
24 | !$OMP THREADPRIVATE(Acoef_V, Bcoef_V) | ||
25 | LOGICAL :: firstcall=.TRUE. | ||
26 | !$OMP THREADPRIVATE(firstcall) | ||
27 | |||
28 | |||
29 | PUBLIC :: climb_wind_down, climb_wind_up | ||
30 | |||
31 | CONTAINS | ||
32 | ! | ||
33 | !**************************************************************************************** | ||
34 | ! | ||
35 | 29963097 | SUBROUTINE climb_wind_init | |
36 | |||
37 | INTEGER :: ierr | ||
38 | CHARACTER(len = 20) :: modname = 'climb_wind_init' | ||
39 | |||
40 | !**************************************************************************************** | ||
41 | ! Allocation of global module variables | ||
42 | ! | ||
43 | !**************************************************************************************** | ||
44 | |||
45 |
4/8✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
|
1 | ALLOCATE(alf1(klon), stat=ierr) |
46 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate alf1',1) |
47 | |||
48 |
4/8✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
|
1 | ALLOCATE(alf2(klon), stat=ierr) |
49 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate alf2',1) |
50 | |||
51 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Kcoefm(klon,klev), stat=ierr) |
52 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate Kcoefm',1) |
53 | |||
54 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Ccoef_U(klon,klev), stat=ierr) |
55 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocate Ccoef_U',1) |
56 | |||
57 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Dcoef_U(klon,klev), stat=ierr) |
58 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocation Dcoef_U',1) |
59 | |||
60 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Ccoef_V(klon,klev), stat=ierr) |
61 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocation Ccoef_V',1) |
62 | |||
63 |
6/12✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
|
2 | ALLOCATE(Dcoef_V(klon,klev), stat=ierr) |
64 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (ierr /= 0) CALL abort_physic(modname,'Pb in allocation Dcoef_V',1) |
65 | |||
66 |
16/32✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✗ Branch 7 not taken.
✓ Branch 8 taken 1 times.
✗ Branch 9 not taken.
✓ Branch 10 taken 1 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✓ Branch 13 taken 1 times.
✓ Branch 14 taken 1 times.
✗ Branch 15 not taken.
✓ Branch 16 taken 1 times.
✗ Branch 17 not taken.
✓ Branch 18 taken 1 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 1 times.
✓ Branch 22 taken 1 times.
✗ Branch 23 not taken.
✓ Branch 24 taken 1 times.
✗ Branch 25 not taken.
✓ Branch 26 taken 1 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 1 times.
✓ Branch 30 taken 1 times.
✗ Branch 31 not taken.
|
1 | ALLOCATE(Acoef_U(klon), Bcoef_U(klon), Acoef_V(klon), Bcoef_V(klon), STAT=ierr) |
67 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF ( ierr /= 0 ) PRINT*,' pb in allloc Acoef_U and Bcoef_U, ierr=', ierr |
68 | |||
69 | 1 | firstcall=.FALSE. | |
70 | |||
71 | 1 | END SUBROUTINE climb_wind_init | |
72 | ! | ||
73 | !**************************************************************************************** | ||
74 | ! | ||
75 | 1920 | SUBROUTINE climb_wind_down(knon, dtime, coef_in, pplay, paprs, temp, delp, u_old, v_old, & | |
76 | !!! nrlmd le 02/05/2011 | ||
77 | Ccoef_U_out, Ccoef_V_out, Dcoef_U_out, Dcoef_V_out, & | ||
78 | Kcoef_m_out, alf_1_out, alf_2_out, & | ||
79 | !!! | ||
80 | Acoef_U_out, Acoef_V_out, Bcoef_U_out, Bcoef_V_out) | ||
81 | ! | ||
82 | ! This routine calculates for the wind components u and v, | ||
83 | ! recursivly the coefficients C and D in equation | ||
84 | ! X(k) = C(k) + D(k)*X(k-1), X=[u,v], k=[1,klev] is the vertical layer. | ||
85 | ! | ||
86 | ! | ||
87 | |||
88 | ! Input arguments | ||
89 | !**************************************************************************************** | ||
90 | INTEGER, INTENT(IN) :: knon | ||
91 | REAL, INTENT(IN) :: dtime | ||
92 | REAL, DIMENSION(klon,klev), INTENT(IN) :: coef_in | ||
93 | REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay ! pres au milieu de couche (Pa) | ||
94 | REAL, DIMENSION(klon,klev+1), INTENT(IN) :: paprs ! pression a inter-couche (Pa) | ||
95 | REAL, DIMENSION(klon,klev), INTENT(IN) :: temp ! temperature | ||
96 | REAL, DIMENSION(klon,klev), INTENT(IN) :: delp | ||
97 | REAL, DIMENSION(klon,klev), INTENT(IN) :: u_old | ||
98 | REAL, DIMENSION(klon,klev), INTENT(IN) :: v_old | ||
99 | |||
100 | ! Output arguments | ||
101 | !**************************************************************************************** | ||
102 | REAL, DIMENSION(klon), INTENT(OUT) :: Acoef_U_out | ||
103 | REAL, DIMENSION(klon), INTENT(OUT) :: Acoef_V_out | ||
104 | REAL, DIMENSION(klon), INTENT(OUT) :: Bcoef_U_out | ||
105 | REAL, DIMENSION(klon), INTENT(OUT) :: Bcoef_V_out | ||
106 | |||
107 | !!! nrlmd le 02/05/2011 | ||
108 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Ccoef_U_out | ||
109 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Ccoef_V_out | ||
110 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Dcoef_U_out | ||
111 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Dcoef_V_out | ||
112 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Kcoef_m_out | ||
113 | REAL, DIMENSION(klon), INTENT(OUT) :: alf_1_out | ||
114 | REAL, DIMENSION(klon), INTENT(OUT) :: alf_2_out | ||
115 | !!! | ||
116 | |||
117 | ! Local variables | ||
118 | !**************************************************************************************** | ||
119 | REAL, DIMENSION(klon) :: u1lay, v1lay | ||
120 | INTEGER :: k, i | ||
121 | |||
122 | ! Include | ||
123 | !**************************************************************************************** | ||
124 | INCLUDE "YOMCST.h" | ||
125 | INCLUDE "compbl.h" | ||
126 | |||
127 | !**************************************************************************************** | ||
128 | ! Initialize module | ||
129 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 1919 times.
|
1920 | IF (firstcall) CALL climb_wind_init |
130 | |||
131 | !**************************************************************************************** | ||
132 | ! Calculate the coefficients C and D in : u(k) = C(k) + D(k)*u(k-1) | ||
133 | ! | ||
134 | !**************************************************************************************** | ||
135 | ! - Define alpha (alf1 and alf2) | ||
136 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | alf1(:) = 1.0 |
137 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | alf2(:) = 1.0 - alf1(:) |
138 | |||
139 | ! - Calculate the coefficients K | ||
140 |
4/4✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | Kcoefm(:,:) = 0.0 |
141 |
2/2✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 1920 times.
|
74880 | DO k = 2, klev |
142 |
2/2✓ Branch 0 taken 29961176 times.
✓ Branch 1 taken 72960 times.
|
30036056 | DO i=1,knon |
143 | Kcoefm(i,k) = coef_in(i,k)*RG*RG*dtime/(pplay(i,k-1)-pplay(i,k)) & | ||
144 | 72960 | *(paprs(i,k)*2/(temp(i,k)+temp(i,k-1))/RD)**2 | |
145 | END DO | ||
146 | END DO | ||
147 | |||
148 | ! - Calculate the coefficients C and D, component "u" | ||
149 | CALL calc_coef(knon, Kcoefm(:,:), delp(:,:), & | ||
150 | u_old(:,:), alf1(:), alf2(:), & | ||
151 | 1920 | Ccoef_U(:,:), Dcoef_U(:,:), Acoef_U(:), Bcoef_U(:)) | |
152 | |||
153 | ! - Calculate the coefficients C and D, component "v" | ||
154 | CALL calc_coef(knon, Kcoefm(:,:), delp(:,:), & | ||
155 | v_old(:,:), alf1(:), alf2(:), & | ||
156 | 1920 | Ccoef_V(:,:), Dcoef_V(:,:), Acoef_V(:), Bcoef_V(:)) | |
157 | |||
158 | !**************************************************************************************** | ||
159 | ! 6) | ||
160 | ! Return the first layer in output variables | ||
161 | ! | ||
162 | !**************************************************************************************** | ||
163 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Acoef_U_out = Acoef_U |
164 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Bcoef_U_out = Bcoef_U |
165 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Acoef_V_out = Acoef_V |
166 |
2/2✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
|
1910400 | Bcoef_V_out = Bcoef_V |
167 | |||
168 | !**************************************************************************************** | ||
169 | ! 7) | ||
170 | ! If Pbl is split, return also the other layers in output variables | ||
171 | ! | ||
172 | !**************************************************************************************** | ||
173 | !!! jyg le 07/02/2012 | ||
174 | !!jyg IF (mod(iflag_pbl_split,2) .eq.1) THEN | ||
175 | 1920 | IF (mod(iflag_pbl_split,10) .ge.1) THEN | |
176 | !!! nrlmd le 02/05/2011 | ||
177 | ✗ | DO k= 1, klev | |
178 | ✗ | DO i= 1, klon | |
179 | ✗ | Ccoef_U_out(i,k) = Ccoef_U(i,k) | |
180 | ✗ | Ccoef_V_out(i,k) = Ccoef_V(i,k) | |
181 | ✗ | Dcoef_U_out(i,k) = Dcoef_U(i,k) | |
182 | ✗ | Dcoef_V_out(i,k) = Dcoef_V(i,k) | |
183 | ✗ | Kcoef_m_out(i,k) = Kcoefm(i,k) | |
184 | ENDDO | ||
185 | ENDDO | ||
186 | ✗ | DO i= 1, klon | |
187 | ✗ | alf_1_out(i) = alf1(i) | |
188 | ✗ | alf_2_out(i) = alf2(i) | |
189 | ENDDO | ||
190 | !!! | ||
191 | ENDIF ! (mod(iflag_pbl_split,2) .ge.1) | ||
192 | !!! | ||
193 | |||
194 | 1920 | END SUBROUTINE climb_wind_down | |
195 | ! | ||
196 | !**************************************************************************************** | ||
197 | ! | ||
198 | 3840 | SUBROUTINE calc_coef(knon, Kcoef, delp, X, alfa1, alfa2, Ccoef, Dcoef, Acoef, Bcoef) | |
199 | ! | ||
200 | ! Find the coefficients C and D in fonction of alfa, K and delp | ||
201 | ! | ||
202 | ! Input arguments | ||
203 | !**************************************************************************************** | ||
204 | INTEGER, INTENT(IN) :: knon | ||
205 | REAL, DIMENSION(klon,klev), INTENT(IN) :: Kcoef, delp | ||
206 | REAL, DIMENSION(klon,klev), INTENT(IN) :: X | ||
207 | REAL, DIMENSION(klon), INTENT(IN) :: alfa1, alfa2 | ||
208 | |||
209 | ! Output arguments | ||
210 | !**************************************************************************************** | ||
211 | REAL, DIMENSION(klon), INTENT(OUT) :: Acoef, Bcoef | ||
212 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: Ccoef, Dcoef | ||
213 | |||
214 | ! local variables | ||
215 | !**************************************************************************************** | ||
216 | INTEGER :: k, i | ||
217 | REAL :: buf | ||
218 | |||
219 | INCLUDE "YOMCST.h" | ||
220 | !**************************************************************************************** | ||
221 | ! | ||
222 | |||
223 | ! Calculate coefficients C and D at top level, k=klev | ||
224 | ! | ||
225 |
4/4✓ Branch 0 taken 3840 times.
✓ Branch 1 taken 149760 times.
✓ Branch 2 taken 148861440 times.
✓ Branch 3 taken 149760 times.
|
149015040 | Ccoef(:,:) = 0.0 |
226 |
4/4✓ Branch 0 taken 149760 times.
✓ Branch 1 taken 3840 times.
✓ Branch 2 taken 148861440 times.
✓ Branch 3 taken 149760 times.
|
149015040 | Dcoef(:,:) = 0.0 |
227 | |||
228 |
2/2✓ Branch 0 taken 1576904 times.
✓ Branch 1 taken 3840 times.
|
1580744 | DO i = 1, knon |
229 | 1576904 | buf = delp(i,klev) + Kcoef(i,klev) | |
230 | |||
231 | 1576904 | Ccoef(i,klev) = X(i,klev)*delp(i,klev)/buf | |
232 | 1580744 | Dcoef(i,klev) = Kcoef(i,klev)/buf | |
233 | END DO | ||
234 | |||
235 | ! | ||
236 | ! Calculate coefficients C and D at top level (klev-1) <= k <= 2 | ||
237 | ! | ||
238 |
2/2✓ Branch 0 taken 142080 times.
✓ Branch 1 taken 3840 times.
|
145920 | DO k=(klev-1),2,-1 |
239 |
2/2✓ Branch 0 taken 58345448 times.
✓ Branch 1 taken 142080 times.
|
58491368 | DO i = 1, knon |
240 | 58345448 | buf = delp(i,k) + Kcoef(i,k) + Kcoef(i,k+1)*(1.-Dcoef(i,k+1)) | |
241 | |||
242 | 58345448 | Ccoef(i,k) = (X(i,k)*delp(i,k) + Kcoef(i,k+1)*Ccoef(i,k+1))/buf | |
243 | 58487528 | Dcoef(i,k) = Kcoef(i,k)/buf | |
244 | END DO | ||
245 | END DO | ||
246 | |||
247 | ! | ||
248 | ! Calculate coeffiecent A and B at surface | ||
249 | ! | ||
250 |
2/2✓ Branch 0 taken 1576904 times.
✓ Branch 1 taken 3840 times.
|
1580744 | DO i = 1, knon |
251 | 1576904 | buf = delp(i,1) + Kcoef(i,2)*(1-Dcoef(i,2)) | |
252 | 1576904 | Acoef(i) = (X(i,1)*delp(i,1) + Kcoef(i,2)*Ccoef(i,2))/buf | |
253 | 3840 | Bcoef(i) = -RG/buf | |
254 | END DO | ||
255 | |||
256 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1920 times.
|
29963096 | END SUBROUTINE calc_coef |
257 | ! | ||
258 | !**************************************************************************************** | ||
259 | ! | ||
260 | |||
261 | 1920 | SUBROUTINE climb_wind_up(knon, dtime, u_old, v_old, flx_u1, flx_v1, & | |
262 | !!! nrlmd le 02/05/2011 | ||
263 | 1920 | Acoef_U_in, Acoef_V_in, Bcoef_U_in, Bcoef_V_in, & | |
264 | 1920 | Ccoef_U_in, Ccoef_V_in, Dcoef_U_in, Dcoef_V_in, & | |
265 | Kcoef_m_in, & | ||
266 | !!! | ||
267 | flx_u_new, flx_v_new, d_u_new, d_v_new) | ||
268 | ! | ||
269 | ! Diffuse the wind components from the surface layer and up to the top layer. | ||
270 | ! Coefficents A, B, C and D are known from before. Start values for the diffusion are the | ||
271 | ! momentum fluxes at surface. | ||
272 | ! | ||
273 | ! u(k=1) = A + B*flx*dtime | ||
274 | ! u(k) = C(k) + D(k)*u(k-1) [2 <= k <= klev] | ||
275 | ! | ||
276 | !**************************************************************************************** | ||
277 | |||
278 | ! Input arguments | ||
279 | !**************************************************************************************** | ||
280 | INTEGER, INTENT(IN) :: knon | ||
281 | REAL, INTENT(IN) :: dtime | ||
282 | REAL, DIMENSION(klon,klev), INTENT(IN) :: u_old | ||
283 | REAL, DIMENSION(klon,klev), INTENT(IN) :: v_old | ||
284 | REAL, DIMENSION(klon), INTENT(IN) :: flx_u1, flx_v1 ! momentum flux | ||
285 | |||
286 | !!! nrlmd le 02/05/2011 | ||
287 | REAL, DIMENSION(klon), INTENT(IN) :: Acoef_U_in,Acoef_V_in, Bcoef_U_in, Bcoef_V_in | ||
288 | REAL, DIMENSION(klon,klev), INTENT(IN) :: Ccoef_U_in, Ccoef_V_in, Dcoef_U_in, Dcoef_V_in | ||
289 | REAL, DIMENSION(klon,klev), INTENT(IN) :: Kcoef_m_in | ||
290 | !!! | ||
291 | |||
292 | ! Output arguments | ||
293 | !**************************************************************************************** | ||
294 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: flx_u_new, flx_v_new | ||
295 | REAL, DIMENSION(klon,klev), INTENT(OUT) :: d_u_new, d_v_new | ||
296 | |||
297 | ! Local variables | ||
298 | !**************************************************************************************** | ||
299 | 3840 | REAL, DIMENSION(klon,klev) :: u_new, v_new | |
300 | INTEGER :: k, i | ||
301 | |||
302 | ! Include | ||
303 | !**************************************************************************************** | ||
304 | INCLUDE "YOMCST.h" | ||
305 | INCLUDE "compbl.h" | ||
306 | |||
307 | ! | ||
308 | !**************************************************************************************** | ||
309 | |||
310 | !!! jyg le 07/02/2012 | ||
311 | !!jyg IF (mod(iflag_pbl_split,2) .eq.1) THEN | ||
312 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1920 times.
|
1920 | IF (mod(iflag_pbl_split,10) .ge.1) THEN |
313 | !!! nrlmd le 02/05/2011 | ||
314 | ✗ | DO i = 1, knon | |
315 | ✗ | Acoef_U(i)=Acoef_U_in(i) | |
316 | ✗ | Acoef_V(i)=Acoef_V_in(i) | |
317 | ✗ | Bcoef_U(i)=Bcoef_U_in(i) | |
318 | ✗ | Bcoef_V(i)=Bcoef_V_in(i) | |
319 | ENDDO | ||
320 | ✗ | DO k = 1, klev | |
321 | ✗ | DO i = 1, knon | |
322 | ✗ | Ccoef_U(i,k)=Ccoef_U_in(i,k) | |
323 | ✗ | Ccoef_V(i,k)=Ccoef_V_in(i,k) | |
324 | ✗ | Dcoef_U(i,k)=Dcoef_U_in(i,k) | |
325 | ✗ | Dcoef_V(i,k)=Dcoef_V_in(i,k) | |
326 | ✗ | Kcoefm(i,k)=Kcoef_m_in(i,k) | |
327 | ENDDO | ||
328 | ENDDO | ||
329 | !!! | ||
330 | ENDIF ! (mod(iflag_pbl_split,2) .ge.1) | ||
331 | !!! | ||
332 | |||
333 | ! Niveau 1 | ||
334 |
2/2✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 788452 times.
|
790372 | DO i = 1, knon |
335 | 788452 | u_new(i,1) = Acoef_U(i) + Bcoef_U(i)*flx_u1(i)*dtime | |
336 | 790372 | v_new(i,1) = Acoef_V(i) + Bcoef_V(i)*flx_v1(i)*dtime | |
337 | END DO | ||
338 | |||
339 | ! Niveau 2 jusqu'au sommet klev | ||
340 |
2/2✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 1920 times.
|
74880 | DO k = 2, klev |
341 |
2/2✓ Branch 0 taken 29961176 times.
✓ Branch 1 taken 72960 times.
|
30036056 | DO i=1, knon |
342 | 29961176 | u_new(i,k) = Ccoef_U(i,k) + Dcoef_U(i,k) * u_new(i,k-1) | |
343 | 30034136 | v_new(i,k) = Ccoef_V(i,k) + Dcoef_V(i,k) * v_new(i,k-1) | |
344 | END DO | ||
345 | END DO | ||
346 | |||
347 | !**************************************************************************************** | ||
348 | ! Calcul flux | ||
349 | ! | ||
350 | !== flux_u/v est le flux de moment angulaire (positif vers bas) | ||
351 | !== dont l'unite est: (kg m/s)/(m**2 s) | ||
352 | ! | ||
353 | !**************************************************************************************** | ||
354 | ! | ||
355 |
4/4✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 74880 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | flx_u_new(:,:) = 0.0 |
356 |
4/4✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | flx_v_new(:,:) = 0.0 |
357 | |||
358 |
2/2✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 788452 times.
|
790372 | flx_u_new(1:knon,1)=flx_u1(1:knon) |
359 |
2/2✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 788452 times.
|
790372 | flx_v_new(1:knon,1)=flx_v1(1:knon) |
360 | |||
361 | ! Niveau 2->klev | ||
362 |
2/2✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 1920 times.
|
74880 | DO k = 2, klev |
363 |
2/2✓ Branch 0 taken 29961176 times.
✓ Branch 1 taken 72960 times.
|
30036056 | DO i = 1, knon |
364 | flx_u_new(i,k) = Kcoefm(i,k)/RG/dtime * & | ||
365 | 29961176 | (u_new(i,k)-u_new(i,k-1)) | |
366 | |||
367 | flx_v_new(i,k) = Kcoefm(i,k)/RG/dtime * & | ||
368 | 30034136 | (v_new(i,k)-v_new(i,k-1)) | |
369 | END DO | ||
370 | END DO | ||
371 | |||
372 | !**************************************************************************************** | ||
373 | ! Calcul tendances | ||
374 | ! | ||
375 | !**************************************************************************************** | ||
376 |
4/4✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 74880 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | d_u_new(:,:) = 0.0 |
377 |
4/4✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 74880 times.
✓ Branch 2 taken 74430720 times.
✓ Branch 3 taken 74880 times.
|
74507520 | d_v_new(:,:) = 0.0 |
378 |
2/2✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
|
76800 | DO k = 1, klev |
379 |
2/2✓ Branch 0 taken 30749628 times.
✓ Branch 1 taken 74880 times.
|
30826428 | DO i = 1, knon |
380 | 30749628 | d_u_new(i,k) = u_new(i,k) - u_old(i,k) | |
381 | 30824508 | d_v_new(i,k) = v_new(i,k) - v_old(i,k) | |
382 | END DO | ||
383 | END DO | ||
384 | |||
385 | 1576904 | END SUBROUTINE climb_wind_up | |
386 | ! | ||
387 | !**************************************************************************************** | ||
388 | ! | ||
389 | END MODULE climb_wind_mod | ||
390 |