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
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
! Figure out the extent of the physical part of a chunk of a grid.
! $Header$
! nx, ny and nz contains the full size of the chunk.
! nx, ny and nz are declared in EHFinder_mod.F90.
nx = cctk_lsh(1)
ny = cctk_lsh(2)
nz = cctk_lsh(3)
! ngx, ngy, ngz contains the size of any ghostzones.
! ngx, ngy and ngz are declared in EHFinder_mod.F90.
ngx = cctk_nghostzones(1)
ngy = cctk_nghostzones(2)
ngz = cctk_nghostzones(3)
! ixl and ixr is set to the minimum and maximum local cell index in the
! x-direction. The corresponding is done for the y- and z-directions.
! ixl, ixr, jyl, jyr, kzl, kzr are declared in EHFinder_mod.F90.
symx = .false.; symy = .false.; symz = .false.;
ixl = 1
ixr = nx
jyl = 1
jyr = ny
kzl = 1
kzr = nz
! I check to see if we have any processor boundaries. If we have, adjust the
! minimum and maximum cell index appropriately with the size of the ghostzone.
if ( cctk_bbox(1) .eq. 0 ) ixl = ixl + ngx
if ( cctk_bbox(2) .eq. 0 ) ixr = ixr - ngx
if ( cctk_bbox(3) .eq. 0 ) jyl = jyl + ngy
if ( cctk_bbox(4) .eq. 0 ) jyr = jyr - ngy
if ( cctk_bbox(5) .eq. 0 ) kzl = kzl + ngz
if ( cctk_bbox(6) .eq. 0 ) kzr = kzr - ngz
! In octant I check to see if the chunk of grid contains any symmetry
! boundaries. This is done by checking if the lowest index of the local grid
! as seen on the global grid is equal to zero. If this is the case adjust the
! minimum and maximum cell index appropriately with the size of the ghostzone.
! Also symx, symy and symz are set to true.
! symx, symy and symz are declared in EHFinder_mod.F90.
if ( CCTK_EQUALS ( domain, 'octant' ) ) then
if (cctk_lbnd(1) .eq. 0 ) then
ixl = ixl + ngx
symx = .true.
end if
if (cctk_lbnd(2) .eq. 0 ) then
jyl = jyl + ngy
symy = .true.
end if
if (cctk_lbnd(3) .eq. 0 ) then
kzl = kzl + ngz
symz = .true.
end if
end if
! In quadrant I do the same, except here I have to check in which direction
! the quadrant is directed.
if ( CCTK_EQUALS ( domain, 'quadrant' ) ) then
if ( CCTK_EQUALS ( quadrant_direction, 'x' ) ) then
if (cctk_lbnd(2) .eq. 0 ) then
jyl = jyl + ngy
symy = .true.
end if
if (cctk_lbnd(3) .eq. 0 ) then
kzl = kzl + ngz
symz = .true.
end if
end if
if ( CCTK_EQUALS ( quadrant_direction, 'y' ) ) then
if (cctk_lbnd(1) .eq. 0 ) then
ixl = ixl + ngx
symx = .true.
end if
if (cctk_lbnd(3) .eq. 0 ) then
kzl = kzl + ngz
symz = .true.
end if
end if
if ( CCTK_EQUALS ( quadrant_direction,'z' ) ) then
if (cctk_lbnd(1) .eq. 0 ) then
ixl = ixl + ngx
symx = .true.
end if
if (cctk_lbnd(2) .eq. 0 ) then
jyl = jyl + ngy
symy = .true.
end if
end if
end if
! Ditto for bitant. Here I have to check which plane is the symmetry-plane.
if ( CCTK_EQUALS ( domain, 'bitant' ) ) then
if ( CCTK_EQUALS ( bitant_plane, 'xy' ) ) then
if (cctk_lbnd(3) .eq. 0 ) then
kzl = kzl + ngz
symz = .true.
end if
endif
if ( CCTK_EQUALS ( bitant_plane, 'xz' ) ) then
if (cctk_lbnd(2) .eq. 0 ) then
jyl = jyl + ngy
symy = .true.
end if
endif
if ( CCTK_EQUALS ( bitant_plane, 'yz' ) ) then
if (cctk_lbnd(1) .eq. 0 ) then
ixl = ixl + ngx
symx = .true.
end if
endif
end if
! At the end ixl, ixr, jyl, jyr, kzl, kzr should contain the minimum index
! and maximum index in each direction that are not ghost cells or
! symmetry cells.
|
