-
Notifications
You must be signed in to change notification settings - Fork 0
/
field_summary_kernel.f90
68 lines (57 loc) · 2.34 KB
/
field_summary_kernel.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
!Crown Copyright 2014 AWE.
!
! This file is part of TeaLeaf.
!
! TeaLeaf is free software: you can redistribute it and/or modify it under
! the terms of the GNU General Public License as published by the
! Free Software Foundation, either version 3 of the License, or (at your option)
! any later version.
!
! TeaLeaf is distributed in the hope that it will be useful, but
! WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
! FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
! details.
!
! You should have received a copy of the GNU General Public License along with
! TeaLeaf. If not, see http://www.gnu.org/licenses/.
!> @brief Fortran field summary kernel
!> @author David Beckingsale, Wayne Gaudin
!> @details The total mass, internal energy, temperature are calculated
MODULE field_summary_kernel_module
CONTAINS
SUBROUTINE field_summary_kernel(x_min,x_max,y_min,y_max,z_min,z_max,halo_exchange_depth,&
volume, &
density, &
energy1, &
u,vol,mass,ie,temp)
IMPLICIT NONE
INTEGER :: x_min,x_max,y_min,y_max,z_min,z_max,halo_exchange_depth
REAL(KIND=8), DIMENSION(x_min-2:x_max+2,y_min-2:y_max+2,z_min-2:z_max+2) :: volume
REAL(KIND=8), DIMENSION(x_min-halo_exchange_depth:x_max+halo_exchange_depth,&
y_min-halo_exchange_depth:y_max+halo_exchange_depth,&
z_min-halo_exchange_depth:z_max+halo_exchange_depth) :: density,energy1,u
REAL(KIND=8) :: vol,mass,ie,temp
REAL(KIND=8) :: cell_vol,cell_mass
INTEGER :: j,k,l
vol=0.0
mass=0.0
ie=0.0
temp=0.0
!$OMP PARALLEL PRIVATE(cell_vol,cell_mass) REDUCTION(+ : vol,mass,ie,temp)
!$OMP DO
DO l=z_min,z_max
DO k=y_min,y_max
DO j=x_min,x_max
cell_vol=volume(j,k,l)
cell_mass=cell_vol*density(j,k,l)
vol=vol+cell_vol
mass=mass+cell_mass
ie=ie+cell_mass*energy1(j,k,l)
temp = temp + cell_mass*u(j,k,l)
ENDDO
ENDDO
ENDDO
!$OMP END DO
!$OMP END PARALLEL
END SUBROUTINE field_summary_kernel
END MODULE field_summary_kernel_module