Neko  0.8.99
A portable framework for high-order spectral element flow simulations
cpu_cdtp.f90
Go to the documentation of this file.
1 ! Copyright (c) 2021, 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 !
34 submodule(opr_cpu) cpu_cdtp
35  implicit none
36 
37 contains
38 
39  module subroutine opr_cpu_cdtp(dtx, x, dr, ds, dt, coef, e_start, e_end)
40  type(coef_t), intent(in) :: coef
41  integer, intent(in) :: e_start, e_end
42  real(kind=rp), intent(inout) :: dtx(coef%Xh%lxyz, e_end - e_start + 1)
43  real(kind=rp), intent(inout) :: x(coef%Xh%lxyz, e_end - e_start + 1)
44  real(kind=rp), intent(in) :: dr(coef%Xh%lxyz, e_end - e_start + 1)
45  real(kind=rp), intent(in) :: ds(coef%Xh%lxyz, e_end - e_start + 1)
46  real(kind=rp), intent(in) :: dt(coef%Xh%lxyz, e_end - e_start + 1)
47  integer :: e_len
48  e_len = e_end - e_start + 1
49 
50  if (e_len .eq. 1) then
51  call opr_cpu_cdtp_single(dtx, x, dr, ds, dt, coef, e_start)
52  else
53  call opr_cpu_cdtp_many(dtx, x, dr, ds, dt, coef, e_start, e_len)
54  end if
55 
56  end subroutine opr_cpu_cdtp
57 
58  subroutine opr_cpu_cdtp_many(dtx, x, dr, ds, dt, coef, e_start, e_len)
59  type(coef_t), intent(in) :: coef
60  integer, intent(in) :: e_start, e_len
61  real(kind=rp), intent(inout) :: dtx(coef%Xh%lxyz, e_len)
62  real(kind=rp), intent(inout) :: x(coef%Xh%lxyz, e_len)
63  real(kind=rp), intent(in) :: dr(coef%Xh%lxyz, e_len)
64  real(kind=rp), intent(in) :: ds(coef%Xh%lxyz, e_len)
65  real(kind=rp), intent(in) :: dt(coef%Xh%lxyz, e_len)
66 
67  associate(xh => coef%Xh)
68  select case (xh%lx)
69  case (14)
70  call cpu_cdtp_lx14(dtx, x, &
71  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
72  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
73  case (13)
74  call cpu_cdtp_lx13(dtx, x, &
75  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
76  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
77  case (12)
78  call cpu_cdtp_lx12(dtx, x, &
79  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
80  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
81  case (11)
82  call cpu_cdtp_lx11(dtx, x, &
83  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
84  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
85  case (10)
86  call cpu_cdtp_lx10(dtx, x, &
87  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
88  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
89  case (9)
90  call cpu_cdtp_lx9(dtx, x, &
91  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
92  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
93  case (8)
94  call cpu_cdtp_lx8(dtx, x, &
95  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
96  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
97  case (7)
98  call cpu_cdtp_lx7(dtx, x, &
99  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
100  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
101  case (6)
102  call cpu_cdtp_lx6(dtx, x, &
103  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
104  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
105  case (5)
106  call cpu_cdtp_lx5(dtx, x, &
107  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
108  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
109  case (4)
110  call cpu_cdtp_lx4(dtx, x, &
111  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
112  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
113  case (3)
114  call cpu_cdtp_lx3(dtx, x, &
115  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
116  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
117  case (2)
118  call cpu_cdtp_lx2(dtx, x, &
119  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
120  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len)
121  case default
122  call cpu_cdtp_lx(dtx, x, &
123  dr(1, e_start), ds(1, e_start), dt(1, e_start), &
124  xh%dxt, xh%dyt, xh%dzt, xh%w3, e_len, xh%lx)
125  end select
126  end associate
127 
128  end subroutine opr_cpu_cdtp_many
129 
130  subroutine opr_cpu_cdtp_single(dtx, x, dr, ds, dt, coef, e)
131  integer, parameter :: e_len = 1
132  type(coef_t), intent(in) :: coef
133  integer, intent(in) :: e
134  real(kind=rp), intent(inout) :: dtx(coef%Xh%lxyz, e_len)
135  real(kind=rp), intent(inout) :: x(coef%Xh%lxyz, e_len)
136  real(kind=rp), intent(in) :: dr(coef%Xh%lxyz, e_len)
137  real(kind=rp), intent(in) :: ds(coef%Xh%lxyz, e_len)
138  real(kind=rp), intent(in) :: dt(coef%Xh%lxyz, e_len)
139 
140  associate(xh => coef%Xh)
141  select case (xh%lx)
142  case (14)
143  call cpu_cdtp_lx14_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
144  xh%dxt, xh%dyt, xh%dzt, xh%w3)
145  case (13)
146  call cpu_cdtp_lx13_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
147  xh%dxt, xh%dyt, xh%dzt, xh%w3)
148  case (12)
149  call cpu_cdtp_lx12_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
150  xh%dxt, xh%dyt, xh%dzt, xh%w3)
151  case (11)
152  call cpu_cdtp_lx11_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
153  xh%dxt, xh%dyt, xh%dzt, xh%w3)
154  case (10)
155  call cpu_cdtp_lx10_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
156  xh%dxt, xh%dyt, xh%dzt, xh%w3)
157  case (9)
158  call cpu_cdtp_lx9_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
159  xh%dxt, xh%dyt, xh%dzt, xh%w3)
160  case (8)
161  call cpu_cdtp_lx8_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
162  xh%dxt, xh%dyt, xh%dzt, xh%w3)
163  case (7)
164  call cpu_cdtp_lx7_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
165  xh%dxt, xh%dyt, xh%dzt, xh%w3)
166  case (6)
167  call cpu_cdtp_lx6_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
168  xh%dxt, xh%dyt, xh%dzt, xh%w3)
169  case (5)
170  call cpu_cdtp_lx5_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
171  xh%dxt, xh%dyt, xh%dzt, xh%w3)
172  case (4)
173  call cpu_cdtp_lx4_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
174  xh%dxt, xh%dyt, xh%dzt, xh%w3)
175  case (3)
176  call cpu_cdtp_lx3_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
177  xh%dxt, xh%dyt, xh%dzt, xh%w3)
178  case (2)
179  call cpu_cdtp_lx2_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
180  xh%dxt, xh%dyt, xh%dzt, xh%w3)
181  case default
182  call cpu_cdtp_lx_single(dtx, x, dr(1,e), ds(1,e), dt(1,e), &
183  xh%dxt, xh%dyt, xh%dzt, xh%w3, xh%lx)
184  end select
185  end associate
186 
187  end subroutine opr_cpu_cdtp_single
188 
189  subroutine cpu_cdtp_lx(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel, lx)
190  integer, intent(in) :: nel, lx
191  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
192  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
193  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
194  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
195  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
196  real(kind=rp) :: tmp
197  integer :: e, i, j, k, l
198 
199  do e = 1, nel
200 
201  do i = 1, lx*lx*lx
202  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
203  end do
204 
205  do i = 1, lx*lx*lx
206  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
207  end do
208 
209  do j = 1, lx * lx
210  do i = 1, lx
211  tmp = 0.0_rp
212  !DIR$ LOOP_INFO MIN_TRIPS(15)
213  do k = 1, lx
214  tmp = tmp + dxt(i,k) * ta1(k,j,1)
215  end do
216  dtx(i,j,1,e) = tmp
217  end do
218  end do
219 
220  do i = 1, lx*lx*lx
221  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
222  end do
223 
224  do k = 1, lx
225  do j = 1, lx
226  do i = 1, lx
227  tmp = 0.0_rp
228  !DIR$ LOOP_INFO MIN_TRIPS(15)
229  do l = 1, lx
230  tmp = tmp + dyt(j,l) * ta1(i,l,k)
231  end do
232  dtx(i,j,k,e) = dtx(i,j,k,e) + tmp
233  end do
234  end do
235  end do
236 
237  do i = 1, lx*lx*lx
238  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
239  end do
240 
241  do k = 1, lx
242  do i = 1, lx*lx
243  tmp = 0.0_rp
244  !DIR$ LOOP_INFO MIN_TRIPS(15)
245  do l = 1, lx
246  tmp = tmp + dzt(k,l) * ta1(i,1,l)
247  end do
248  dtx(i,1,k,e) = dtx(i,1,k,e) + tmp
249  end do
250  end do
251 
252  end do
253  end subroutine cpu_cdtp_lx
254 
255  subroutine cpu_cdtp_lx14(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
256  integer, parameter :: lx = 14
257  integer, intent(in) :: nel
258  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
259  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
260  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
261  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
262  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
263  integer :: e, i, j, k
264 
265  do e = 1, nel
266 
267  do i = 1, lx*lx*lx
268  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
269  end do
270 
271  do i = 1, lx*lx*lx
272  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
273  end do
274 
275  do j = 1, lx * lx
276  do i = 1, lx
277  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
278  + dxt(i,2) * ta1(2,j,1) &
279  + dxt(i,3) * ta1(3,j,1) &
280  + dxt(i,4) * ta1(4,j,1) &
281  + dxt(i,5) * ta1(5,j,1) &
282  + dxt(i,6) * ta1(6,j,1) &
283  + dxt(i,7) * ta1(7,j,1) &
284  + dxt(i,8) * ta1(8,j,1) &
285  + dxt(i,9) * ta1(9,j,1) &
286  + dxt(i,10) * ta1(10,j,1) &
287  + dxt(i,11) * ta1(11,j,1) &
288  + dxt(i,12) * ta1(12,j,1) &
289  + dxt(i,13) * ta1(13,j,1) &
290  + dxt(i,14) * ta1(14,j,1)
291  end do
292  end do
293 
294  do i = 1, lx*lx*lx
295  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
296  end do
297 
298  do k = 1, lx
299  do j = 1, lx
300  do i = 1, lx
301  dtx(i,j,k,e) = dtx(i,j,k,e) &
302  + dyt(j,1) * ta1(i,1,k) &
303  + dyt(j,2) * ta1(i,2,k) &
304  + dyt(j,3) * ta1(i,3,k) &
305  + dyt(j,4) * ta1(i,4,k) &
306  + dyt(j,5) * ta1(i,5,k) &
307  + dyt(j,6) * ta1(i,6,k) &
308  + dyt(j,7) * ta1(i,7,k) &
309  + dyt(j,8) * ta1(i,8,k) &
310  + dyt(j,9) * ta1(i,9,k) &
311  + dyt(j,10) * ta1(i,10,k) &
312  + dyt(j,11) * ta1(i,11,k) &
313  + dyt(j,12) * ta1(i,12,k) &
314  + dyt(j,13) * ta1(i,13,k) &
315  + dyt(j,14) * ta1(i,14,k)
316  end do
317  end do
318  end do
319 
320  do i = 1, lx*lx*lx
321  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
322  end do
323 
324  do k = 1, lx
325  do i = 1, lx*lx
326  dtx(i,1,k,e) = dtx(i,1,k,e) &
327  + dzt(k,1) * ta1(i,1,1) &
328  + dzt(k,2) * ta1(i,1,2) &
329  + dzt(k,3) * ta1(i,1,3) &
330  + dzt(k,4) * ta1(i,1,4) &
331  + dzt(k,5) * ta1(i,1,5) &
332  + dzt(k,6) * ta1(i,1,6) &
333  + dzt(k,7) * ta1(i,1,7) &
334  + dzt(k,8) * ta1(i,1,8) &
335  + dzt(k,9) * ta1(i,1,9) &
336  + dzt(k,10) * ta1(i,1,10) &
337  + dzt(k,11) * ta1(i,1,11) &
338  + dzt(k,12) * ta1(i,1,12) &
339  + dzt(k,13) * ta1(i,1,13) &
340  + dzt(k,14) * ta1(i,1,14)
341  end do
342  end do
343 
344  end do
345  end subroutine cpu_cdtp_lx14
346 
347  subroutine cpu_cdtp_lx13(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
348  integer, parameter :: lx = 13
349  integer, intent(in) :: nel
350  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
351  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
352  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
353  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
354  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
355  integer :: e, i, j, k
356 
357  do e = 1, nel
358 
359  do i = 1, lx*lx*lx
360  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
361  end do
362 
363  do i = 1, lx*lx*lx
364  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
365  end do
366 
367  do j = 1, lx * lx
368  do i = 1, lx
369  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
370  + dxt(i,2) * ta1(2,j,1) &
371  + dxt(i,3) * ta1(3,j,1) &
372  + dxt(i,4) * ta1(4,j,1) &
373  + dxt(i,5) * ta1(5,j,1) &
374  + dxt(i,6) * ta1(6,j,1) &
375  + dxt(i,7) * ta1(7,j,1) &
376  + dxt(i,8) * ta1(8,j,1) &
377  + dxt(i,9) * ta1(9,j,1) &
378  + dxt(i,10) * ta1(10,j,1) &
379  + dxt(i,11) * ta1(11,j,1) &
380  + dxt(i,12) * ta1(12,j,1) &
381  + dxt(i,13) * ta1(13,j,1)
382  end do
383  end do
384 
385  do i = 1, lx*lx*lx
386  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
387  end do
388 
389  do k = 1, lx
390  do j = 1, lx
391  do i = 1, lx
392  dtx(i,j,k,e) = dtx(i,j,k,e) &
393  + dyt(j,1) * ta1(i,1,k) &
394  + dyt(j,2) * ta1(i,2,k) &
395  + dyt(j,3) * ta1(i,3,k) &
396  + dyt(j,4) * ta1(i,4,k) &
397  + dyt(j,5) * ta1(i,5,k) &
398  + dyt(j,6) * ta1(i,6,k) &
399  + dyt(j,7) * ta1(i,7,k) &
400  + dyt(j,8) * ta1(i,8,k) &
401  + dyt(j,9) * ta1(i,9,k) &
402  + dyt(j,10) * ta1(i,10,k) &
403  + dyt(j,11) * ta1(i,11,k) &
404  + dyt(j,12) * ta1(i,12,k) &
405  + dyt(j,13) * ta1(i,13,k)
406  end do
407  end do
408  end do
409 
410  do i = 1, lx*lx*lx
411  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
412  end do
413 
414  do k = 1, lx
415  do i = 1, lx*lx
416  dtx(i,1,k,e) = dtx(i,1,k,e) &
417  + dzt(k,1) * ta1(i,1,1) &
418  + dzt(k,2) * ta1(i,1,2) &
419  + dzt(k,3) * ta1(i,1,3) &
420  + dzt(k,4) * ta1(i,1,4) &
421  + dzt(k,5) * ta1(i,1,5) &
422  + dzt(k,6) * ta1(i,1,6) &
423  + dzt(k,7) * ta1(i,1,7) &
424  + dzt(k,8) * ta1(i,1,8) &
425  + dzt(k,9) * ta1(i,1,9) &
426  + dzt(k,10) * ta1(i,1,10) &
427  + dzt(k,11) * ta1(i,1,11) &
428  + dzt(k,12) * ta1(i,1,12) &
429  + dzt(k,13) * ta1(i,1,13)
430  end do
431  end do
432 
433  end do
434  end subroutine cpu_cdtp_lx13
435 
436  subroutine cpu_cdtp_lx12(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
437  integer, parameter :: lx = 12
438  integer, intent(in) :: nel
439  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
440  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
441  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
442  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
443  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
444  integer :: e, i, j, k
445 
446  do e = 1, nel
447 
448  do i = 1, lx*lx*lx
449  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
450  end do
451 
452  do i = 1, lx*lx*lx
453  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
454  end do
455 
456  do j = 1, lx * lx
457  do i = 1, lx
458  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
459  + dxt(i,2) * ta1(2,j,1) &
460  + dxt(i,3) * ta1(3,j,1) &
461  + dxt(i,4) * ta1(4,j,1) &
462  + dxt(i,5) * ta1(5,j,1) &
463  + dxt(i,6) * ta1(6,j,1) &
464  + dxt(i,7) * ta1(7,j,1) &
465  + dxt(i,8) * ta1(8,j,1) &
466  + dxt(i,9) * ta1(9,j,1) &
467  + dxt(i,10) * ta1(10,j,1) &
468  + dxt(i,11) * ta1(11,j,1) &
469  + dxt(i,12) * ta1(12,j,1)
470  end do
471  end do
472 
473  do i = 1, lx*lx*lx
474  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
475  end do
476 
477  do k = 1, lx
478  do j = 1, lx
479  do i = 1, lx
480  dtx(i,j,k,e) = dtx(i,j,k,e) &
481  + dyt(j,1) * ta1(i,1,k) &
482  + dyt(j,2) * ta1(i,2,k) &
483  + dyt(j,3) * ta1(i,3,k) &
484  + dyt(j,4) * ta1(i,4,k) &
485  + dyt(j,5) * ta1(i,5,k) &
486  + dyt(j,6) * ta1(i,6,k) &
487  + dyt(j,7) * ta1(i,7,k) &
488  + dyt(j,8) * ta1(i,8,k) &
489  + dyt(j,9) * ta1(i,9,k) &
490  + dyt(j,10) * ta1(i,10,k) &
491  + dyt(j,11) * ta1(i,11,k) &
492  + dyt(j,12) * ta1(i,12,k)
493  end do
494  end do
495  end do
496 
497  do i = 1, lx*lx*lx
498  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
499  end do
500 
501  do k = 1, lx
502  do i = 1, lx*lx
503  dtx(i,1,k,e) = dtx(i,1,k,e) &
504  + dzt(k,1) * ta1(i,1,1) &
505  + dzt(k,2) * ta1(i,1,2) &
506  + dzt(k,3) * ta1(i,1,3) &
507  + dzt(k,4) * ta1(i,1,4) &
508  + dzt(k,5) * ta1(i,1,5) &
509  + dzt(k,6) * ta1(i,1,6) &
510  + dzt(k,7) * ta1(i,1,7) &
511  + dzt(k,8) * ta1(i,1,8) &
512  + dzt(k,9) * ta1(i,1,9) &
513  + dzt(k,10) * ta1(i,1,10) &
514  + dzt(k,11) * ta1(i,1,11) &
515  + dzt(k,12) * ta1(i,1,12)
516  end do
517  end do
518 
519  end do
520  end subroutine cpu_cdtp_lx12
521 
522  subroutine cpu_cdtp_lx11(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
523  integer, parameter :: lx = 11
524  integer, intent(in) :: nel
525  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
526  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
527  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
528  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
529  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
530  integer :: e, i, j, k
531 
532  do e = 1, nel
533 
534  do i = 1, lx*lx*lx
535  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
536  end do
537 
538  do i = 1, lx*lx*lx
539  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
540  end do
541 
542  do j = 1, lx * lx
543  do i = 1, lx
544  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
545  + dxt(i,2) * ta1(2,j,1) &
546  + dxt(i,3) * ta1(3,j,1) &
547  + dxt(i,4) * ta1(4,j,1) &
548  + dxt(i,5) * ta1(5,j,1) &
549  + dxt(i,6) * ta1(6,j,1) &
550  + dxt(i,7) * ta1(7,j,1) &
551  + dxt(i,8) * ta1(8,j,1) &
552  + dxt(i,9) * ta1(9,j,1) &
553  + dxt(i,10) * ta1(10,j,1) &
554  + dxt(i,11) * ta1(11,j,1)
555  end do
556  end do
557 
558  do i = 1, lx*lx*lx
559  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
560  end do
561 
562  do k = 1, lx
563  do j = 1, lx
564  do i = 1, lx
565  dtx(i,j,k,e) = dtx(i,j,k,e) &
566  + dyt(j,1) * ta1(i,1,k) &
567  + dyt(j,2) * ta1(i,2,k) &
568  + dyt(j,3) * ta1(i,3,k) &
569  + dyt(j,4) * ta1(i,4,k) &
570  + dyt(j,5) * ta1(i,5,k) &
571  + dyt(j,6) * ta1(i,6,k) &
572  + dyt(j,7) * ta1(i,7,k) &
573  + dyt(j,8) * ta1(i,8,k) &
574  + dyt(j,9) * ta1(i,9,k) &
575  + dyt(j,10) * ta1(i,10,k) &
576  + dyt(j,11) * ta1(i,11,k)
577  end do
578  end do
579  end do
580 
581  do i = 1, lx*lx*lx
582  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
583  end do
584 
585  do k = 1, lx
586  do i = 1, lx*lx
587  dtx(i,1,k,e) = dtx(i,1,k,e) &
588  + dzt(k,1) * ta1(i,1,1) &
589  + dzt(k,2) * ta1(i,1,2) &
590  + dzt(k,3) * ta1(i,1,3) &
591  + dzt(k,4) * ta1(i,1,4) &
592  + dzt(k,5) * ta1(i,1,5) &
593  + dzt(k,6) * ta1(i,1,6) &
594  + dzt(k,7) * ta1(i,1,7) &
595  + dzt(k,8) * ta1(i,1,8) &
596  + dzt(k,9) * ta1(i,1,9) &
597  + dzt(k,10) * ta1(i,1,10) &
598  + dzt(k,11) * ta1(i,1,11)
599  end do
600  end do
601 
602  end do
603  end subroutine cpu_cdtp_lx11
604 
605  subroutine cpu_cdtp_lx10(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
606  integer, parameter :: lx = 10
607  integer, intent(in) :: nel
608  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
609  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
610  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
611  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
612  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
613  integer :: e, i, j, k
614 
615  do e = 1, nel
616 
617  do i = 1, lx*lx*lx
618  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
619  end do
620 
621  do i = 1, lx*lx*lx
622  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
623  end do
624 
625  do j = 1, lx * lx
626  do i = 1, lx
627  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
628  + dxt(i,2) * ta1(2,j,1) &
629  + dxt(i,3) * ta1(3,j,1) &
630  + dxt(i,4) * ta1(4,j,1) &
631  + dxt(i,5) * ta1(5,j,1) &
632  + dxt(i,6) * ta1(6,j,1) &
633  + dxt(i,7) * ta1(7,j,1) &
634  + dxt(i,8) * ta1(8,j,1) &
635  + dxt(i,9) * ta1(9,j,1) &
636  + dxt(i,10) * ta1(10,j,1)
637  end do
638  end do
639 
640  do i = 1, lx*lx*lx
641  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
642  end do
643 
644  do k = 1, lx
645  do j = 1, lx
646  do i = 1, lx
647  dtx(i,j,k,e) = dtx(i,j,k,e) &
648  + dyt(j,1) * ta1(i,1,k) &
649  + dyt(j,2) * ta1(i,2,k) &
650  + dyt(j,3) * ta1(i,3,k) &
651  + dyt(j,4) * ta1(i,4,k) &
652  + dyt(j,5) * ta1(i,5,k) &
653  + dyt(j,6) * ta1(i,6,k) &
654  + dyt(j,7) * ta1(i,7,k) &
655  + dyt(j,8) * ta1(i,8,k) &
656  + dyt(j,9) * ta1(i,9,k) &
657  + dyt(j,10) * ta1(i,10,k)
658 
659  end do
660  end do
661  end do
662 
663  do i = 1, lx*lx*lx
664  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
665  end do
666 
667  do k = 1, lx
668  do i = 1, lx*lx
669  dtx(i,1,k,e) = dtx(i,1,k,e) &
670  + dzt(k,1) * ta1(i,1,1) &
671  + dzt(k,2) * ta1(i,1,2) &
672  + dzt(k,3) * ta1(i,1,3) &
673  + dzt(k,4) * ta1(i,1,4) &
674  + dzt(k,5) * ta1(i,1,5) &
675  + dzt(k,6) * ta1(i,1,6) &
676  + dzt(k,7) * ta1(i,1,7) &
677  + dzt(k,8) * ta1(i,1,8) &
678  + dzt(k,9) * ta1(i,1,9) &
679  + dzt(k,10) * ta1(i,1,10)
680  end do
681  end do
682 
683  end do
684  end subroutine cpu_cdtp_lx10
685 
686  subroutine cpu_cdtp_lx9(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
687  integer, parameter :: lx = 9
688  integer, intent(in) :: nel
689  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
690  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
691  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
692  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
693  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
694  integer :: e, i, j, k
695 
696  do e = 1, nel
697 
698  do i = 1, lx*lx*lx
699  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
700  end do
701 
702  do i = 1, lx*lx*lx
703  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
704  end do
705 
706  do j = 1, lx * lx
707  do i = 1, lx
708  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
709  + dxt(i,2) * ta1(2,j,1) &
710  + dxt(i,3) * ta1(3,j,1) &
711  + dxt(i,4) * ta1(4,j,1) &
712  + dxt(i,5) * ta1(5,j,1) &
713  + dxt(i,6) * ta1(6,j,1) &
714  + dxt(i,7) * ta1(7,j,1) &
715  + dxt(i,8) * ta1(8,j,1) &
716  + dxt(i,9) * ta1(9,j,1)
717  end do
718  end do
719 
720  do i = 1, lx*lx*lx
721  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
722  end do
723 
724  do k = 1, lx
725  do j = 1, lx
726  do i = 1, lx
727  dtx(i,j,k,e) = dtx(i,j,k,e) &
728  + dyt(j,1) * ta1(i,1,k) &
729  + dyt(j,2) * ta1(i,2,k) &
730  + dyt(j,3) * ta1(i,3,k) &
731  + dyt(j,4) * ta1(i,4,k) &
732  + dyt(j,5) * ta1(i,5,k) &
733  + dyt(j,6) * ta1(i,6,k) &
734  + dyt(j,7) * ta1(i,7,k) &
735  + dyt(j,8) * ta1(i,8,k) &
736  + dyt(j,9) * ta1(i,9,k)
737  end do
738  end do
739  end do
740 
741  do i = 1, lx*lx*lx
742  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
743  end do
744 
745  do k = 1, lx
746  do i = 1, lx*lx
747  dtx(i,1,k,e) = dtx(i,1,k,e) &
748  + dzt(k,1) * ta1(i,1,1) &
749  + dzt(k,2) * ta1(i,1,2) &
750  + dzt(k,3) * ta1(i,1,3) &
751  + dzt(k,4) * ta1(i,1,4) &
752  + dzt(k,5) * ta1(i,1,5) &
753  + dzt(k,6) * ta1(i,1,6) &
754  + dzt(k,7) * ta1(i,1,7) &
755  + dzt(k,8) * ta1(i,1,8) &
756  + dzt(k,9) * ta1(i,1,9)
757  end do
758  end do
759 
760  end do
761  end subroutine cpu_cdtp_lx9
762 
763  subroutine cpu_cdtp_lx8(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
764  integer, parameter :: lx = 8
765  integer, intent(in) :: nel
766  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
767  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
768  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
769  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
770  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
771  integer :: e, i, j, k
772 
773  do e = 1, nel
774 
775  do i = 1, lx*lx*lx
776  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
777  end do
778 
779  do i = 1, lx*lx*lx
780  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
781  end do
782 
783  do j = 1, lx * lx
784  do i = 1, lx
785  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
786  + dxt(i,2) * ta1(2,j,1) &
787  + dxt(i,3) * ta1(3,j,1) &
788  + dxt(i,4) * ta1(4,j,1) &
789  + dxt(i,5) * ta1(5,j,1) &
790  + dxt(i,6) * ta1(6,j,1) &
791  + dxt(i,7) * ta1(7,j,1) &
792  + dxt(i,8) * ta1(8,j,1)
793  end do
794  end do
795 
796  do i = 1, lx*lx*lx
797  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
798  end do
799 
800  do k = 1, lx
801  do j = 1, lx
802  do i = 1, lx
803  dtx(i,j,k,e) = dtx(i,j,k,e) &
804  + dyt(j,1) * ta1(i,1,k) &
805  + dyt(j,2) * ta1(i,2,k) &
806  + dyt(j,3) * ta1(i,3,k) &
807  + dyt(j,4) * ta1(i,4,k) &
808  + dyt(j,5) * ta1(i,5,k) &
809  + dyt(j,6) * ta1(i,6,k) &
810  + dyt(j,7) * ta1(i,7,k) &
811  + dyt(j,8) * ta1(i,8,k)
812  end do
813  end do
814  end do
815 
816  do i = 1, lx*lx*lx
817  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
818  end do
819 
820  do k = 1, lx
821  do i = 1, lx*lx
822  dtx(i,1,k,e) = dtx(i,1,k,e) &
823  + dzt(k,1) * ta1(i,1,1) &
824  + dzt(k,2) * ta1(i,1,2) &
825  + dzt(k,3) * ta1(i,1,3) &
826  + dzt(k,4) * ta1(i,1,4) &
827  + dzt(k,5) * ta1(i,1,5) &
828  + dzt(k,6) * ta1(i,1,6) &
829  + dzt(k,7) * ta1(i,1,7) &
830  + dzt(k,8) * ta1(i,1,8)
831  end do
832  end do
833 
834  end do
835  end subroutine cpu_cdtp_lx8
836 
837  subroutine cpu_cdtp_lx7(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
838  integer, parameter :: lx = 7
839  integer, intent(in) :: nel
840  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
841  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
842  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
843  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
844  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
845  integer :: e, i, j, k
846 
847  do e = 1, nel
848 
849  do i = 1, lx*lx*lx
850  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
851  end do
852 
853  do i = 1, lx*lx*lx
854  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
855  end do
856 
857  do j = 1, lx * lx
858  do i = 1, lx
859  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
860  + dxt(i,2) * ta1(2,j,1) &
861  + dxt(i,3) * ta1(3,j,1) &
862  + dxt(i,4) * ta1(4,j,1) &
863  + dxt(i,5) * ta1(5,j,1) &
864  + dxt(i,6) * ta1(6,j,1) &
865  + dxt(i,7) * ta1(7,j,1)
866  end do
867  end do
868 
869  do i = 1, lx*lx*lx
870  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
871  end do
872 
873  do k = 1, lx
874  do j = 1, lx
875  do i = 1, lx
876  dtx(i,j,k,e) = dtx(i,j,k,e) &
877  + dyt(j,1) * ta1(i,1,k) &
878  + dyt(j,2) * ta1(i,2,k) &
879  + dyt(j,3) * ta1(i,3,k) &
880  + dyt(j,4) * ta1(i,4,k) &
881  + dyt(j,5) * ta1(i,5,k) &
882  + dyt(j,6) * ta1(i,6,k) &
883  + dyt(j,7) * ta1(i,7,k)
884  end do
885  end do
886  end do
887 
888  do i = 1, lx*lx*lx
889  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
890  end do
891 
892  do k = 1, lx
893  do i = 1, lx*lx
894  dtx(i,1,k,e) = dtx(i,1,k,e) &
895  + dzt(k,1) * ta1(i,1,1) &
896  + dzt(k,2) * ta1(i,1,2) &
897  + dzt(k,3) * ta1(i,1,3) &
898  + dzt(k,4) * ta1(i,1,4) &
899  + dzt(k,5) * ta1(i,1,5) &
900  + dzt(k,6) * ta1(i,1,6) &
901  + dzt(k,7) * ta1(i,1,7)
902  end do
903  end do
904 
905  end do
906  end subroutine cpu_cdtp_lx7
907 
908  subroutine cpu_cdtp_lx6(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
909  integer, parameter :: lx = 6
910  integer, intent(in) :: nel
911  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
912  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
913  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
914  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
915  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
916  integer :: e, i, j, k
917 
918  do e = 1, nel
919 
920  do i = 1, lx*lx*lx
921  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
922  end do
923 
924  do i = 1, lx*lx*lx
925  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
926  end do
927 
928  do j = 1, lx * lx
929  do i = 1, lx
930  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
931  + dxt(i,2) * ta1(2,j,1) &
932  + dxt(i,3) * ta1(3,j,1) &
933  + dxt(i,4) * ta1(4,j,1) &
934  + dxt(i,5) * ta1(5,j,1) &
935  + dxt(i,6) * ta1(6,j,1)
936  end do
937  end do
938 
939  do i = 1, lx*lx*lx
940  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
941  end do
942 
943  do k = 1, lx
944  do j = 1, lx
945  do i = 1, lx
946  dtx(i,j,k,e) = dtx(i,j,k,e) &
947  + dyt(j,1) * ta1(i,1,k) &
948  + dyt(j,2) * ta1(i,2,k) &
949  + dyt(j,3) * ta1(i,3,k) &
950  + dyt(j,4) * ta1(i,4,k) &
951  + dyt(j,5) * ta1(i,5,k) &
952  + dyt(j,6) * ta1(i,6,k)
953  end do
954  end do
955  end do
956 
957  do i = 1, lx*lx*lx
958  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
959  end do
960 
961  do k = 1, lx
962  do i = 1, lx*lx
963  dtx(i,1,k,e) = dtx(i,1,k,e) &
964  + dzt(k,1) * ta1(i,1,1) &
965  + dzt(k,2) * ta1(i,1,2) &
966  + dzt(k,3) * ta1(i,1,3) &
967  + dzt(k,4) * ta1(i,1,4) &
968  + dzt(k,5) * ta1(i,1,5) &
969  + dzt(k,6) * ta1(i,1,6)
970  end do
971  end do
972 
973  end do
974  end subroutine cpu_cdtp_lx6
975 
976  subroutine cpu_cdtp_lx5(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
977  integer, parameter :: lx = 5
978  integer, intent(in) :: nel
979  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
980  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
981  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
982  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
983  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
984  integer :: e, i, j, k
985 
986  do e = 1, nel
987 
988  do i = 1, lx*lx*lx
989  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
990  end do
991 
992  do i = 1, lx*lx*lx
993  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
994  end do
995 
996  do j = 1, lx * lx
997  do i = 1, lx
998  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
999  + dxt(i,2) * ta1(2,j,1) &
1000  + dxt(i,3) * ta1(3,j,1) &
1001  + dxt(i,4) * ta1(4,j,1) &
1002  + dxt(i,5) * ta1(5,j,1)
1003  end do
1004  end do
1005 
1006  do i = 1, lx*lx*lx
1007  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
1008  end do
1009 
1010  do k = 1, lx
1011  do j = 1, lx
1012  do i = 1, lx
1013  dtx(i,j,k,e) = dtx(i,j,k,e) &
1014  + dyt(j,1) * ta1(i,1,k) &
1015  + dyt(j,2) * ta1(i,2,k) &
1016  + dyt(j,3) * ta1(i,3,k) &
1017  + dyt(j,4) * ta1(i,4,k) &
1018  + dyt(j,5) * ta1(i,5,k)
1019  end do
1020  end do
1021  end do
1022 
1023  do i = 1, lx*lx*lx
1024  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
1025  end do
1026 
1027  do k = 1, lx
1028  do i = 1, lx*lx
1029  dtx(i,1,k,e) = dtx(i,1,k,e) &
1030  + dzt(k,1) * ta1(i,1,1) &
1031  + dzt(k,2) * ta1(i,1,2) &
1032  + dzt(k,3) * ta1(i,1,3) &
1033  + dzt(k,4) * ta1(i,1,4) &
1034  + dzt(k,5) * ta1(i,1,5)
1035  end do
1036  end do
1037 
1038  end do
1039  end subroutine cpu_cdtp_lx5
1040 
1041  subroutine cpu_cdtp_lx4(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
1042  integer, parameter :: lx = 4
1043  integer, intent(in) :: nel
1044  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
1045  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
1046  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1047  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1048  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1049  integer :: e, i, j, k
1050 
1051  do e = 1, nel
1052 
1053  do i = 1, lx*lx*lx
1054  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
1055  end do
1056 
1057  do i = 1, lx*lx*lx
1058  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
1059  end do
1060 
1061  do j = 1, lx * lx
1062  do i = 1, lx
1063  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
1064  + dxt(i,2) * ta1(2,j,1) &
1065  + dxt(i,3) * ta1(3,j,1) &
1066  + dxt(i,4) * ta1(4,j,1)
1067  end do
1068  end do
1069 
1070  do i = 1, lx*lx*lx
1071  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
1072  end do
1073 
1074  do k = 1, lx
1075  do j = 1, lx
1076  do i = 1, lx
1077  dtx(i,j,k,e) = dtx(i,j,k,e) &
1078  + dyt(j,1) * ta1(i,1,k) &
1079  + dyt(j,2) * ta1(i,2,k) &
1080  + dyt(j,3) * ta1(i,3,k) &
1081  + dyt(j,4) * ta1(i,4,k)
1082  end do
1083  end do
1084  end do
1085 
1086  do i = 1, lx*lx*lx
1087  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
1088  end do
1089 
1090  do k = 1, lx
1091  do i = 1, lx*lx
1092  dtx(i,1,k,e) = dtx(i,1,k,e) &
1093  + dzt(k,1) * ta1(i,1,1) &
1094  + dzt(k,2) * ta1(i,1,2) &
1095  + dzt(k,3) * ta1(i,1,3) &
1096  + dzt(k,4) * ta1(i,1,4)
1097  end do
1098  end do
1099 
1100  end do
1101  end subroutine cpu_cdtp_lx4
1102 
1103  subroutine cpu_cdtp_lx3(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
1104  integer, parameter :: lx = 3
1105  integer, intent(in) :: nel
1106  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
1107  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
1108  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1109  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1110  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1111  integer :: e, i, j, k
1112 
1113  do e = 1, nel
1114 
1115  do i = 1, lx*lx*lx
1116  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
1117  end do
1118 
1119  do i = 1, lx*lx*lx
1120  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
1121  end do
1122 
1123  do j = 1, lx * lx
1124  do i = 1, lx
1125  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
1126  + dxt(i,2) * ta1(2,j,1) &
1127  + dxt(i,3) * ta1(3,j,1)
1128  end do
1129  end do
1130 
1131  do i = 1, lx*lx*lx
1132  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
1133  end do
1134 
1135  do k = 1, lx
1136  do j = 1, lx
1137  do i = 1, lx
1138  dtx(i,j,k,e) = dtx(i,j,k,e) &
1139  + dyt(j,1) * ta1(i,1,k) &
1140  + dyt(j,2) * ta1(i,2,k) &
1141  + dyt(j,3) * ta1(i,3,k)
1142  end do
1143  end do
1144  end do
1145 
1146  do i = 1, lx*lx*lx
1147  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
1148  end do
1149 
1150  do k = 1, lx
1151  do i = 1, lx*lx
1152  dtx(i,1,k,e) = dtx(i,1,k,e) &
1153  + dzt(k,1) * ta1(i,1,1) &
1154  + dzt(k,2) * ta1(i,1,2) &
1155  + dzt(k,3) * ta1(i,1,3)
1156  end do
1157  end do
1158 
1159  end do
1160  end subroutine cpu_cdtp_lx3
1161 
1162  subroutine cpu_cdtp_lx2(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, nel)
1163  integer, parameter :: lx = 2
1164  integer, intent(in) :: nel
1165  real(kind=rp), dimension(lx, lx, lx, nel), intent(inout) :: dtx
1166  real(kind=rp), dimension(lx, lx, lx, nel), intent(in) :: x, dr, ds, dt
1167  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1168  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1169  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1170  integer :: e, i, j, k
1171 
1172  do e = 1, nel
1173 
1174  do i = 1, lx*lx*lx
1175  wx(i,1,1) = x(i,1,1,e) * w3(i,1,1)
1176  end do
1177 
1178  do i = 1, lx*lx*lx
1179  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1,e)
1180  end do
1181 
1182  do j = 1, lx * lx
1183  do i = 1, lx
1184  dtx(i,j,1,e) = dxt(i,1) * ta1(1,j,1) &
1185  + dxt(i,2) * ta1(2,j,1)
1186  end do
1187  end do
1188 
1189  do i = 1, lx*lx*lx
1190  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1,e)
1191  end do
1192 
1193  do k = 1, lx
1194  do j = 1, lx
1195  do i = 1, lx
1196  dtx(i,j,k,e) = dtx(i,j,k,e) &
1197  + dyt(j,1) * ta1(i,1,k) &
1198  + dyt(j,2) * ta1(i,2,k)
1199  end do
1200  end do
1201  end do
1202 
1203  do i = 1, lx*lx*lx
1204  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1,e)
1205  end do
1206 
1207  do k = 1, lx
1208  do i = 1, lx*lx
1209  dtx(i,1,k,e) = dtx(i,1,k,e) &
1210  + dzt(k,1) * ta1(i,1,1) &
1211  + dzt(k,2) * ta1(i,1,2)
1212  end do
1213  end do
1214 
1215  end do
1216  end subroutine cpu_cdtp_lx2
1217 
1218  subroutine cpu_cdtp_lx_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3, lx)
1219  integer, intent(in) :: lx
1220  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1221  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1222  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1223  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1224  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1225  real(kind=rp) :: tmp
1226  integer :: i, j, k, l
1227 
1228  do i = 1, lx*lx*lx
1229  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1230  end do
1231 
1232  do i = 1, lx*lx*lx
1233  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1234  end do
1235 
1236  do j = 1, lx * lx
1237  do i = 1, lx
1238  tmp = 0.0_rp
1239  !DIR$ LOOP_INFO MIN_TRIPS(15)
1240  do k = 1, lx
1241  tmp = tmp + dxt(i,k) * ta1(k,j,1)
1242  end do
1243  dtx(i,j,1) = tmp
1244  end do
1245  end do
1246 
1247  do i = 1, lx*lx*lx
1248  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1249  end do
1250 
1251  do k = 1, lx
1252  do j = 1, lx
1253  do i = 1, lx
1254  tmp = 0.0_rp
1255  !DIR$ LOOP_INFO MIN_TRIPS(15)
1256  do l = 1, lx
1257  tmp = tmp + dyt(j,l) * ta1(i,l,k)
1258  end do
1259  dtx(i,j,k) = dtx(i,j,k) + tmp
1260  end do
1261  end do
1262  end do
1263 
1264  do i = 1, lx*lx*lx
1265  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1266  end do
1267 
1268  do k = 1, lx
1269  do i = 1, lx*lx
1270  tmp = 0.0_rp
1271  !DIR$ LOOP_INFO MIN_TRIPS(15)
1272  do l = 1, lx
1273  tmp = tmp + dzt(k,l) * ta1(i,1,l)
1274  end do
1275  dtx(i,1,k) = dtx(i,1,k) + tmp
1276  end do
1277  end do
1278 
1279  end subroutine cpu_cdtp_lx_single
1280 
1281  subroutine cpu_cdtp_lx14_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1282  integer, parameter :: lx = 14
1283  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1284  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1285  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1286  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1287  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1288  integer :: i, j, k
1289 
1290  do i = 1, lx*lx*lx
1291  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1292  end do
1293 
1294  do i = 1, lx*lx*lx
1295  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1296  end do
1297 
1298  do j = 1, lx * lx
1299  do i = 1, lx
1300  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1301  + dxt(i,2) * ta1(2,j,1) &
1302  + dxt(i,3) * ta1(3,j,1) &
1303  + dxt(i,4) * ta1(4,j,1) &
1304  + dxt(i,5) * ta1(5,j,1) &
1305  + dxt(i,6) * ta1(6,j,1) &
1306  + dxt(i,7) * ta1(7,j,1) &
1307  + dxt(i,8) * ta1(8,j,1) &
1308  + dxt(i,9) * ta1(9,j,1) &
1309  + dxt(i,10) * ta1(10,j,1) &
1310  + dxt(i,11) * ta1(11,j,1) &
1311  + dxt(i,12) * ta1(12,j,1) &
1312  + dxt(i,13) * ta1(13,j,1) &
1313  + dxt(i,14) * ta1(14,j,1)
1314  end do
1315  end do
1316 
1317  do i = 1, lx*lx*lx
1318  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1319  end do
1320 
1321  do k = 1, lx
1322  do j = 1, lx
1323  do i = 1, lx
1324  dtx(i,j,k) = dtx(i,j,k) &
1325  + dyt(j,1) * ta1(i,1,k) &
1326  + dyt(j,2) * ta1(i,2,k) &
1327  + dyt(j,3) * ta1(i,3,k) &
1328  + dyt(j,4) * ta1(i,4,k) &
1329  + dyt(j,5) * ta1(i,5,k) &
1330  + dyt(j,6) * ta1(i,6,k) &
1331  + dyt(j,7) * ta1(i,7,k) &
1332  + dyt(j,8) * ta1(i,8,k) &
1333  + dyt(j,9) * ta1(i,9,k) &
1334  + dyt(j,10) * ta1(i,10,k) &
1335  + dyt(j,11) * ta1(i,11,k) &
1336  + dyt(j,12) * ta1(i,12,k) &
1337  + dyt(j,13) * ta1(i,13,k) &
1338  + dyt(j,14) * ta1(i,14,k)
1339  end do
1340  end do
1341  end do
1342 
1343  do i = 1, lx*lx*lx
1344  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1345  end do
1346 
1347  do k = 1, lx
1348  do i = 1, lx*lx
1349  dtx(i,1,k) = dtx(i,1,k) &
1350  + dzt(k,1) * ta1(i,1,1) &
1351  + dzt(k,2) * ta1(i,1,2) &
1352  + dzt(k,3) * ta1(i,1,3) &
1353  + dzt(k,4) * ta1(i,1,4) &
1354  + dzt(k,5) * ta1(i,1,5) &
1355  + dzt(k,6) * ta1(i,1,6) &
1356  + dzt(k,7) * ta1(i,1,7) &
1357  + dzt(k,8) * ta1(i,1,8) &
1358  + dzt(k,9) * ta1(i,1,9) &
1359  + dzt(k,10) * ta1(i,1,10) &
1360  + dzt(k,11) * ta1(i,1,11) &
1361  + dzt(k,12) * ta1(i,1,12) &
1362  + dzt(k,13) * ta1(i,1,13) &
1363  + dzt(k,14) * ta1(i,1,14)
1364  end do
1365  end do
1366 
1367  end subroutine cpu_cdtp_lx14_single
1368 
1369  subroutine cpu_cdtp_lx13_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1370  integer, parameter :: lx = 13
1371  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1372  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1373  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1374  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1375  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1376  integer :: i, j, k
1377 
1378  do i = 1, lx*lx*lx
1379  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1380  end do
1381 
1382  do i = 1, lx*lx*lx
1383  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1384  end do
1385 
1386  do j = 1, lx * lx
1387  do i = 1, lx
1388  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1389  + dxt(i,2) * ta1(2,j,1) &
1390  + dxt(i,3) * ta1(3,j,1) &
1391  + dxt(i,4) * ta1(4,j,1) &
1392  + dxt(i,5) * ta1(5,j,1) &
1393  + dxt(i,6) * ta1(6,j,1) &
1394  + dxt(i,7) * ta1(7,j,1) &
1395  + dxt(i,8) * ta1(8,j,1) &
1396  + dxt(i,9) * ta1(9,j,1) &
1397  + dxt(i,10) * ta1(10,j,1) &
1398  + dxt(i,11) * ta1(11,j,1) &
1399  + dxt(i,12) * ta1(12,j,1) &
1400  + dxt(i,13) * ta1(13,j,1)
1401  end do
1402  end do
1403 
1404  do i = 1, lx*lx*lx
1405  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1406  end do
1407 
1408  do k = 1, lx
1409  do j = 1, lx
1410  do i = 1, lx
1411  dtx(i,j,k) = dtx(i,j,k) &
1412  + dyt(j,1) * ta1(i,1,k) &
1413  + dyt(j,2) * ta1(i,2,k) &
1414  + dyt(j,3) * ta1(i,3,k) &
1415  + dyt(j,4) * ta1(i,4,k) &
1416  + dyt(j,5) * ta1(i,5,k) &
1417  + dyt(j,6) * ta1(i,6,k) &
1418  + dyt(j,7) * ta1(i,7,k) &
1419  + dyt(j,8) * ta1(i,8,k) &
1420  + dyt(j,9) * ta1(i,9,k) &
1421  + dyt(j,10) * ta1(i,10,k) &
1422  + dyt(j,11) * ta1(i,11,k) &
1423  + dyt(j,12) * ta1(i,12,k) &
1424  + dyt(j,13) * ta1(i,13,k)
1425  end do
1426  end do
1427  end do
1428 
1429  do i = 1, lx*lx*lx
1430  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1431  end do
1432 
1433  do k = 1, lx
1434  do i = 1, lx*lx
1435  dtx(i,1,k) = dtx(i,1,k) &
1436  + dzt(k,1) * ta1(i,1,1) &
1437  + dzt(k,2) * ta1(i,1,2) &
1438  + dzt(k,3) * ta1(i,1,3) &
1439  + dzt(k,4) * ta1(i,1,4) &
1440  + dzt(k,5) * ta1(i,1,5) &
1441  + dzt(k,6) * ta1(i,1,6) &
1442  + dzt(k,7) * ta1(i,1,7) &
1443  + dzt(k,8) * ta1(i,1,8) &
1444  + dzt(k,9) * ta1(i,1,9) &
1445  + dzt(k,10) * ta1(i,1,10) &
1446  + dzt(k,11) * ta1(i,1,11) &
1447  + dzt(k,12) * ta1(i,1,12) &
1448  + dzt(k,13) * ta1(i,1,13)
1449  end do
1450  end do
1451 
1452  end subroutine cpu_cdtp_lx13_single
1453 
1454  subroutine cpu_cdtp_lx12_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1455  integer, parameter :: lx = 12
1456  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1457  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1458  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1459  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1460  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1461  integer :: i, j, k
1462 
1463  do i = 1, lx*lx*lx
1464  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1465  end do
1466 
1467  do i = 1, lx*lx*lx
1468  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1469  end do
1470 
1471  do j = 1, lx * lx
1472  do i = 1, lx
1473  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1474  + dxt(i,2) * ta1(2,j,1) &
1475  + dxt(i,3) * ta1(3,j,1) &
1476  + dxt(i,4) * ta1(4,j,1) &
1477  + dxt(i,5) * ta1(5,j,1) &
1478  + dxt(i,6) * ta1(6,j,1) &
1479  + dxt(i,7) * ta1(7,j,1) &
1480  + dxt(i,8) * ta1(8,j,1) &
1481  + dxt(i,9) * ta1(9,j,1) &
1482  + dxt(i,10) * ta1(10,j,1) &
1483  + dxt(i,11) * ta1(11,j,1) &
1484  + dxt(i,12) * ta1(12,j,1)
1485  end do
1486  end do
1487 
1488  do i = 1, lx*lx*lx
1489  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1490  end do
1491 
1492  do k = 1, lx
1493  do j = 1, lx
1494  do i = 1, lx
1495  dtx(i,j,k) = dtx(i,j,k) &
1496  + dyt(j,1) * ta1(i,1,k) &
1497  + dyt(j,2) * ta1(i,2,k) &
1498  + dyt(j,3) * ta1(i,3,k) &
1499  + dyt(j,4) * ta1(i,4,k) &
1500  + dyt(j,5) * ta1(i,5,k) &
1501  + dyt(j,6) * ta1(i,6,k) &
1502  + dyt(j,7) * ta1(i,7,k) &
1503  + dyt(j,8) * ta1(i,8,k) &
1504  + dyt(j,9) * ta1(i,9,k) &
1505  + dyt(j,10) * ta1(i,10,k) &
1506  + dyt(j,11) * ta1(i,11,k) &
1507  + dyt(j,12) * ta1(i,12,k)
1508  end do
1509  end do
1510  end do
1511 
1512  do i = 1, lx*lx*lx
1513  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1514  end do
1515 
1516  do k = 1, lx
1517  do i = 1, lx*lx
1518  dtx(i,1,k) = dtx(i,1,k) &
1519  + dzt(k,1) * ta1(i,1,1) &
1520  + dzt(k,2) * ta1(i,1,2) &
1521  + dzt(k,3) * ta1(i,1,3) &
1522  + dzt(k,4) * ta1(i,1,4) &
1523  + dzt(k,5) * ta1(i,1,5) &
1524  + dzt(k,6) * ta1(i,1,6) &
1525  + dzt(k,7) * ta1(i,1,7) &
1526  + dzt(k,8) * ta1(i,1,8) &
1527  + dzt(k,9) * ta1(i,1,9) &
1528  + dzt(k,10) * ta1(i,1,10) &
1529  + dzt(k,11) * ta1(i,1,11) &
1530  + dzt(k,12) * ta1(i,1,12)
1531  end do
1532  end do
1533 
1534  end subroutine cpu_cdtp_lx12_single
1535 
1536  subroutine cpu_cdtp_lx11_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1537  integer, parameter :: lx = 11
1538  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1539  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1540  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1541  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1542  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1543  integer :: i, j, k
1544 
1545  do i = 1, lx*lx*lx
1546  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1547  end do
1548 
1549  do i = 1, lx*lx*lx
1550  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1551  end do
1552 
1553  do j = 1, lx * lx
1554  do i = 1, lx
1555  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1556  + dxt(i,2) * ta1(2,j,1) &
1557  + dxt(i,3) * ta1(3,j,1) &
1558  + dxt(i,4) * ta1(4,j,1) &
1559  + dxt(i,5) * ta1(5,j,1) &
1560  + dxt(i,6) * ta1(6,j,1) &
1561  + dxt(i,7) * ta1(7,j,1) &
1562  + dxt(i,8) * ta1(8,j,1) &
1563  + dxt(i,9) * ta1(9,j,1) &
1564  + dxt(i,10) * ta1(10,j,1) &
1565  + dxt(i,11) * ta1(11,j,1)
1566  end do
1567  end do
1568 
1569  do i = 1, lx*lx*lx
1570  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1571  end do
1572 
1573  do k = 1, lx
1574  do j = 1, lx
1575  do i = 1, lx
1576  dtx(i,j,k) = dtx(i,j,k) &
1577  + dyt(j,1) * ta1(i,1,k) &
1578  + dyt(j,2) * ta1(i,2,k) &
1579  + dyt(j,3) * ta1(i,3,k) &
1580  + dyt(j,4) * ta1(i,4,k) &
1581  + dyt(j,5) * ta1(i,5,k) &
1582  + dyt(j,6) * ta1(i,6,k) &
1583  + dyt(j,7) * ta1(i,7,k) &
1584  + dyt(j,8) * ta1(i,8,k) &
1585  + dyt(j,9) * ta1(i,9,k) &
1586  + dyt(j,10) * ta1(i,10,k) &
1587  + dyt(j,11) * ta1(i,11,k)
1588  end do
1589  end do
1590  end do
1591 
1592  do i = 1, lx*lx*lx
1593  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1594  end do
1595 
1596  do k = 1, lx
1597  do i = 1, lx*lx
1598  dtx(i,1,k) = dtx(i,1,k) &
1599  + dzt(k,1) * ta1(i,1,1) &
1600  + dzt(k,2) * ta1(i,1,2) &
1601  + dzt(k,3) * ta1(i,1,3) &
1602  + dzt(k,4) * ta1(i,1,4) &
1603  + dzt(k,5) * ta1(i,1,5) &
1604  + dzt(k,6) * ta1(i,1,6) &
1605  + dzt(k,7) * ta1(i,1,7) &
1606  + dzt(k,8) * ta1(i,1,8) &
1607  + dzt(k,9) * ta1(i,1,9) &
1608  + dzt(k,10) * ta1(i,1,10) &
1609  + dzt(k,11) * ta1(i,1,11)
1610  end do
1611  end do
1612 
1613  end subroutine cpu_cdtp_lx11_single
1614 
1615  subroutine cpu_cdtp_lx10_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1616  integer, parameter :: lx = 10
1617  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1618  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1619  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1620  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1621  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1622  integer :: i, j, k
1623 
1624  do i = 1, lx*lx*lx
1625  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1626  end do
1627 
1628  do i = 1, lx*lx*lx
1629  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1630  end do
1631 
1632  do j = 1, lx * lx
1633  do i = 1, lx
1634  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1635  + dxt(i,2) * ta1(2,j,1) &
1636  + dxt(i,3) * ta1(3,j,1) &
1637  + dxt(i,4) * ta1(4,j,1) &
1638  + dxt(i,5) * ta1(5,j,1) &
1639  + dxt(i,6) * ta1(6,j,1) &
1640  + dxt(i,7) * ta1(7,j,1) &
1641  + dxt(i,8) * ta1(8,j,1) &
1642  + dxt(i,9) * ta1(9,j,1) &
1643  + dxt(i,10) * ta1(10,j,1)
1644  end do
1645  end do
1646 
1647  do i = 1, lx*lx*lx
1648  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1649  end do
1650 
1651  do k = 1, lx
1652  do j = 1, lx
1653  do i = 1, lx
1654  dtx(i,j,k) = dtx(i,j,k) &
1655  + dyt(j,1) * ta1(i,1,k) &
1656  + dyt(j,2) * ta1(i,2,k) &
1657  + dyt(j,3) * ta1(i,3,k) &
1658  + dyt(j,4) * ta1(i,4,k) &
1659  + dyt(j,5) * ta1(i,5,k) &
1660  + dyt(j,6) * ta1(i,6,k) &
1661  + dyt(j,7) * ta1(i,7,k) &
1662  + dyt(j,8) * ta1(i,8,k) &
1663  + dyt(j,9) * ta1(i,9,k) &
1664  + dyt(j,10) * ta1(i,10,k)
1665 
1666  end do
1667  end do
1668  end do
1669 
1670  do i = 1, lx*lx*lx
1671  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1672  end do
1673 
1674  do k = 1, lx
1675  do i = 1, lx*lx
1676  dtx(i,1,k) = dtx(i,1,k) &
1677  + dzt(k,1) * ta1(i,1,1) &
1678  + dzt(k,2) * ta1(i,1,2) &
1679  + dzt(k,3) * ta1(i,1,3) &
1680  + dzt(k,4) * ta1(i,1,4) &
1681  + dzt(k,5) * ta1(i,1,5) &
1682  + dzt(k,6) * ta1(i,1,6) &
1683  + dzt(k,7) * ta1(i,1,7) &
1684  + dzt(k,8) * ta1(i,1,8) &
1685  + dzt(k,9) * ta1(i,1,9) &
1686  + dzt(k,10) * ta1(i,1,10)
1687  end do
1688  end do
1689 
1690  end subroutine cpu_cdtp_lx10_single
1691 
1692  subroutine cpu_cdtp_lx9_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1693  integer, parameter :: lx = 9
1694  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1695  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1696  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1697  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1698  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1699  integer :: i, j, k
1700 
1701  do i = 1, lx*lx*lx
1702  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1703  end do
1704 
1705  do i = 1, lx*lx*lx
1706  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1707  end do
1708 
1709  do j = 1, lx * lx
1710  do i = 1, lx
1711  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1712  + dxt(i,2) * ta1(2,j,1) &
1713  + dxt(i,3) * ta1(3,j,1) &
1714  + dxt(i,4) * ta1(4,j,1) &
1715  + dxt(i,5) * ta1(5,j,1) &
1716  + dxt(i,6) * ta1(6,j,1) &
1717  + dxt(i,7) * ta1(7,j,1) &
1718  + dxt(i,8) * ta1(8,j,1) &
1719  + dxt(i,9) * ta1(9,j,1)
1720  end do
1721  end do
1722 
1723  do i = 1, lx*lx*lx
1724  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1725  end do
1726 
1727  do k = 1, lx
1728  do j = 1, lx
1729  do i = 1, lx
1730  dtx(i,j,k) = dtx(i,j,k) &
1731  + dyt(j,1) * ta1(i,1,k) &
1732  + dyt(j,2) * ta1(i,2,k) &
1733  + dyt(j,3) * ta1(i,3,k) &
1734  + dyt(j,4) * ta1(i,4,k) &
1735  + dyt(j,5) * ta1(i,5,k) &
1736  + dyt(j,6) * ta1(i,6,k) &
1737  + dyt(j,7) * ta1(i,7,k) &
1738  + dyt(j,8) * ta1(i,8,k) &
1739  + dyt(j,9) * ta1(i,9,k)
1740  end do
1741  end do
1742  end do
1743 
1744  do i = 1, lx*lx*lx
1745  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1746  end do
1747 
1748  do k = 1, lx
1749  do i = 1, lx*lx
1750  dtx(i,1,k) = dtx(i,1,k) &
1751  + dzt(k,1) * ta1(i,1,1) &
1752  + dzt(k,2) * ta1(i,1,2) &
1753  + dzt(k,3) * ta1(i,1,3) &
1754  + dzt(k,4) * ta1(i,1,4) &
1755  + dzt(k,5) * ta1(i,1,5) &
1756  + dzt(k,6) * ta1(i,1,6) &
1757  + dzt(k,7) * ta1(i,1,7) &
1758  + dzt(k,8) * ta1(i,1,8) &
1759  + dzt(k,9) * ta1(i,1,9)
1760  end do
1761  end do
1762 
1763  end subroutine cpu_cdtp_lx9_single
1764 
1765  subroutine cpu_cdtp_lx8_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1766  integer, parameter :: lx = 8
1767  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1768  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1769  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1770  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1771  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1772  integer :: i, j, k
1773 
1774  do i = 1, lx*lx*lx
1775  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1776  end do
1777 
1778  do i = 1, lx*lx*lx
1779  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1780  end do
1781 
1782  do j = 1, lx * lx
1783  do i = 1, lx
1784  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1785  + dxt(i,2) * ta1(2,j,1) &
1786  + dxt(i,3) * ta1(3,j,1) &
1787  + dxt(i,4) * ta1(4,j,1) &
1788  + dxt(i,5) * ta1(5,j,1) &
1789  + dxt(i,6) * ta1(6,j,1) &
1790  + dxt(i,7) * ta1(7,j,1) &
1791  + dxt(i,8) * ta1(8,j,1)
1792  end do
1793  end do
1794 
1795  do i = 1, lx*lx*lx
1796  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1797  end do
1798 
1799  do k = 1, lx
1800  do j = 1, lx
1801  do i = 1, lx
1802  dtx(i,j,k) = dtx(i,j,k) &
1803  + dyt(j,1) * ta1(i,1,k) &
1804  + dyt(j,2) * ta1(i,2,k) &
1805  + dyt(j,3) * ta1(i,3,k) &
1806  + dyt(j,4) * ta1(i,4,k) &
1807  + dyt(j,5) * ta1(i,5,k) &
1808  + dyt(j,6) * ta1(i,6,k) &
1809  + dyt(j,7) * ta1(i,7,k) &
1810  + dyt(j,8) * ta1(i,8,k)
1811  end do
1812  end do
1813  end do
1814 
1815  do i = 1, lx*lx*lx
1816  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1817  end do
1818 
1819  do k = 1, lx
1820  do i = 1, lx*lx
1821  dtx(i,1,k) = dtx(i,1,k) &
1822  + dzt(k,1) * ta1(i,1,1) &
1823  + dzt(k,2) * ta1(i,1,2) &
1824  + dzt(k,3) * ta1(i,1,3) &
1825  + dzt(k,4) * ta1(i,1,4) &
1826  + dzt(k,5) * ta1(i,1,5) &
1827  + dzt(k,6) * ta1(i,1,6) &
1828  + dzt(k,7) * ta1(i,1,7) &
1829  + dzt(k,8) * ta1(i,1,8)
1830  end do
1831  end do
1832 
1833  end subroutine cpu_cdtp_lx8_single
1834 
1835  subroutine cpu_cdtp_lx7_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1836  integer, parameter :: lx = 7
1837  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1838  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1839  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1840  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1841  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1842  integer :: i, j, k
1843 
1844  do i = 1, lx*lx*lx
1845  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1846  end do
1847 
1848  do i = 1, lx*lx*lx
1849  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1850  end do
1851 
1852  do j = 1, lx * lx
1853  do i = 1, lx
1854  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1855  + dxt(i,2) * ta1(2,j,1) &
1856  + dxt(i,3) * ta1(3,j,1) &
1857  + dxt(i,4) * ta1(4,j,1) &
1858  + dxt(i,5) * ta1(5,j,1) &
1859  + dxt(i,6) * ta1(6,j,1) &
1860  + dxt(i,7) * ta1(7,j,1)
1861  end do
1862  end do
1863 
1864  do i = 1, lx*lx*lx
1865  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1866  end do
1867 
1868  do k = 1, lx
1869  do j = 1, lx
1870  do i = 1, lx
1871  dtx(i,j,k) = dtx(i,j,k) &
1872  + dyt(j,1) * ta1(i,1,k) &
1873  + dyt(j,2) * ta1(i,2,k) &
1874  + dyt(j,3) * ta1(i,3,k) &
1875  + dyt(j,4) * ta1(i,4,k) &
1876  + dyt(j,5) * ta1(i,5,k) &
1877  + dyt(j,6) * ta1(i,6,k) &
1878  + dyt(j,7) * ta1(i,7,k)
1879  end do
1880  end do
1881  end do
1882 
1883  do i = 1, lx*lx*lx
1884  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1885  end do
1886 
1887  do k = 1, lx
1888  do i = 1, lx*lx
1889  dtx(i,1,k) = dtx(i,1,k) &
1890  + dzt(k,1) * ta1(i,1,1) &
1891  + dzt(k,2) * ta1(i,1,2) &
1892  + dzt(k,3) * ta1(i,1,3) &
1893  + dzt(k,4) * ta1(i,1,4) &
1894  + dzt(k,5) * ta1(i,1,5) &
1895  + dzt(k,6) * ta1(i,1,6) &
1896  + dzt(k,7) * ta1(i,1,7)
1897  end do
1898  end do
1899 
1900  end subroutine cpu_cdtp_lx7_single
1901 
1902  subroutine cpu_cdtp_lx6_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1903  integer, parameter :: lx = 6
1904  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1905  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1906  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1907  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1908  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1909  integer :: i, j, k
1910 
1911  do i = 1, lx*lx*lx
1912  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1913  end do
1914 
1915  do i = 1, lx*lx*lx
1916  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1917  end do
1918 
1919  do j = 1, lx * lx
1920  do i = 1, lx
1921  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1922  + dxt(i,2) * ta1(2,j,1) &
1923  + dxt(i,3) * ta1(3,j,1) &
1924  + dxt(i,4) * ta1(4,j,1) &
1925  + dxt(i,5) * ta1(5,j,1) &
1926  + dxt(i,6) * ta1(6,j,1)
1927  end do
1928  end do
1929 
1930  do i = 1, lx*lx*lx
1931  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1932  end do
1933 
1934  do k = 1, lx
1935  do j = 1, lx
1936  do i = 1, lx
1937  dtx(i,j,k) = dtx(i,j,k) &
1938  + dyt(j,1) * ta1(i,1,k) &
1939  + dyt(j,2) * ta1(i,2,k) &
1940  + dyt(j,3) * ta1(i,3,k) &
1941  + dyt(j,4) * ta1(i,4,k) &
1942  + dyt(j,5) * ta1(i,5,k) &
1943  + dyt(j,6) * ta1(i,6,k)
1944  end do
1945  end do
1946  end do
1947 
1948  do i = 1, lx*lx*lx
1949  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
1950  end do
1951 
1952  do k = 1, lx
1953  do i = 1, lx*lx
1954  dtx(i,1,k) = dtx(i,1,k) &
1955  + dzt(k,1) * ta1(i,1,1) &
1956  + dzt(k,2) * ta1(i,1,2) &
1957  + dzt(k,3) * ta1(i,1,3) &
1958  + dzt(k,4) * ta1(i,1,4) &
1959  + dzt(k,5) * ta1(i,1,5) &
1960  + dzt(k,6) * ta1(i,1,6)
1961  end do
1962  end do
1963 
1964  end subroutine cpu_cdtp_lx6_single
1965 
1966  subroutine cpu_cdtp_lx5_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
1967  integer, parameter :: lx = 5
1968  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
1969  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
1970  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
1971  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
1972  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
1973  integer :: i, j, k
1974 
1975  do i = 1, lx*lx*lx
1976  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
1977  end do
1978 
1979  do i = 1, lx*lx*lx
1980  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
1981  end do
1982 
1983  do j = 1, lx * lx
1984  do i = 1, lx
1985  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
1986  + dxt(i,2) * ta1(2,j,1) &
1987  + dxt(i,3) * ta1(3,j,1) &
1988  + dxt(i,4) * ta1(4,j,1) &
1989  + dxt(i,5) * ta1(5,j,1)
1990  end do
1991  end do
1992 
1993  do i = 1, lx*lx*lx
1994  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
1995  end do
1996 
1997  do k = 1, lx
1998  do j = 1, lx
1999  do i = 1, lx
2000  dtx(i,j,k) = dtx(i,j,k) &
2001  + dyt(j,1) * ta1(i,1,k) &
2002  + dyt(j,2) * ta1(i,2,k) &
2003  + dyt(j,3) * ta1(i,3,k) &
2004  + dyt(j,4) * ta1(i,4,k) &
2005  + dyt(j,5) * ta1(i,5,k)
2006  end do
2007  end do
2008  end do
2009 
2010  do i = 1, lx*lx*lx
2011  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
2012  end do
2013 
2014  do k = 1, lx
2015  do i = 1, lx*lx
2016  dtx(i,1,k) = dtx(i,1,k) &
2017  + dzt(k,1) * ta1(i,1,1) &
2018  + dzt(k,2) * ta1(i,1,2) &
2019  + dzt(k,3) * ta1(i,1,3) &
2020  + dzt(k,4) * ta1(i,1,4) &
2021  + dzt(k,5) * ta1(i,1,5)
2022  end do
2023  end do
2024 
2025  end subroutine cpu_cdtp_lx5_single
2026 
2027  subroutine cpu_cdtp_lx4_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
2028  integer, parameter :: lx = 4
2029  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
2030  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
2031  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
2032  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
2033  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
2034  integer :: i, j, k
2035 
2036  do i = 1, lx*lx*lx
2037  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
2038  end do
2039 
2040  do i = 1, lx*lx*lx
2041  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
2042  end do
2043 
2044  do j = 1, lx * lx
2045  do i = 1, lx
2046  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
2047  + dxt(i,2) * ta1(2,j,1) &
2048  + dxt(i,3) * ta1(3,j,1) &
2049  + dxt(i,4) * ta1(4,j,1)
2050  end do
2051  end do
2052 
2053  do i = 1, lx*lx*lx
2054  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
2055  end do
2056 
2057  do k = 1, lx
2058  do j = 1, lx
2059  do i = 1, lx
2060  dtx(i,j,k) = dtx(i,j,k) &
2061  + dyt(j,1) * ta1(i,1,k) &
2062  + dyt(j,2) * ta1(i,2,k) &
2063  + dyt(j,3) * ta1(i,3,k) &
2064  + dyt(j,4) * ta1(i,4,k)
2065  end do
2066  end do
2067  end do
2068 
2069  do i = 1, lx*lx*lx
2070  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
2071  end do
2072 
2073  do k = 1, lx
2074  do i = 1, lx*lx
2075  dtx(i,1,k) = dtx(i,1,k) &
2076  + dzt(k,1) * ta1(i,1,1) &
2077  + dzt(k,2) * ta1(i,1,2) &
2078  + dzt(k,3) * ta1(i,1,3) &
2079  + dzt(k,4) * ta1(i,1,4)
2080  end do
2081  end do
2082 
2083  end subroutine cpu_cdtp_lx4_single
2084 
2085  subroutine cpu_cdtp_lx3_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
2086  integer, parameter :: lx = 3
2087  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
2088  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
2089  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
2090  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
2091  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
2092  integer :: i, j, k
2093 
2094  do i = 1, lx*lx*lx
2095  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
2096  end do
2097 
2098  do i = 1, lx*lx*lx
2099  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
2100  end do
2101 
2102  do j = 1, lx * lx
2103  do i = 1, lx
2104  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
2105  + dxt(i,2) * ta1(2,j,1) &
2106  + dxt(i,3) * ta1(3,j,1)
2107  end do
2108  end do
2109 
2110  do i = 1, lx*lx*lx
2111  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
2112  end do
2113 
2114  do k = 1, lx
2115  do j = 1, lx
2116  do i = 1, lx
2117  dtx(i,j,k) = dtx(i,j,k) &
2118  + dyt(j,1) * ta1(i,1,k) &
2119  + dyt(j,2) * ta1(i,2,k) &
2120  + dyt(j,3) * ta1(i,3,k)
2121  end do
2122  end do
2123  end do
2124 
2125  do i = 1, lx*lx*lx
2126  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
2127  end do
2128 
2129  do k = 1, lx
2130  do i = 1, lx*lx
2131  dtx(i,1,k) = dtx(i,1,k) &
2132  + dzt(k,1) * ta1(i,1,1) &
2133  + dzt(k,2) * ta1(i,1,2) &
2134  + dzt(k,3) * ta1(i,1,3)
2135  end do
2136  end do
2137 
2138  end subroutine cpu_cdtp_lx3_single
2139 
2140  subroutine cpu_cdtp_lx2_single(dtx, x, dr, ds, dt, dxt, dyt, dzt, w3)
2141  integer, parameter :: lx = 2
2142  real(kind=rp), dimension(lx, lx, lx), intent(inout) :: dtx
2143  real(kind=rp), dimension(lx, lx, lx), intent(in) :: x, dr, ds, dt
2144  real(kind=rp), dimension(lx, lx, lx), intent(in) :: w3
2145  real(kind=rp), intent(in), dimension(lx, lx) :: dxt, dyt, dzt
2146  real(kind=rp), dimension(lx, lx, lx) :: wx, ta1
2147  integer :: i, j, k
2148 
2149  do i = 1, lx*lx*lx
2150  wx(i,1,1) = x(i,1,1) * w3(i,1,1)
2151  end do
2152 
2153  do i = 1, lx*lx*lx
2154  ta1(i,1,1) = wx(i,1,1) * dr(i,1,1)
2155  end do
2156 
2157  do j = 1, lx * lx
2158  do i = 1, lx
2159  dtx(i,j,1) = dxt(i,1) * ta1(1,j,1) &
2160  + dxt(i,2) * ta1(2,j,1)
2161  end do
2162  end do
2163 
2164  do i = 1, lx*lx*lx
2165  ta1(i,1,1) = wx(i,1,1) * ds(i,1,1)
2166  end do
2167 
2168  do k = 1, lx
2169  do j = 1, lx
2170  do i = 1, lx
2171  dtx(i,j,k) = dtx(i,j,k) &
2172  + dyt(j,1) * ta1(i,1,k) &
2173  + dyt(j,2) * ta1(i,2,k)
2174  end do
2175  end do
2176  end do
2177 
2178  do i = 1, lx*lx*lx
2179  ta1(i,1,1) = wx(i,1,1) * dt(i,1,1)
2180  end do
2181 
2182  do k = 1, lx
2183  do i = 1, lx*lx
2184  dtx(i,1,k) = dtx(i,1,k) &
2185  + dzt(k,1) * ta1(i,1,1) &
2186  + dzt(k,2) * ta1(i,1,2)
2187  end do
2188  end do
2189 
2190  end subroutine cpu_cdtp_lx2_single
2191 end submodule cpu_cdtp
opr_cpu
Operators CPU backend.
Definition: opr_cpu.f90:34