Neko  0.9.99
A portable framework for high-order spectral element flow simulations
coriolis_source_term_cpu.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 !
35 
37  use num_types, only : rp
38  use field_list, only : field_list_t
39  use math, only : vcross
40  use field, only : field_t
41  implicit none
42  private
43 
45 
46 contains
47 
52  subroutine coriolis_source_term_compute_cpu(fields, omega, u_geo)
53  type(field_list_t), intent(inout) :: fields
54  real(kind=rp), intent(in) :: omega(3)
55  real(kind=rp), intent(in) :: u_geo(3)
56  integer :: i, n
57  type(field_t), pointer :: u, v, w
58  real(kind=rp) :: ui, vi, wi
59 
60 
61  n = fields%item_size(1)
62 
63  u => fields%get_by_index(1)
64  v => fields%get_by_index(2)
65  w => fields%get_by_index(3)
66 
67  do concurrent(i = 1:n)
68  ui = u%x(i,1,1,1) - u_geo(1)
69  vi = v%x(i,1,1,1) - u_geo(2)
70  wi = w%x(i,1,1,1) - u_geo(3)
71 
72  u%x(i,1,1,1) = u%x(i,1,1,1) - 2.0_rp * (omega(2) * wi - omega(3) * vi)
73  v%x(i,1,1,1) = v%x(i,1,1,1) - 2.0_rp * (omega(3) * ui - omega(1) * wi)
74  w%x(i,1,1,1) = w%x(i,1,1,1) - 2.0_rp * (omega(1) * vi - omega(2) * ui)
75  end do
76 
78 
79 end module coriolis_source_term_cpu
Implements the cpu kernel for the coriolis_source_term_t type. Maintainer: Timofey Mukha.
subroutine, public coriolis_source_term_compute_cpu(fields, omega, u_geo)
Computes the generic Coriolis source term on the cpu.
Defines a field.
Definition: field.f90:34
Definition: math.f90:60
subroutine, public vcross(u1, u2, u3, v1, v2, v3, w1, w2, w3, n)
Compute a cross product assuming vector components etc.
Definition: math.f90:514
integer, parameter, public rp
Global precision used in computations.
Definition: num_types.f90:12
field_list_t, To be able to group fields together
Definition: field_list.f90:13