Neko 0.9.99
A portable framework for high-order spectral element flow simulations
Loading...
Searching...
No Matches
sx_set_convect_rst.f90
Go to the documentation of this file.
1! Copyright (c) 2024, The Neko Authors
2! All rights reserved.
3!
4! Redistribution and use in source and binary forms, with or without
5! modification, are permitted provided that the following conditions
6! are met:
7!
8! * Redistributions of source code must retain the above copyright
9! notice, this list of conditions and the following disclaimer.
10!
11! * Redistributions in binary form must reproduce the above
12! copyright notice, this list of conditions and the following
13! disclaimer in the documentation and/or other materials provided
14! with the distribution.
15!
16! * Neither the name of the authors nor the names of its
17! contributors may be used to endorse or promote products derived
18! from this software without specific prior written permission.
19!
20! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
21! "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
22! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
23! FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
24! COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
25! INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
26! BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
27! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
28! CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
29! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
30! ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
31! POSSIBILITY OF SUCH DAMAGE.
32!
33submodule(opr_sx) sx_set_convect_rst
34 use num_types, only : rp
35 implicit none
36
37contains
38
39 module subroutine opr_sx_set_convect_rst(cr, cs, ct, cx, cy, cz, xh, coef)
40 type(space_t), intent(inout) :: Xh
41 type(coef_t), intent(inout) :: coef
42 real(kind=rp), dimension(Xh%lxyz, coef%msh%nelv), &
43 intent(inout) :: cr, cs, ct
44 real(kind=rp), dimension(Xh%lxyz, coef%msh%nelv), &
45 intent(in) :: cx, cy, cz
46 associate(drdx => coef%drdx, drdy => coef%drdy, drdz => coef%drdz, &
47 dsdx => coef%dsdx, dsdy => coef%dsdy, dsdz => coef%dsdz, &
48 dtdx => coef%dtdx, dtdy => coef%dtdy, dtdz => coef%dtdz, &
49 nelv => coef%msh%nelv, lx => xh%lx, w3 => xh%w3)
50
51 select case (lx)
52 case (18)
53 call sx_set_convect_rst_lx18(cr, cs, ct, cx, cy, cz, drdx, &
54 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
55 case (17)
56 call sx_set_convect_rst_lx17(cr, cs, ct, cx, cy, cz, drdx, &
57 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
58 case (16)
59 call sx_set_convect_rst_lx16(cr, cs, ct, cx, cy, cz, drdx, &
60 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
61 case (15)
62 call sx_set_convect_rst_lx15(cr, cs, ct, cx, cy, cz, drdx, &
63 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
64 case (14)
65 call sx_set_convect_rst_lx14(cr, cs, ct, cx, cy, cz, drdx, &
66 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
67 case (13)
68 call sx_set_convect_rst_lx13(cr, cs, ct, cx, cy, cz, drdx, &
69 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
70 case (12)
71 call sx_set_convect_rst_lx12(cr, cs, ct, cx, cy, cz, drdx, &
72 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
73 case (11)
74 call sx_set_convect_rst_lx11(cr, cs, ct, cx, cy, cz, drdx, &
75 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
76 case (10)
77 call sx_set_convect_rst_lx10(cr, cs, ct, cx, cy, cz, drdx, &
78 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
79 case (9)
80 call sx_set_convect_rst_lx9(cr, cs, ct, cx, cy, cz, drdx, &
81 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
82 case (8)
83 call sx_set_convect_rst_lx8(cr, cs, ct, cx, cy, cz, drdx, &
84 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
85 case (7)
86 call sx_set_convect_rst_lx7(cr, cs, ct, cx, cy, cz, drdx, &
87 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
88 case (6)
89 call sx_set_convect_rst_lx6(cr, cs, ct, cx, cy, cz, drdx, &
90 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
91 case (5)
92 call sx_set_convect_rst_lx5(cr, cs, ct, cx, cy, cz, drdx, &
93 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
94 case (4)
95 call sx_set_convect_rst_lx4(cr, cs, ct, cx, cy, cz, drdx, &
96 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
97 case (3)
98 call sx_set_convect_rst_lx3(cr, cs, ct, cx, cy, cz, drdx, &
99 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
100 case (2)
101 call sx_set_convect_rst_lx2(cr, cs, ct, cx, cy, cz, drdx, &
102 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv)
103 case default
104 call sx_set_convect_rst_lx(cr, cs, ct, cx, cy, cz, drdx, &
105 dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, nelv, lx)
106 end select
107 end associate
108 end subroutine opr_sx_set_convect_rst
109
110 subroutine sx_set_convect_rst_lx(cr, cs, ct, cx, cy, cz, &
111 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n, lx)
112 integer, intent(in) :: n, lx
113 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
114 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
115 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
116 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
117 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
118 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
119 integer :: e, i
120
121 do i = 1, lx * lx * lx
122 do e = 1, n
123 cr(i,1,1,e) = w3(i,1,1) &
124 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
125 + cy(i,1,1,e) * drdy(i,1,1,e) &
126 + cz(i,1,1,e) * drdz(i,1,1,e) )
127 cs(i,1,1,e) = w3(i,1,1) &
128 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
129 + cy(i,1,1,e) * dsdy(i,1,1,e) &
130 + cz(i,1,1,e) * dsdz(i,1,1,e))
131 ct(i,1,1,e) = w3(i,1,1) &
132 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
133 + cy(i,1,1,e) * dtdy(i,1,1,e) &
134 + cz(i,1,1,e) * dtdz(i,1,1,e))
135 end do
136 end do
137
138 end subroutine sx_set_convect_rst_lx
139
140 subroutine sx_set_convect_rst_lx18(cr, cs, ct, cx, cy, cz, &
141 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
142 integer, parameter :: lx = 18
143 integer, intent(in) :: n
144 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
145 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
146 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
147 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
148 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
149 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
150 integer :: e, i
151
152 do i = 1, lx * lx * lx
153 do e = 1, n
154 cr(i,1,1,e) = w3(i,1,1) &
155 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
156 + cy(i,1,1,e) * drdy(i,1,1,e) &
157 + cz(i,1,1,e) * drdz(i,1,1,e) )
158 cs(i,1,1,e) = w3(i,1,1) &
159 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
160 + cy(i,1,1,e) * dsdy(i,1,1,e) &
161 + cz(i,1,1,e) * dsdz(i,1,1,e))
162 ct(i,1,1,e) = w3(i,1,1) &
163 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
164 + cy(i,1,1,e) * dtdy(i,1,1,e) &
165 + cz(i,1,1,e) * dtdz(i,1,1,e))
166 end do
167 end do
168
169 end subroutine sx_set_convect_rst_lx18
170
171 subroutine sx_set_convect_rst_lx17(cr, cs, ct, cx, cy, cz, &
172 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
173 integer, parameter :: lx = 17
174 integer, intent(in) :: n
175 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
176 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
177 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
178 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
179 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
180 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
181 integer :: e, i
182
183 do i = 1, lx * lx * lx
184 do e = 1, n
185 cr(i,1,1,e) = w3(i,1,1) &
186 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
187 + cy(i,1,1,e) * drdy(i,1,1,e) &
188 + cz(i,1,1,e) * drdz(i,1,1,e) )
189 cs(i,1,1,e) = w3(i,1,1) &
190 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
191 + cy(i,1,1,e) * dsdy(i,1,1,e) &
192 + cz(i,1,1,e) * dsdz(i,1,1,e))
193 ct(i,1,1,e) = w3(i,1,1) &
194 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
195 + cy(i,1,1,e) * dtdy(i,1,1,e) &
196 + cz(i,1,1,e) * dtdz(i,1,1,e))
197 end do
198 end do
199
200 end subroutine sx_set_convect_rst_lx17
201
202 subroutine sx_set_convect_rst_lx16(cr, cs, ct, cx, cy, cz, &
203 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
204 integer, parameter :: lx = 16
205 integer, intent(in) :: n
206 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
207 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
208 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
209 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
210 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
211 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
212 integer :: e, i
213
214 do i = 1, lx * lx * lx
215 do e = 1, n
216 cr(i,1,1,e) = w3(i,1,1) &
217 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
218 + cy(i,1,1,e) * drdy(i,1,1,e) &
219 + cz(i,1,1,e) * drdz(i,1,1,e) )
220 cs(i,1,1,e) = w3(i,1,1) &
221 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
222 + cy(i,1,1,e) * dsdy(i,1,1,e) &
223 + cz(i,1,1,e) * dsdz(i,1,1,e))
224 ct(i,1,1,e) = w3(i,1,1) &
225 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
226 + cy(i,1,1,e) * dtdy(i,1,1,e) &
227 + cz(i,1,1,e) * dtdz(i,1,1,e))
228 end do
229 end do
230
231 end subroutine sx_set_convect_rst_lx16
232
233 subroutine sx_set_convect_rst_lx15(cr, cs, ct, cx, cy, cz, &
234 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
235 integer, parameter :: lx = 15
236 integer, intent(in) :: n
237 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
238 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
239 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
240 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
241 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
242 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
243 integer :: e, i
244
245 do i = 1, lx * lx * lx
246 do e = 1, n
247 cr(i,1,1,e) = w3(i,1,1) &
248 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
249 + cy(i,1,1,e) * drdy(i,1,1,e) &
250 + cz(i,1,1,e) * drdz(i,1,1,e) )
251 cs(i,1,1,e) = w3(i,1,1) &
252 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
253 + cy(i,1,1,e) * dsdy(i,1,1,e) &
254 + cz(i,1,1,e) * dsdz(i,1,1,e))
255 ct(i,1,1,e) = w3(i,1,1) &
256 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
257 + cy(i,1,1,e) * dtdy(i,1,1,e) &
258 + cz(i,1,1,e) * dtdz(i,1,1,e))
259 end do
260 end do
261
262 end subroutine sx_set_convect_rst_lx15
263
264 subroutine sx_set_convect_rst_lx14(cr, cs, ct, cx, cy, cz, &
265 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
266 integer, parameter :: lx = 14
267 integer, intent(in) :: n
268 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
269 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
270 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
271 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
272 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
273 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
274 integer :: e, i
275
276 do i = 1, lx * lx * lx
277 do e = 1, n
278 cr(i,1,1,e) = w3(i,1,1) &
279 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
280 + cy(i,1,1,e) * drdy(i,1,1,e) &
281 + cz(i,1,1,e) * drdz(i,1,1,e) )
282 cs(i,1,1,e) = w3(i,1,1) &
283 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
284 + cy(i,1,1,e) * dsdy(i,1,1,e) &
285 + cz(i,1,1,e) * dsdz(i,1,1,e))
286 ct(i,1,1,e) = w3(i,1,1) &
287 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
288 + cy(i,1,1,e) * dtdy(i,1,1,e) &
289 + cz(i,1,1,e) * dtdz(i,1,1,e))
290 end do
291 end do
292
293 end subroutine sx_set_convect_rst_lx14
294
295 subroutine sx_set_convect_rst_lx13(cr, cs, ct, cx, cy, cz, &
296 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
297 integer, parameter :: lx = 13
298 integer, intent(in) :: n
299 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
300 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
301 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
302 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
303 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
304 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
305 integer :: e, i
306
307 do i = 1, lx * lx * lx
308 do e = 1, n
309 cr(i,1,1,e) = w3(i,1,1) &
310 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
311 + cy(i,1,1,e) * drdy(i,1,1,e) &
312 + cz(i,1,1,e) * drdz(i,1,1,e) )
313 cs(i,1,1,e) = w3(i,1,1) &
314 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
315 + cy(i,1,1,e) * dsdy(i,1,1,e) &
316 + cz(i,1,1,e) * dsdz(i,1,1,e))
317 ct(i,1,1,e) = w3(i,1,1) &
318 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
319 + cy(i,1,1,e) * dtdy(i,1,1,e) &
320 + cz(i,1,1,e) * dtdz(i,1,1,e))
321 end do
322 end do
323
324 end subroutine sx_set_convect_rst_lx13
325
326 subroutine sx_set_convect_rst_lx12(cr, cs, ct, cx, cy, cz, &
327 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
328 integer, parameter :: lx = 12
329 integer, intent(in) :: n
330 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
331 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
332 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
333 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
334 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
335 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
336 integer :: e, i
337
338 do i = 1, lx * lx * lx
339 do e = 1, n
340 cr(i,1,1,e) = w3(i,1,1) &
341 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
342 + cy(i,1,1,e) * drdy(i,1,1,e) &
343 + cz(i,1,1,e) * drdz(i,1,1,e) )
344 cs(i,1,1,e) = w3(i,1,1) &
345 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
346 + cy(i,1,1,e) * dsdy(i,1,1,e) &
347 + cz(i,1,1,e) * dsdz(i,1,1,e))
348 ct(i,1,1,e) = w3(i,1,1) &
349 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
350 + cy(i,1,1,e) * dtdy(i,1,1,e) &
351 + cz(i,1,1,e) * dtdz(i,1,1,e))
352 end do
353 end do
354
355 end subroutine sx_set_convect_rst_lx12
356
357 subroutine sx_set_convect_rst_lx11(cr, cs, ct, cx, cy, cz, &
358 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
359 integer, parameter :: lx = 11
360 integer, intent(in) :: n
361 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
362 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
363 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
364 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
365 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
366 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
367 integer :: e, i
368
369 do i = 1, lx * lx * lx
370 do e = 1, n
371 cr(i,1,1,e) = w3(i,1,1) &
372 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
373 + cy(i,1,1,e) * drdy(i,1,1,e) &
374 + cz(i,1,1,e) * drdz(i,1,1,e) )
375 cs(i,1,1,e) = w3(i,1,1) &
376 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
377 + cy(i,1,1,e) * dsdy(i,1,1,e) &
378 + cz(i,1,1,e) * dsdz(i,1,1,e))
379 ct(i,1,1,e) = w3(i,1,1) &
380 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
381 + cy(i,1,1,e) * dtdy(i,1,1,e) &
382 + cz(i,1,1,e) * dtdz(i,1,1,e))
383 end do
384 end do
385
386 end subroutine sx_set_convect_rst_lx11
387
388 subroutine sx_set_convect_rst_lx10(cr, cs, ct, cx, cy, cz, &
389 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
390 integer, parameter :: lx = 10
391 integer, intent(in) :: n
392 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
393 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
394 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
395 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
396 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
397 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
398 integer :: e, i
399
400 do i = 1, lx * lx * lx
401 do e = 1, n
402 cr(i,1,1,e) = w3(i,1,1) &
403 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
404 + cy(i,1,1,e) * drdy(i,1,1,e) &
405 + cz(i,1,1,e) * drdz(i,1,1,e) )
406 cs(i,1,1,e) = w3(i,1,1) &
407 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
408 + cy(i,1,1,e) * dsdy(i,1,1,e) &
409 + cz(i,1,1,e) * dsdz(i,1,1,e))
410 ct(i,1,1,e) = w3(i,1,1) &
411 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
412 + cy(i,1,1,e) * dtdy(i,1,1,e) &
413 + cz(i,1,1,e) * dtdz(i,1,1,e))
414 end do
415 end do
416
417 end subroutine sx_set_convect_rst_lx10
418
419 subroutine sx_set_convect_rst_lx9(cr, cs, ct, cx, cy, cz, &
420 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
421 integer, parameter :: lx = 9
422 integer, intent(in) :: n
423 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
424 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
425 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
426 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
427 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
428 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
429 integer :: e, i
430
431 do i = 1, lx * lx * lx
432 do e = 1, n
433 cr(i,1,1,e) = w3(i,1,1) &
434 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
435 + cy(i,1,1,e) * drdy(i,1,1,e) &
436 + cz(i,1,1,e) * drdz(i,1,1,e) )
437 cs(i,1,1,e) = w3(i,1,1) &
438 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
439 + cy(i,1,1,e) * dsdy(i,1,1,e) &
440 + cz(i,1,1,e) * dsdz(i,1,1,e))
441 ct(i,1,1,e) = w3(i,1,1) &
442 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
443 + cy(i,1,1,e) * dtdy(i,1,1,e) &
444 + cz(i,1,1,e) * dtdz(i,1,1,e))
445 end do
446 end do
447
448 end subroutine sx_set_convect_rst_lx9
449
450 subroutine sx_set_convect_rst_lx8(cr, cs, ct, cx, cy, cz, &
451 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
452 integer, parameter :: lx = 8
453 integer, intent(in) :: n
454 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
455 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
456 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
457 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
458 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
459 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
460 integer :: e, i
461
462 do i = 1, lx * lx * lx
463 do e = 1, n
464 cr(i,1,1,e) = w3(i,1,1) &
465 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
466 + cy(i,1,1,e) * drdy(i,1,1,e) &
467 + cz(i,1,1,e) * drdz(i,1,1,e) )
468 cs(i,1,1,e) = w3(i,1,1) &
469 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
470 + cy(i,1,1,e) * dsdy(i,1,1,e) &
471 + cz(i,1,1,e) * dsdz(i,1,1,e))
472 ct(i,1,1,e) = w3(i,1,1) &
473 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
474 + cy(i,1,1,e) * dtdy(i,1,1,e) &
475 + cz(i,1,1,e) * dtdz(i,1,1,e))
476 end do
477 end do
478
479 end subroutine sx_set_convect_rst_lx8
480
481 subroutine sx_set_convect_rst_lx7(cr, cs, ct, cx, cy, cz, &
482 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
483 integer, parameter :: lx = 7
484 integer, intent(in) :: n
485 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
486 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
487 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
488 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
489 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
490 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
491 integer :: e, i
492
493 do i = 1, lx * lx * lx
494 do e = 1, n
495 cr(i,1,1,e) = w3(i,1,1) &
496 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
497 + cy(i,1,1,e) * drdy(i,1,1,e) &
498 + cz(i,1,1,e) * drdz(i,1,1,e) )
499 cs(i,1,1,e) = w3(i,1,1) &
500 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
501 + cy(i,1,1,e) * dsdy(i,1,1,e) &
502 + cz(i,1,1,e) * dsdz(i,1,1,e))
503 ct(i,1,1,e) = w3(i,1,1) &
504 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
505 + cy(i,1,1,e) * dtdy(i,1,1,e) &
506 + cz(i,1,1,e) * dtdz(i,1,1,e))
507 end do
508 end do
509
510 end subroutine sx_set_convect_rst_lx7
511
512 subroutine sx_set_convect_rst_lx6(cr, cs, ct, cx, cy, cz, &
513 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
514 integer, parameter :: lx = 6
515 integer, intent(in) :: n
516 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
517 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
518 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
519 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
520 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
521 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
522 integer :: e, i
523
524 do i = 1, lx * lx * lx
525 do e = 1, n
526 cr(i,1,1,e) = w3(i,1,1) &
527 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
528 + cy(i,1,1,e) * drdy(i,1,1,e) &
529 + cz(i,1,1,e) * drdz(i,1,1,e) )
530 cs(i,1,1,e) = w3(i,1,1) &
531 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
532 + cy(i,1,1,e) * dsdy(i,1,1,e) &
533 + cz(i,1,1,e) * dsdz(i,1,1,e))
534 ct(i,1,1,e) = w3(i,1,1) &
535 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
536 + cy(i,1,1,e) * dtdy(i,1,1,e) &
537 + cz(i,1,1,e) * dtdz(i,1,1,e))
538 end do
539 end do
540
541end subroutine sx_set_convect_rst_lx6
542
543subroutine sx_set_convect_rst_lx5(cr, cs, ct, cx, cy, cz, &
544 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
545 integer, parameter :: lx = 5
546 integer, intent(in) :: n
547 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
548 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
549 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
550 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
551 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
552 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
553 integer :: e, i
554
555 do i = 1, lx * lx * lx
556 do e = 1, n
557 cr(i,1,1,e) = w3(i,1,1) &
558 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
559 + cy(i,1,1,e) * drdy(i,1,1,e) &
560 + cz(i,1,1,e) * drdz(i,1,1,e) )
561 cs(i,1,1,e) = w3(i,1,1) &
562 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
563 + cy(i,1,1,e) * dsdy(i,1,1,e) &
564 + cz(i,1,1,e) * dsdz(i,1,1,e))
565 ct(i,1,1,e) = w3(i,1,1) &
566 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
567 + cy(i,1,1,e) * dtdy(i,1,1,e) &
568 + cz(i,1,1,e) * dtdz(i,1,1,e))
569 end do
570 end do
571
572 end subroutine sx_set_convect_rst_lx5
573
574 subroutine sx_set_convect_rst_lx4(cr, cs, ct, cx, cy, cz, &
575 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
576 integer, parameter :: lx = 4
577 integer, intent(in) :: n
578 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
579 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
580 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
581 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
582 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
583 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
584 integer :: e, i
585
586 do i = 1, lx * lx * lx
587 do e = 1, n
588 cr(i,1,1,e) = w3(i,1,1) &
589 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
590 + cy(i,1,1,e) * drdy(i,1,1,e) &
591 + cz(i,1,1,e) * drdz(i,1,1,e) )
592 cs(i,1,1,e) = w3(i,1,1) &
593 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
594 + cy(i,1,1,e) * dsdy(i,1,1,e) &
595 + cz(i,1,1,e) * dsdz(i,1,1,e))
596 ct(i,1,1,e) = w3(i,1,1) &
597 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
598 + cy(i,1,1,e) * dtdy(i,1,1,e) &
599 + cz(i,1,1,e) * dtdz(i,1,1,e))
600 end do
601 end do
602
603 end subroutine sx_set_convect_rst_lx4
604
605 subroutine sx_set_convect_rst_lx3(cr, cs, ct, cx, cy, cz, &
606 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
607 integer, parameter :: lx = 3
608 integer, intent(in) :: n
609 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
610 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
611 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
612 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
613 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
614 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
615 integer :: e, i
616
617 do i = 1, lx * lx * lx
618 do e = 1, n
619 cr(i,1,1,e) = w3(i,1,1) &
620 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
621 + cy(i,1,1,e) * drdy(i,1,1,e) &
622 + cz(i,1,1,e) * drdz(i,1,1,e) )
623 cs(i,1,1,e) = w3(i,1,1) &
624 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
625 + cy(i,1,1,e) * dsdy(i,1,1,e) &
626 + cz(i,1,1,e) * dsdz(i,1,1,e))
627 ct(i,1,1,e) = w3(i,1,1) &
628 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
629 + cy(i,1,1,e) * dtdy(i,1,1,e) &
630 + cz(i,1,1,e) * dtdz(i,1,1,e))
631 end do
632 end do
633
634 end subroutine sx_set_convect_rst_lx3
635
636 subroutine sx_set_convect_rst_lx2(cr, cs, ct, cx, cy, cz, &
637 drdx, dsdx, dtdx, drdy, dsdy, dtdy, drdz, dsdz, dtdz, w3, n)
638 integer, parameter :: lx = 2
639 integer, intent(in) :: n
640 real(kind=rp), dimension(lx, lx, lx, n), intent(inout) :: cr, cs, ct
641 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: cx, cy, cz
642 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdx, dsdx, dtdx
643 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdy, dsdy, dtdy
644 real(kind=rp), dimension(lx, lx, lx, n), intent(in) :: drdz, dsdz, dtdz
645 real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
646 integer :: e, i
647
648 do i = 1, lx * lx * lx
649 do e = 1, n
650 cr(i,1,1,e) = w3(i,1,1) &
651 * ( cx(i,1,1,e) * drdx(i,1,1,e) &
652 + cy(i,1,1,e) * drdy(i,1,1,e) &
653 + cz(i,1,1,e) * drdz(i,1,1,e) )
654 cs(i,1,1,e) = w3(i,1,1) &
655 * ( cx(i,1,1,e) * dsdx(i,1,1,e) &
656 + cy(i,1,1,e) * dsdy(i,1,1,e) &
657 + cz(i,1,1,e) * dsdz(i,1,1,e))
658 ct(i,1,1,e) = w3(i,1,1) &
659 * ( cx(i,1,1,e) * dtdx(i,1,1,e) &
660 + cy(i,1,1,e) * dtdy(i,1,1,e) &
661 + cz(i,1,1,e) * dtdz(i,1,1,e))
662 end do
663 end do
664
665 end subroutine sx_set_convect_rst_lx2
666
667end submodule sx_set_convect_rst
num_types
Definition num_types.f90:1
num_types::rp
integer, parameter, public rp
Global precision used in computations.
Definition num_types.f90:12
opr_sx
Operators SX-Aurora backend.
Definition opr_sx.f90:2