-
Notifications
You must be signed in to change notification settings - Fork 0
/
initialise_chunk.f90
70 lines (57 loc) · 2.81 KB
/
initialise_chunk.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
69
70
!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 Driver for chunk initialisation.
!> @author David Beckingsale, Wayne Gaudin
!> @details Invokes the user specified chunk initialisation kernel.
SUBROUTINE initialise_chunk(chunk)
USE tea_module
USE initialise_chunk_kernel_module
IMPLICIT NONE
INTEGER :: chunk
REAL(KIND=8) :: xmin,ymin,zmin,dx,dy,dz
dx=(grid%xmax-grid%xmin)/REAL(grid%x_cells)
dy=(grid%ymax-grid%ymin)/REAL(grid%y_cells)
dz=(grid%zmax-grid%zmin)/REAL(grid%z_cells)
xmin=grid%xmin+dx*REAL(chunks(chunk)%field%left-1)
ymin=grid%ymin+dy*REAL(chunks(chunk)%field%bottom-1)
zmin=grid%zmin+dz*REAL(chunks(chunk)%field%back-1)
IF(use_fortran_kernels) THEN
CALL initialise_chunk_kernel(chunks(chunk)%field%x_min, &
chunks(chunk)%field%x_max, &
chunks(chunk)%field%y_min, &
chunks(chunk)%field%y_max, &
chunks(chunk)%field%z_min, &
chunks(chunk)%field%z_max, &
xmin,ymin,zmin,dx,dy,dz, &
chunks(chunk)%field%vertexx, &
chunks(chunk)%field%vertexdx, &
chunks(chunk)%field%vertexy, &
chunks(chunk)%field%vertexdy, &
chunks(chunk)%field%vertexz, &
chunks(chunk)%field%vertexdz, &
chunks(chunk)%field%cellx, &
chunks(chunk)%field%celldx, &
chunks(chunk)%field%celly, &
chunks(chunk)%field%celldy, &
chunks(chunk)%field%cellz, &
chunks(chunk)%field%celldz, &
chunks(chunk)%field%volume, &
chunks(chunk)%field%xarea, &
chunks(chunk)%field%yarea, &
chunks(chunk)%field%zarea )
ENDIF
END SUBROUTINE initialise_chunk