Directory: | ./ |
---|---|
File: | phys/flott_gwd_rando_m.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 122 | 129 | 94.6% |
Branches: | 179 | 196 | 91.3% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | ! | ||
2 | ! $Id: flott_gwd_rando_m.F90 3531 2019-06-06 15:08:45Z fairhead $ | ||
3 | ! | ||
4 | module FLOTT_GWD_rando_m | ||
5 | |||
6 | implicit none | ||
7 | |||
8 | contains | ||
9 | |||
10 |
7/14✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 480 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 480 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 480 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 480 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 480 times.
✗ Branch 12 not taken.
✓ Branch 13 taken 480 times.
|
480 | SUBROUTINE FLOTT_GWD_rando(DTIME, pp, tt, uu, vv, prec, zustr, zvstr, d_u, & |
11 |
4/8✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 480 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 480 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 480 times.
|
480 | d_v,east_gwstress,west_gwstress) |
12 | |||
13 | ! Parametrization of the momentum flux deposition due to a discrete | ||
14 | ! number of gravity waves. | ||
15 | ! Author: F. Lott | ||
16 | ! July, 12th, 2012 | ||
17 | ! Gaussian distribution of the source, source is precipitation | ||
18 | ! Reference: Lott (JGR, vol 118, page 8897, 2013) | ||
19 | |||
20 | !ONLINE: | ||
21 | use dimphy, only: klon, klev | ||
22 | use assert_m, only: assert | ||
23 | USE ioipsl_getin_p_mod, ONLY : getin_p | ||
24 | USE vertical_layers_mod, ONLY : presnivs | ||
25 | CHARACTER (LEN=20) :: modname='flott_gwd_rando' | ||
26 | CHARACTER (LEN=80) :: abort_message | ||
27 | |||
28 | include "YOMCST.h" | ||
29 | include "clesphys.h" | ||
30 | ! OFFLINE: | ||
31 | ! include "dimensions.h" | ||
32 | ! include "dimphy.h" | ||
33 | ! END OF DIFFERENCE ONLINE-OFFLINE | ||
34 | include "YOEGWD.h" | ||
35 | |||
36 | ! 0. DECLARATIONS: | ||
37 | |||
38 | ! 0.1 INPUTS | ||
39 | REAL, intent(in)::DTIME ! Time step of the Physics | ||
40 | REAL, intent(in):: pp(:, :) ! (KLON, KLEV) Pressure at full levels | ||
41 | REAL, intent(in):: prec(:) ! (klon) Precipitation (kg/m^2/s) | ||
42 | REAL, intent(in):: TT(:, :) ! (KLON, KLEV) Temp at full levels | ||
43 | REAL, intent(in):: UU(:, :) ! (KLON, KLEV) Zonal wind at full levels | ||
44 | REAL, intent(in):: VV(:, :) ! (KLON, KLEV) Merid wind at full levels | ||
45 | |||
46 | ! 0.2 OUTPUTS | ||
47 | REAL, intent(out):: zustr(:), zvstr(:) ! (KLON) Surface Stresses | ||
48 | |||
49 | REAL, intent(inout):: d_u(:, :), d_v(:, :) | ||
50 | REAL, intent(inout):: east_gwstress(:, :) ! Profile of eastward stress | ||
51 | REAL, intent(inout):: west_gwstress(:, :) ! Profile of westward stress | ||
52 | |||
53 | ! (KLON, KLEV) tendencies on winds | ||
54 | |||
55 | ! O.3 INTERNAL ARRAYS | ||
56 | 960 | REAL BVLOW(klon) | |
57 | REAL DZ ! Characteristic depth of the Source | ||
58 | |||
59 | INTEGER II, JJ, LL | ||
60 | |||
61 | ! 0.3.0 TIME SCALE OF THE LIFE CYCLE OF THE WAVES PARAMETERIZED | ||
62 | |||
63 | REAL DELTAT | ||
64 | |||
65 | ! 0.3.1 GRAVITY-WAVES SPECIFICATIONS | ||
66 | |||
67 | INTEGER, PARAMETER:: NK = 2, NP = 2, NO = 2, NW = NK * NP * NO | ||
68 | INTEGER JK, JP, JO, JW | ||
69 | INTEGER, PARAMETER:: NA = 5 !number of realizations to get the phase speed | ||
70 | REAL KMIN, KMAX ! Min and Max horizontal wavenumbers | ||
71 | REAL CMAX ! standard deviation of the phase speed distribution | ||
72 | REAL RUWMAX,SAT ! ONLINE SPECIFIED IN run.def | ||
73 | REAL CPHA ! absolute PHASE VELOCITY frequency | ||
74 | 960 | REAL ZK(NW, KLON) ! Horizontal wavenumber amplitude | |
75 | 960 | REAL ZP(NW, KLON) ! Horizontal wavenumber angle | |
76 | 960 | REAL ZO(NW, KLON) ! Absolute frequency ! | |
77 | |||
78 | ! Waves Intr. freq. at the 1/2 lev surrounding the full level | ||
79 | 960 | REAL ZOM(NW, KLON), ZOP(NW, KLON) | |
80 | |||
81 | ! Wave EP-fluxes at the 2 semi levels surrounding the full level | ||
82 | 960 | REAL WWM(NW, KLON), WWP(NW, KLON) | |
83 | |||
84 | 960 | REAL RUW0(NW, KLON) ! Fluxes at launching level | |
85 | |||
86 | 960 | REAL RUWP(NW, KLON), RVWP(NW, KLON) | |
87 | ! Fluxes X and Y for each waves at 1/2 Levels | ||
88 | |||
89 | INTEGER LAUNCH, LTROP ! Launching altitude and tropo altitude | ||
90 | |||
91 | REAL XLAUNCH ! Controle the launching altitude | ||
92 | REAL XTROP ! SORT of Tropopause altitude | ||
93 | 960 | REAL RUW(KLON, KLEV + 1) ! Flux x at semi levels | |
94 | 960 | REAL RVW(KLON, KLEV + 1) ! Flux y at semi levels | |
95 | |||
96 | REAL PRMAX ! Maximum value of PREC, and for which our linear formula | ||
97 | ! for GWs parameterisation apply | ||
98 | |||
99 | ! 0.3.2 PARAMETERS OF WAVES DISSIPATIONS | ||
100 | |||
101 | REAL RDISS, ZOISEC ! COEFF DE DISSIPATION, SECURITY FOR INTRINSIC FREQ | ||
102 | |||
103 | ! 0.3.3 BACKGROUND FLOW AT 1/2 LEVELS AND VERTICAL COORDINATE | ||
104 | |||
105 | REAL H0 ! Characteristic Height of the atmosphere | ||
106 | REAL PR, TR ! Reference Pressure and Temperature | ||
107 | |||
108 | 960 | REAL ZH(KLON, KLEV + 1) ! Log-pressure altitude | |
109 | |||
110 | 960 | REAL UH(KLON, KLEV + 1), VH(KLON, KLEV + 1) ! Winds at 1/2 levels | |
111 | 960 | REAL PH(KLON, KLEV + 1) ! Pressure at 1/2 levels | |
112 | REAL PSEC ! Security to avoid division by 0 pressure | ||
113 | 960 | REAL BV(KLON, KLEV + 1) ! Brunt Vaisala freq. (BVF) at 1/2 levels | |
114 | REAL BVSEC ! Security to avoid negative BVF | ||
115 | REAL RAN_NUM_1,RAN_NUM_2,RAN_NUM_3 | ||
116 | |||
117 | 960 | REAL, DIMENSION(klev+1) ::HREF | |
118 | |||
119 | LOGICAL, SAVE :: gwd_reproductibilite_mpiomp=.true. | ||
120 | LOGICAL, SAVE :: firstcall = .TRUE. | ||
121 | !$OMP THREADPRIVATE(firstcall,gwd_reproductibilite_mpiomp) | ||
122 | |||
123 | |||
124 |
2/2✓ Branch 0 taken 1 times.
✓ Branch 1 taken 479 times.
|
480 | IF (firstcall) THEN |
125 | ! Cle introduite pour resoudre un probleme de non reproductibilite | ||
126 | ! Le but est de pouvoir tester de revenir a la version precedenete | ||
127 | ! A eliminer rapidement | ||
128 | 1 | CALL getin_p('gwd_reproductibilite_mpiomp',gwd_reproductibilite_mpiomp) | |
129 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (NW+3*NA>=KLEV) THEN |
130 | ✗ | abort_message = 'NW+3*NA>=KLEV Probleme pour generation des ondes' | |
131 | ✗ | CALL abort_physic (modname,abort_message,1) | |
132 | ENDIF | ||
133 | 1 | firstcall=.false. | |
134 | ENDIF | ||
135 | |||
136 | |||
137 | !----------------------------------------------------------------- | ||
138 | |||
139 | ! 1. INITIALISATIONS | ||
140 | |||
141 | ! 1.1 Basic parameter | ||
142 | |||
143 | ! Are provided from elsewhere (latent heat of vaporization, dry | ||
144 | ! gaz constant for air, gravity constant, heat capacity of dry air | ||
145 | ! at constant pressure, earth rotation rate, pi). | ||
146 | |||
147 | ! 1.2 Tuning parameters of V14 | ||
148 | |||
149 | |||
150 | RDISS = 0.5 ! Diffusion parameter | ||
151 | ! ONLINE | ||
152 | 480 | RUWMAX=GWD_RANDO_RUWMAX | |
153 | 480 | SAT=gwd_rando_sat | |
154 | !END ONLINE | ||
155 | ! OFFLINE | ||
156 | ! RUWMAX= 1.75 ! Launched flux | ||
157 | ! SAT=0.25 ! Saturation parameter | ||
158 | ! END OFFLINE | ||
159 | |||
160 | PRMAX = 20. / 24. /3600. | ||
161 | ! maximum of rain for which our theory applies (in kg/m^2/s) | ||
162 | |||
163 | ! Characteristic depth of the source | ||
164 | DZ = 1000. | ||
165 | XLAUNCH=0.5 ! Parameter that control launching altitude | ||
166 | XTROP=0.2 ! Parameter that control tropopause altitude | ||
167 | DELTAT=24.*3600. ! Time scale of the waves (first introduced in 9b) | ||
168 | ! OFFLINE | ||
169 | ! DELTAT=DTIME | ||
170 | ! END OFFLINE | ||
171 | |||
172 | KMIN = 2.E-5 | ||
173 | ! minimum horizontal wavenumber (inverse of the subgrid scale resolution) | ||
174 | |||
175 | KMAX = 1.E-3 ! Max horizontal wavenumber | ||
176 | CMAX = 30. ! Max phase speed velocity | ||
177 | |||
178 | TR = 240. ! Reference Temperature | ||
179 | PR = 101300. ! Reference pressure | ||
180 | 480 | H0 = RD * TR / RG ! Characteristic vertical scale height | |
181 | |||
182 | BVSEC = 5.E-3 ! Security to avoid negative BVF | ||
183 | PSEC = 1.E-6 ! Security to avoid division by 0 pressure | ||
184 | ZOISEC = 1.E-6 ! Security FOR 0 INTRINSIC FREQ | ||
185 | |||
186 | IF (1==0) THEN | ||
187 | !ONLINE | ||
188 | call assert(klon == (/size(pp, 1), size(tt, 1), size(uu, 1), & | ||
189 | size(vv, 1), size(zustr), size(zvstr), size(d_u, 1), & | ||
190 | size(d_v, 1), & | ||
191 | size(east_gwstress, 1), size(west_gwstress, 1) /), & | ||
192 | "FLOTT_GWD_RANDO klon") | ||
193 | call assert(klev == (/size(pp, 2), size(tt, 2), size(uu, 2), & | ||
194 | size(vv, 2), size(d_u, 2), size(d_v, 2), & | ||
195 | size(east_gwstress,2), size(west_gwstress,2) /), & | ||
196 | "FLOTT_GWD_RANDO klev") | ||
197 | !END ONLINE | ||
198 | ENDIF | ||
199 | |||
200 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | IF(DELTAT < DTIME)THEN |
201 | ✗ | abort_message='flott_gwd_rando: deltat < dtime!' | |
202 | ✗ | CALL abort_physic(modname,abort_message,1) | |
203 | ENDIF | ||
204 | |||
205 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | IF (KLEV < NW) THEN |
206 | ✗ | abort_message='flott_gwd_rando: you will have problem with random numbers' | |
207 | ✗ | CALL abort_physic(modname,abort_message,1) | |
208 | ENDIF | ||
209 | |||
210 | ! 2. EVALUATION OF THE BACKGROUND FLOW AT SEMI-LEVELS | ||
211 | |||
212 | ! Pressure and Inv of pressure | ||
213 |
2/2✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 480 times.
|
18720 | DO LL = 2, KLEV |
214 |
2/2✓ Branch 0 taken 18130560 times.
✓ Branch 1 taken 18240 times.
|
18149280 | PH(:, LL) = EXP((LOG(PP(:, LL)) + LOG(PP(:, LL - 1))) / 2.) |
215 | end DO | ||
216 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | PH(:, KLEV + 1) = 0. |
217 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | PH(:, 1) = 2. * PP(:, 1) - PH(:, 2) |
218 | |||
219 | ! Launching altitude | ||
220 | |||
221 | !Pour revenir a la version non reproductible en changeant le nombre de process | ||
222 |
1/2✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
|
480 | IF (gwd_reproductibilite_mpiomp) THEN |
223 | ! Reprend la formule qui calcule PH en fonction de PP=play | ||
224 |
2/2✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 480 times.
|
18720 | DO LL = 2, KLEV |
225 | 18720 | HREF(LL) = EXP((LOG(presnivs(LL)) + LOG(presnivs(LL - 1))) / 2.) | |
226 | end DO | ||
227 | 480 | HREF(KLEV + 1) = 0. | |
228 | 480 | HREF(1) = 2. * presnivs(1) - HREF(2) | |
229 | ELSE | ||
230 | ✗ | HREF(1:KLEV)=PH(KLON/2,1:KLEV) | |
231 | ENDIF | ||
232 | |||
233 | LAUNCH=0 | ||
234 | LTROP =0 | ||
235 |
2/2✓ Branch 0 taken 480 times.
✓ Branch 1 taken 18720 times.
|
19200 | DO LL = 1, KLEV |
236 |
2/2✓ Branch 0 taken 5760 times.
✓ Branch 1 taken 12960 times.
|
19200 | IF (HREF(LL) / HREF(1) > XLAUNCH) LAUNCH = LL |
237 | ENDDO | ||
238 |
2/2✓ Branch 0 taken 480 times.
✓ Branch 1 taken 18720 times.
|
19200 | DO LL = 1, KLEV |
239 |
2/2✓ Branch 0 taken 8160 times.
✓ Branch 1 taken 10560 times.
|
19200 | IF (HREF(LL) / HREF(1) > XTROP) LTROP = LL |
240 | ENDDO | ||
241 | !LAUNCH=22 ; LTROP=33 | ||
242 | ! print*,'LAUNCH=',LAUNCH,'LTROP=',LTROP | ||
243 | |||
244 | ! Log pressure vert. coordinate | ||
245 |
2/2✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
|
19680 | DO LL = 1, KLEV + 1 |
246 |
2/2✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
|
19104480 | ZH(:, LL) = H0 * LOG(PR / (PH(:, LL) + PSEC)) |
247 | end DO | ||
248 | |||
249 | ! BV frequency | ||
250 |
2/2✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 480 times.
|
18720 | DO LL = 2, KLEV |
251 | ! BVSEC: BV Frequency (UH USED IS AS A TEMPORARY ARRAY DOWN TO WINDS) | ||
252 | UH(:, LL) = 0.5 * (TT(:, LL) + TT(:, LL - 1)) & | ||
253 | * RD**2 / RCPD / H0**2 + (TT(:, LL) & | ||
254 |
2/2✓ Branch 0 taken 18130560 times.
✓ Branch 1 taken 18240 times.
|
18149280 | - TT(:, LL - 1)) / (ZH(:, LL) - ZH(:, LL - 1)) * RD / H0 |
255 | end DO | ||
256 | BVLOW(:) = 0.5 * (TT(:, LTROP )+ TT(:, LAUNCH)) & | ||
257 | * RD**2 / RCPD / H0**2 + (TT(:, LTROP ) & | ||
258 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | - TT(:, LAUNCH))/(ZH(:, LTROP )- ZH(:, LAUNCH)) * RD / H0 |
259 | |||
260 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | UH(:, 1) = UH(:, 2) |
261 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | UH(:, KLEV + 1) = UH(:, KLEV) |
262 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | BV(:, 1) = UH(:, 2) |
263 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | BV(:, KLEV + 1) = UH(:, KLEV) |
264 | ! SMOOTHING THE BV HELPS | ||
265 |
2/2✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 480 times.
|
18720 | DO LL = 2, KLEV |
266 |
2/2✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 18130560 times.
|
18149280 | BV(:, LL)=(UH(:, LL+1)+2.*UH(:, LL)+UH(:, LL-1))/4. |
267 | end DO | ||
268 | |||
269 |
4/4✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 19084800 times.
✓ Branch 3 taken 19200 times.
|
19104480 | BV=MAX(SQRT(MAX(BV, 0.)), BVSEC) |
270 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | BVLOW=MAX(SQRT(MAX(BVLOW, 0.)), BVSEC) |
271 | |||
272 | |||
273 | ! WINDS | ||
274 |
2/2✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 480 times.
|
18720 | DO LL = 2, KLEV |
275 |
2/2✓ Branch 0 taken 18130560 times.
✓ Branch 1 taken 18240 times.
|
18148800 | UH(:, LL) = 0.5 * (UU(:, LL) + UU(:, LL - 1)) ! Zonal wind |
276 |
2/2✓ Branch 0 taken 18130560 times.
✓ Branch 1 taken 18240 times.
|
18149280 | VH(:, LL) = 0.5 * (VV(:, LL) + VV(:, LL - 1)) ! Meridional wind |
277 | end DO | ||
278 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | UH(:, 1) = 0. |
279 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | VH(:, 1) = 0. |
280 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | UH(:, KLEV + 1) = UU(:, KLEV) |
281 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | VH(:, KLEV + 1) = VV(:, KLEV) |
282 | |||
283 | ! 3 WAVES CHARACTERISTICS CHOSEN RANDOMLY AT THE LAUNCH ALTITUDE | ||
284 | |||
285 | ! The mod functions of weird arguments are used to produce the | ||
286 | ! waves characteristics in an almost stochastic way | ||
287 | |||
288 |
2/2✓ Branch 0 taken 3840 times.
✓ Branch 1 taken 480 times.
|
4320 | DO JW = 1, NW |
289 | ! Angle | ||
290 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3821280 | DO II = 1, KLON |
291 | ! Angle (0 or PI so far) | ||
292 | 3816960 | RAN_NUM_1=MOD(TT(II, JW) * 10., 1.) | |
293 | 3816960 | RAN_NUM_2= MOD(TT(II, JW) * 100., 1.) | |
294 | ZP(JW, II) = (SIGN(1., 0.5 - RAN_NUM_1) + 1.) & | ||
295 | 3816960 | * RPI / 2. | |
296 | ! Horizontal wavenumber amplitude | ||
297 | 3816960 | ZK(JW, II) = KMIN + (KMAX - KMIN) *RAN_NUM_2 | |
298 | ! Horizontal phase speed | ||
299 | CPHA = 0. | ||
300 |
2/2✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 3816960 times.
|
22901760 | DO JJ = 1, NA |
301 | 19084800 | RAN_NUM_3=MOD(TT(II, JW+3*JJ)**2, 1.) | |
302 | CPHA = CPHA + & | ||
303 | 22901760 | CMAX*2.*(RAN_NUM_3 -0.5)*SQRT(3.)/SQRT(NA*1.) | |
304 | END DO | ||
305 |
2/2✓ Branch 0 taken 1910847 times.
✓ Branch 1 taken 1906113 times.
|
3816960 | IF (CPHA.LT.0.) THEN |
306 | 1910847 | CPHA = -1.*CPHA | |
307 | 1910847 | ZP(JW,II) = ZP(JW,II) + RPI | |
308 | ENDIF | ||
309 | ! Absolute frequency is imposed | ||
310 | 3816960 | ZO(JW, II) = CPHA * ZK(JW, II) | |
311 | ! Intrinsic frequency is imposed | ||
312 | ZO(JW, II) = ZO(JW, II) & | ||
313 | + ZK(JW, II) * COS(ZP(JW, II)) * UH(II, LAUNCH) & | ||
314 | 3816960 | + ZK(JW, II) * SIN(ZP(JW, II)) * VH(II, LAUNCH) | |
315 | ! Momentum flux at launch lev | ||
316 | 3820800 | RUW0(JW, II) = RUWMAX | |
317 | ENDDO | ||
318 | ENDDO | ||
319 | |||
320 | ! 4. COMPUTE THE FLUXES | ||
321 | |||
322 | ! 4.1 Vertical velocity at launching altitude to ensure | ||
323 | ! the correct value to the imposed fluxes. | ||
324 | |||
325 |
2/2✓ Branch 0 taken 3840 times.
✓ Branch 1 taken 480 times.
|
4320 | DO JW = 1, NW |
326 | |||
327 | ! Evaluate intrinsic frequency at launching altitude: | ||
328 | ZOP(JW, :) = ZO(JW, :) & | ||
329 | - ZK(JW, :) * COS(ZP(JW, :)) * UH(:, LAUNCH) & | ||
330 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3820800 | - ZK(JW, :) * SIN(ZP(JW, :)) * VH(:, LAUNCH) |
331 | |||
332 | ! VERSION WITH CONVECTIVE SOURCE | ||
333 | |||
334 | ! Vertical velocity at launch level, value to ensure the | ||
335 | ! imposed factor related to the convective forcing: | ||
336 | ! precipitations. | ||
337 | |||
338 | ! tanh limitation to values above prmax: | ||
339 | WWP(JW, :) = RUW0(JW, :) & | ||
340 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3820800 | * (RD / RCPD / H0 * RLVTT * PRMAX * TANH(PREC(:) / PRMAX))**2 |
341 | |||
342 | ! Factor related to the characteristics of the waves: | ||
343 | WWP(JW, :) = WWP(JW, :) * ZK(JW, :)**3 / KMIN / BVLOW(:) & | ||
344 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3820800 | / MAX(ABS(ZOP(JW, :)), ZOISEC)**3 |
345 | |||
346 | ! Moderation by the depth of the source (dz here): | ||
347 | WWP(JW, :) = WWP(JW, :) & | ||
348 | * EXP(- BVLOW(:)**2 / MAX(ABS(ZOP(JW, :)), ZOISEC)**2 * ZK(JW, :)**2 & | ||
349 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3820800 | * DZ**2) |
350 | |||
351 | ! Put the stress in the right direction: | ||
352 | RUWP(JW, :) = ZOP(JW, :) / MAX(ABS(ZOP(JW, :)), ZOISEC)**2 & | ||
353 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3820800 | * BV(:, LAUNCH) * COS(ZP(JW, :)) * WWP(JW, :)**2 |
354 | RVWP(JW, :) = ZOP(JW, :) / MAX(ABS(ZOP(JW, :)), ZOISEC)**2 & | ||
355 |
2/2✓ Branch 0 taken 3816960 times.
✓ Branch 1 taken 3840 times.
|
3821280 | * BV(:, LAUNCH) * SIN(ZP(JW, :)) * WWP(JW, :)**2 |
356 | end DO | ||
357 | |||
358 | |||
359 | ! 4.2 Uniform values below the launching altitude | ||
360 | |||
361 |
2/2✓ Branch 0 taken 5760 times.
✓ Branch 1 taken 480 times.
|
6240 | DO LL = 1, LAUNCH |
362 |
2/2✓ Branch 0 taken 5725440 times.
✓ Branch 1 taken 5760 times.
|
5731200 | RUW(:, LL) = 0 |
363 |
2/2✓ Branch 0 taken 5725440 times.
✓ Branch 1 taken 5760 times.
|
5731200 | RVW(:, LL) = 0 |
364 |
2/2✓ Branch 0 taken 46080 times.
✓ Branch 1 taken 5760 times.
|
52320 | DO JW = 1, NW |
365 |
2/2✓ Branch 0 taken 45803520 times.
✓ Branch 1 taken 46080 times.
|
45849600 | RUW(:, LL) = RUW(:, LL) + RUWP(JW, :) |
366 |
2/2✓ Branch 0 taken 45803520 times.
✓ Branch 1 taken 46080 times.
|
45855360 | RVW(:, LL) = RVW(:, LL) + RVWP(JW, :) |
367 | end DO | ||
368 | end DO | ||
369 | |||
370 | ! 4.3 Loop over altitudes, with passage from one level to the next | ||
371 | ! done by i) conserving the EP flux, ii) dissipating a little, | ||
372 | ! iii) testing critical levels, and vi) testing the breaking. | ||
373 | |||
374 |
2/2✓ Branch 0 taken 12960 times.
✓ Branch 1 taken 480 times.
|
13440 | DO LL = LAUNCH, KLEV - 1 |
375 | ! Warning: all the physics is here (passage from one level | ||
376 | ! to the next) | ||
377 |
2/2✓ Branch 0 taken 103680 times.
✓ Branch 1 taken 12960 times.
|
116640 | DO JW = 1, NW |
378 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | ZOM(JW, :) = ZOP(JW, :) |
379 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | WWM(JW, :) = WWP(JW, :) |
380 | ! Intrinsic Frequency | ||
381 | ZOP(JW, :) = ZO(JW, :) - ZK(JW, :) * COS(ZP(JW, :)) * UH(:, LL + 1) & | ||
382 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | - ZK(JW, :) * SIN(ZP(JW, :)) * VH(:, LL + 1) |
383 | |||
384 | ! No breaking (Eq.6) | ||
385 | ! Dissipation (Eq. 8) | ||
386 | WWP(JW, :) = WWM(JW, :) * EXP(- 4. * RDISS * PR / (PH(:, LL + 1) & | ||
387 | + PH(:, LL)) * ((BV(:, LL + 1) + BV(:, LL)) / 2.)**3 & | ||
388 | / MAX(ABS(ZOP(JW, :) + ZOM(JW, :)) / 2., ZOISEC)**4 & | ||
389 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | * ZK(JW, :)**3 * (ZH(:, LL + 1) - ZH(:, LL))) |
390 | |||
391 | ! Critical levels (forced to zero if intrinsic frequency changes sign) | ||
392 | ! Saturation (Eq. 12) | ||
393 | WWP(JW, :) = min(WWP(JW, :), MAX(0., & | ||
394 | SIGN(1., ZOP(JW, :) * ZOM(JW, :))) * ABS(ZOP(JW, :))**3 & | ||
395 | / BV(:, LL + 1) * EXP(- ZH(:, LL + 1) / H0) * KMIN**2 & | ||
396 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103174560 | * SAT**2 / ZK(JW, :)**4) |
397 | end DO | ||
398 | |||
399 | ! Evaluate EP-flux from Eq. 7 and give the right orientation to | ||
400 | ! the stress | ||
401 | |||
402 |
2/2✓ Branch 0 taken 103680 times.
✓ Branch 1 taken 12960 times.
|
116640 | DO JW = 1, NW |
403 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | RUWP(JW, :) = SIGN(1., ZOP(JW, :))*COS(ZP(JW, :)) * WWP(JW, :) |
404 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103174560 | RVWP(JW, :) = SIGN(1., ZOP(JW, :))*SIN(ZP(JW, :)) * WWP(JW, :) |
405 | end DO | ||
406 | |||
407 |
2/2✓ Branch 0 taken 12882240 times.
✓ Branch 1 taken 12960 times.
|
12895200 | RUW(:, LL + 1) = 0. |
408 |
2/2✓ Branch 0 taken 12882240 times.
✓ Branch 1 taken 12960 times.
|
12895200 | RVW(:, LL + 1) = 0. |
409 | |||
410 |
2/2✓ Branch 0 taken 12960 times.
✓ Branch 1 taken 103680 times.
|
117120 | DO JW = 1, NW |
411 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | RUW(:, LL + 1) = RUW(:, LL + 1) + RUWP(JW, :) |
412 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | RVW(:, LL + 1) = RVW(:, LL + 1) + RVWP(JW, :) |
413 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103161600 | EAST_GWSTRESS(:, LL)=EAST_GWSTRESS(:, LL)+MAX(0.,RUWP(JW,:))/FLOAT(NW) |
414 |
2/2✓ Branch 0 taken 103057920 times.
✓ Branch 1 taken 103680 times.
|
103174560 | WEST_GWSTRESS(:, LL)=WEST_GWSTRESS(:, LL)+MIN(0.,RUWP(JW,:))/FLOAT(NW) |
415 | end DO | ||
416 | end DO | ||
417 | ! OFFLINE ONLY | ||
418 | ! PRINT *,'SAT PROFILE:' | ||
419 | ! DO LL=1,KLEV | ||
420 | ! PRINT *,ZH(KLON/2,LL)/1000.,SAT*(2.+TANH(ZH(KLON/2,LL)/H0-8.)) | ||
421 | ! ENDDO | ||
422 | |||
423 | ! 5 CALCUL DES TENDANCES: | ||
424 | |||
425 | ! 5.1 Rectification des flux au sommet et dans les basses couches | ||
426 | |||
427 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | RUW(:, KLEV + 1) = 0. |
428 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | RVW(:, KLEV + 1) = 0. |
429 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | RUW(:, 1) = RUW(:, LAUNCH) |
430 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | RVW(:, 1) = RVW(:, LAUNCH) |
431 |
2/2✓ Branch 0 taken 5760 times.
✓ Branch 1 taken 480 times.
|
6240 | DO LL = 1, LAUNCH |
432 |
2/2✓ Branch 0 taken 5725440 times.
✓ Branch 1 taken 5760 times.
|
5731200 | RUW(:, LL) = RUW(:, LAUNCH+1) |
433 |
2/2✓ Branch 0 taken 5725440 times.
✓ Branch 1 taken 5760 times.
|
5731200 | RVW(:, LL) = RVW(:, LAUNCH+1) |
434 |
2/2✓ Branch 0 taken 5725440 times.
✓ Branch 1 taken 5760 times.
|
5731200 | EAST_GWSTRESS(:, LL) = EAST_GWSTRESS(:, LAUNCH) |
435 |
2/2✓ Branch 0 taken 5725440 times.
✓ Branch 1 taken 5760 times.
|
5731680 | WEST_GWSTRESS(:, LL) = WEST_GWSTRESS(:, LAUNCH) |
436 | end DO | ||
437 | |||
438 | ! AR-1 RECURSIVE FORMULA (13) IN VERSION 4 | ||
439 |
2/2✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
|
19200 | DO LL = 1, KLEV |
440 | D_U(:, LL) = (1.-DTIME/DELTAT) * D_U(:, LL) + DTIME/DELTAT/REAL(NW) * & | ||
441 | RG * (RUW(:, LL + 1) - RUW(:, LL)) & | ||
442 |
2/2✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
|
18626400 | / (PH(:, LL + 1) - PH(:, LL)) * DTIME |
443 | ! NO AR-1 FOR MERIDIONAL TENDENCIES | ||
444 | D_V(:, LL) = 1./REAL(NW) * & | ||
445 | RG * (RVW(:, LL + 1) - RVW(:, LL)) & | ||
446 |
2/2✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 18607680 times.
|
18626880 | / (PH(:, LL + 1) - PH(:, LL)) * DTIME |
447 | ENDDO | ||
448 | |||
449 | ! Cosmetic: evaluation of the cumulated stress | ||
450 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | ZUSTR = 0. |
451 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | ZVSTR = 0. |
452 |
2/2✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
|
19200 | DO LL = 1, KLEV |
453 |
2/2✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
|
18626400 | ZUSTR = ZUSTR + D_U(:, LL) / RG * (PH(:, LL + 1) - PH(:, LL))/DTIME |
454 |
2/2✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 18607680 times.
|
18626880 | ZVSTR = ZVSTR + D_V(:, LL) / RG * (PH(:, LL + 1) - PH(:, LL))/DTIME |
455 | ENDDO | ||
456 | |||
457 | |||
458 | 480 | END SUBROUTINE FLOTT_GWD_RANDO | |
459 | |||
460 | end module FLOTT_GWD_rando_m | ||
461 |