MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_body_forces.fpp.f90
Go to the documentation of this file.
1# 1 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
2!>
3!! @file
4!! @brief Contains module m_body_forces
5
6# 1 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 1
7# 1 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 1
8# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
9# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
10# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
11# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
12# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
13# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
14
15# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
16# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
17# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
18
19# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
20
21# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
22
23# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
24
25# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
26
27# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
28
29# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
30
31# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
32! New line at end of file is required for FYPP
33# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
34# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
35# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
36# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
37# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
38# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
39# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
40# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
41
42# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
43# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
44# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
45
46# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
47
48# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
49
50# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
51
52# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
53
54# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
55
56# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
57
58# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
59! New line at end of file is required for FYPP
60# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
61
62# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
63# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
64# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
65# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
66# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
67
68# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
69
70# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
71
72# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
73
74# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
75
76# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
77
78# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
79
80# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
81
82# 76 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
83
84# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
85
86# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
87
88# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
89
90# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
91
92# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
93
94# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
95
96# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
97
98# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
99
100# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
101
102# 151 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
103
104# 192 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
105
106# 207 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
107
108# 232 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
109
110# 243 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
111
112# 245 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
113# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
114
115# 283 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
116
117# 293 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
118
119# 303 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
120
121# 312 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
122
123# 329 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
124
125# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
126
127# 346 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
128
129# 352 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
130
131# 358 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
132
133# 364 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
134
135# 370 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
136
137# 376 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
138! New line at end of file is required for FYPP
139# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
140# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
141# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
142# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
143# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
144# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
145# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
146# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
147
148# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
149# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
150# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
151
152# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
153
154# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
155
156# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
157
158# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
159
160# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
161
162# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
163
164# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
165! New line at end of file is required for FYPP
166# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
167
168# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
169
170# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
171
172# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
173
174# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
175
176# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
177
178# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
179
180# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
181
182# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
183
184# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
185
186# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
187
188# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
189
190# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
191
192# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
193
194# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
195
196# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
197
198# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
199
200# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
201
202# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
203
204# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
205
206# 192 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
207
208# 213 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
209
210# 241 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
211
212# 256 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
213
214# 266 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
215
216# 275 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
217
218# 291 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
219
220# 301 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
221
222# 308 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
223! New line at end of file is required for FYPP
224# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
225
226# 21 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
227
228# 37 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
229
230# 50 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
231
232# 76 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
233
234# 91 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
235
236# 102 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
237
238# 115 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
239
240# 143 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
241
242# 154 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
243
244# 165 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
245
246# 176 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
247
248# 187 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
249
250# 198 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
251
252# 208 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
253
254# 214 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
255
256# 220 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
257
258# 226 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
259
260# 232 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
261
262# 234 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
263# 235 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
264! New line at end of file is required for FYPP
265# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
266
267# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
268
269! Caution:
270! This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI rank.
271! That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0.
272! For an example see misc/nvidia_uvm/bind.sh.
273# 63 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
274
275# 81 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
276
277# 88 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
278
279# 111 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
280
281# 127 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
282
283# 153 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
284
285# 159 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
286
287# 167 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
288! New line at end of file is required for FYPP
289# 6 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp" 2
290
291!> @brief Computes gravitational and user-defined body force source terms for the momentum equations
292module m_body_forces
293
294 use m_derived_types !< definitions of the derived types
295
296 use m_global_parameters !< definitions of the global parameters
297
298 use m_variables_conversion
299
300 use m_nvtx
301
302! $:USE_GPU_MODULE()
303
304 implicit none
305
306 private;
307 public :: s_compute_body_forces_rhs, &
308 s_initialize_body_forces_module, &
309 s_finalize_body_forces_module
310
311 real(wp), allocatable, dimension(:, :, :) :: rhom
312
313# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
314#if defined(MFC_OpenACC)
315# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
316!$acc declare create(rhoM)
317# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
318#elif defined(MFC_OpenMP)
319# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
320!$omp declare target (rhoM)
321# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
322#endif
323
324contains
325
326 !> This subroutine initializes the module global array of mixture
327 !! densities in each grid cell
328 impure subroutine s_initialize_body_forces_module
329
330 ! Simulation is at least 2D
331 if (n > 0) then
332 ! Simulation is 3D
333 if (p > 0) then
334#ifdef MFC_DEBUG
335# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
336 block
337# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
338 use iso_fortran_env, only: output_unit
339# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
340
341# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
342 print *, 'm_body_forces.fpp:40: ', '@:ALLOCATE(rhoM(-buff_size:buff_size + m, -buff_size:buff_size + n, -buff_size:buff_size + p))'
343# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
344
345# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
346 call flush (output_unit)
347# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
348 end block
349# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
350#endif
351# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
352 allocate (rhom(-buff_size:buff_size + m, -buff_size:buff_size + n, -buff_size:buff_size + p))
353# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
354
355# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
356
357# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
358#if defined(MFC_OpenACC)
359# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
360!$acc enter data create(rhoM)
361# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
362#elif defined(MFC_OpenMP)
363# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
364!$omp target enter data map(always,alloc:rhoM)
365# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
366#endif
367# 43 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
368 ! Simulation is 2D
369 else
370#ifdef MFC_DEBUG
371# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
372 block
373# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
374 use iso_fortran_env, only: output_unit
375# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
376
377# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
378 print *, 'm_body_forces.fpp:45: ', '@:ALLOCATE(rhoM(-buff_size:buff_size + m, -buff_size:buff_size + n, 0:0))'
379# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
380
381# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
382 call flush (output_unit)
383# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
384 end block
385# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
386#endif
387# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
388 allocate (rhom(-buff_size:buff_size + m, -buff_size:buff_size + n, 0:0))
389# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
390
391# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
392
393# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
394#if defined(MFC_OpenACC)
395# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
396!$acc enter data create(rhoM)
397# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
398#elif defined(MFC_OpenMP)
399# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
400!$omp target enter data map(always,alloc:rhoM)
401# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
402#endif
403# 48 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
404 end if
405 ! Simulation is 1D
406 else
407#ifdef MFC_DEBUG
408# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
409 block
410# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
411 use iso_fortran_env, only: output_unit
412# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
413
414# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
415 print *, 'm_body_forces.fpp:51: ', '@:ALLOCATE(rhoM(-buff_size:buff_size + m, 0:0, 0:0))'
416# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
417
418# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
419 call flush (output_unit)
420# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
421 end block
422# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
423#endif
424# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
425 allocate (rhom(-buff_size:buff_size + m, 0:0, 0:0))
426# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
427
428# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
429
430# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
431#if defined(MFC_OpenACC)
432# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
433!$acc enter data create(rhoM)
434# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
435#elif defined(MFC_OpenMP)
436# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
437!$omp target enter data map(always,alloc:rhoM)
438# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
439#endif
440# 54 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
441 end if
442
443 end subroutine s_initialize_body_forces_module
444
445 !> This subroutine computes the acceleration at time t
446 subroutine s_compute_acceleration(t)
447
448 real(wp), intent(in) :: t
449
450# 64 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
451 if (bf_x) then
452 accel_bf(1) = g_x+k_x*sin(w_x*t - p_x)
453 end if
454# 64 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
455 if (bf_y) then
456 accel_bf(2) = g_y+k_y*sin(w_y*t - p_y)
457 end if
458# 64 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
459 if (bf_z) then
460 accel_bf(3) = g_z+k_z*sin(w_z*t - p_z)
461 end if
462# 68 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
463
464
465# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
466#if defined(MFC_OpenACC)
467# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
468!$acc update device(accel_bf)
469# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
470#elif defined(MFC_OpenMP)
471# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
472!$omp target update to(accel_bf)
473# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
474#endif
475
476 end subroutine s_compute_acceleration
477
478 !> This subroutine calculates the mixture density at each cell
479 !! center
480 !! param q_cons_vf Conservative variable
481 subroutine s_compute_mixture_density(q_cons_vf)
482
483 type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
484 integer :: i, j, k, l !< standard iterators
485
486
487# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
488
489# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
490#if defined(MFC_OpenACC)
491# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
492!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l)
493# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
494#elif defined(MFC_OpenMP)
495# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
496
497# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
498
499# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
500
501# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
502!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l)
503# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
504#endif
505# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
506
507 do l = 0, p
508 do k = 0, n
509 do j = 0, m
510 rhom(j, k, l) = 0._wp
511 do i = 1, num_fluids
512 rhom(j, k, l) = rhom(j, k, l) + &
513 q_cons_vf(contxb + i - 1)%sf(j, k, l)
514 end do
515 end do
516 end do
517 end do
518
519# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
520
521# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
522#if defined(MFC_OpenACC)
523# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
524!$acc end parallel loop
525# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
526#elif defined(MFC_OpenMP)
527# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
528
529# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
530
531# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
532!$omp end target teams loop
533# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
534#endif
535# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
536
537
538 end subroutine s_compute_mixture_density
539
540 !> This subroutine calculates the source term due to body forces
541 !! so the system can be advanced in time
542 !! @param q_cons_vf Conservative variables
543 !! @param q_prim_vf Primitive variables
544 !! @param rhs_vf Right-hand side accumulator
545 subroutine s_compute_body_forces_rhs(q_prim_vf, q_cons_vf, rhs_vf)
546
547 type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
548 type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
549 type(scalar_field), dimension(sys_size), intent(inout) :: rhs_vf
550
551 integer :: i, j, k, l !< Loop variables
552
553 call s_compute_acceleration(mytime)
554 call s_compute_mixture_density(q_cons_vf)
555
556
557# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
558
559# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
560#if defined(MFC_OpenACC)
561# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
562!$acc parallel loop collapse(4) gang vector default(present) private(i, j, k, l)
563# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
564#elif defined(MFC_OpenMP)
565# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
566
567# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
568
569# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
570
571# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
572!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(4) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k, l)
573# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
574#endif
575# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
576
577 do i = momxb, e_idx
578 do l = 0, p
579 do k = 0, n
580 do j = 0, m
581 rhs_vf(i)%sf(j, k, l) = 0._wp
582 end do
583 end do
584 end do
585 end do
586
587# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
588
589# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
590#if defined(MFC_OpenACC)
591# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
592!$acc end parallel loop
593# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
594#elif defined(MFC_OpenMP)
595# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
596
597# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
598
599# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
600!$omp end target teams loop
601# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
602#endif
603# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
604
605
606 if (bf_x) then ! x-direction body forces
607
608
609# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
610
611# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
612#if defined(MFC_OpenACC)
613# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
614!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l)
615# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
616#elif defined(MFC_OpenMP)
617# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
618
619# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
620
621# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
622
623# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
624!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l)
625# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
626#endif
627# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
628
629 do l = 0, p
630 do k = 0, n
631 do j = 0, m
632 rhs_vf(momxb)%sf(j, k, l) = rhs_vf(momxb)%sf(j, k, l) + &
633 rhom(j, k, l)*accel_bf(1)
634 rhs_vf(e_idx)%sf(j, k, l) = rhs_vf(e_idx)%sf(j, k, l) + &
635 q_cons_vf(momxb)%sf(j, k, l)*accel_bf(1)
636 end do
637 end do
638 end do
639
640# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
641
642# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
643#if defined(MFC_OpenACC)
644# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
645!$acc end parallel loop
646# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
647#elif defined(MFC_OpenMP)
648# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
649
650# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
651
652# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
653!$omp end target teams loop
654# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
655#endif
656# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
657
658 end if
659
660 if (bf_y) then ! y-direction body forces
661
662
663# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
664
665# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
666#if defined(MFC_OpenACC)
667# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
668!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l)
669# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
670#elif defined(MFC_OpenMP)
671# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
672
673# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
674
675# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
676
677# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
678!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l)
679# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
680#endif
681# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
682
683 do l = 0, p
684 do k = 0, n
685 do j = 0, m
686 rhs_vf(momxb + 1)%sf(j, k, l) = rhs_vf(momxb + 1)%sf(j, k, l) + &
687 rhom(j, k, l)*accel_bf(2)
688 rhs_vf(e_idx)%sf(j, k, l) = rhs_vf(e_idx)%sf(j, k, l) + &
689 q_cons_vf(momxb + 1)%sf(j, k, l)*accel_bf(2)
690 end do
691 end do
692 end do
693
694# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
695
696# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
697#if defined(MFC_OpenACC)
698# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
699!$acc end parallel loop
700# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
701#elif defined(MFC_OpenMP)
702# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
703
704# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
705
706# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
707!$omp end target teams loop
708# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
709#endif
710# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
711
712 end if
713
714 if (bf_z) then ! z-direction body forces
715
716
717# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
718
719# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
720#if defined(MFC_OpenACC)
721# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
722!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l)
723# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
724#elif defined(MFC_OpenMP)
725# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
726
727# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
728
729# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
730
731# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
732!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l)
733# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
734#endif
735# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
736
737 do l = 0, p
738 do k = 0, n
739 do j = 0, m
740 rhs_vf(momxe)%sf(j, k, l) = rhs_vf(momxe)%sf(j, k, l) + &
741 rhom(j, k, l)*accel_bf(3)
742 rhs_vf(e_idx)%sf(j, k, l) = rhs_vf(e_idx)%sf(j, k, l) + &
743 q_cons_vf(momxe)%sf(j, k, l)*accel_bf(3)
744 end do
745 end do
746 end do
747
748# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
749
750# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
751#if defined(MFC_OpenACC)
752# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
753!$acc end parallel loop
754# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
755#elif defined(MFC_OpenMP)
756# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
757
758# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
759
760# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
761!$omp end target teams loop
762# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
763#endif
764# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
765
766
767 end if
768
769 end subroutine s_compute_body_forces_rhs
770
771 !> @brief Deallocates module variables used for body force computations.
772 impure subroutine s_finalize_body_forces_module
773
774#ifdef MFC_DEBUG
775# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
776 block
777# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
778 use iso_fortran_env, only: output_unit
779# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
780
781# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
782 print *, 'm_body_forces.fpp:179: ', '@:DEALLOCATE(rhoM)'
783# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
784
785# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
786 call flush (output_unit)
787# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
788 end block
789# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
790#endif
791# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
792
793# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
794#if defined(MFC_OpenACC)
795# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
796!$acc exit data delete(rhoM)
797# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
798#elif defined(MFC_OpenMP)
799# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
800!$omp target exit data map(release:rhoM)
801# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
802#endif
803# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
804 deallocate (rhom)
805
806 end subroutine s_finalize_body_forces_module
807
808end module m_body_forces
q_cons_vf
type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf
Definition m_phase_change.fpp.f90:1119
k
integer, intent(in) k
Definition m_phase_change.fpp.f90:1121
j
integer, intent(in) j
Definition m_phase_change.fpp.f90:1121
l
integer, intent(in) l
Definition m_phase_change.fpp.f90:1121
m_body_forces
Computes gravitational and user-defined body force source terms for the momentum equations.
Definition m_body_forces.fpp.f90:292
m_body_forces::s_compute_body_forces_rhs
subroutine, public s_compute_body_forces_rhs(q_prim_vf, q_cons_vf, rhs_vf)
This subroutine calculates the source term due to body forces so the system can be advanced in time.
Definition m_body_forces.fpp.f90:546
m_body_forces::s_initialize_body_forces_module
impure subroutine, public s_initialize_body_forces_module
This subroutine initializes the module global array of mixture densities in each grid cell.
Definition m_body_forces.fpp.f90:329
m_body_forces::s_compute_mixture_density
subroutine s_compute_mixture_density(q_cons_vf)
This subroutine calculates the mixture density at each cell center param q_cons_vf Conservative varia...
Definition m_body_forces.fpp.f90:482
m_body_forces::s_compute_acceleration
subroutine s_compute_acceleration(t)
This subroutine computes the acceleration at time t.
Definition m_body_forces.fpp.f90:447
m_body_forces::rhom
real(wp), dimension(:, :, :), allocatable rhom
Definition m_body_forces.fpp.f90:311
m_body_forces::s_finalize_body_forces_module
impure subroutine, public s_finalize_body_forces_module
Deallocates module variables used for body force computations.
Definition m_body_forces.fpp.f90:773
m_derived_types
Shared derived types for field data, patch geometry, bubble dynamics, and MPI I/O structures.
Definition m_derived_types.fpp.f90:292
m_global_parameters
Global parameters for the computational domain, fluid properties, and simulation algorithm configurat...
Definition m_global_parameters.fpp.f90:303
m_global_parameters::buff_size
integer buff_size
The number of cells that are necessary to be able to store enough boundary conditions data to march t...
Definition m_global_parameters.fpp.f90:906
m_global_parameters::p
integer p
Definition m_global_parameters.fpp.f90:331
m_global_parameters::n
integer n
Definition m_global_parameters.fpp.f90:331
m_global_parameters::m
integer m
Definition m_global_parameters.fpp.f90:331
m_nvtx
NVIDIA NVTX profiling API bindings for GPU performance instrumentation.
Definition m_nvtx.f90:6
m_variables_conversion
Conservative-to-primitive variable conversion, mixture property evaluation, and pressure computation.
Definition m_variables_conversion.fpp.f90:303