10 SUBROUTINE ddxdn_16 (ds, a, b)
11 real,
dimension(16,16,16),
intent(in):: a
12 real,
dimension(16,16,16):: b
21 b(ix,iy,iz) = c*(a(ix ,iy,iz) - a(ix-1,iy,iz))
25 END SUBROUTINE ddxdn_16
28 SUBROUTINE ddydn_16 (ds, a, b)
29 real,
dimension(16,16,16),
intent(in):: a
30 real,
dimension(16,16,16):: b
39 b(ix,iy,iz) = c*(a(ix,iy ,iz) - a(ix,iy-1,iz))
43 END SUBROUTINE ddydn_16
46 SUBROUTINE ddzdn_16 (ds, a, b)
47 real,
dimension(16,16,16),
intent(in):: a
48 real,
dimension(16,16,16):: b
57 b(ix,iy,iz) = c*(a(ix,iy,iz ) - a(ix,iy,iz-1))
61 END SUBROUTINE ddzdn_16
64 SUBROUTINE ddxup_16 (ds, a, b)
65 real,
dimension(16,16,16),
intent(in):: a
66 real,
dimension(16,16,16):: b
75 b(ix,iy,iz) = c*(a(ix+1,iy,iz) - a(ix ,iy,iz))
79 END SUBROUTINE ddxup_16
82 SUBROUTINE ddyup_16 (ds, a, b)
83 real,
dimension(16,16,16),
intent(in):: a
84 real,
dimension(16,16,16):: b
93 b(ix,iy,iz) = c*(a(ix,iy+1,iz) - a(ix,iy ,iz))
97 END SUBROUTINE ddyup_16
100 SUBROUTINE ddzup_16 (ds, a, b)
101 real,
dimension(16,16,16),
intent(in):: a
102 real,
dimension(16,16,16):: b
103 integer :: ix, iy, iz
111 b(ix,iy,iz) = c*(a(ix,iy,iz+1) - a(ix,iy,iz ))
115 END SUBROUTINE ddzup_16
118 SUBROUTINE xdn_16 (a, b)
119 real,
dimension(16,16,16),
intent(in):: a
120 real,
dimension(16,16,16):: b
122 integer :: ix, iy, iz
129 b(ix,iy,iz) = c*(a(ix ,iy,iz) + a(ix-1,iy,iz))
133 END SUBROUTINE xdn_16
136 SUBROUTINE ydn_16 (a, b)
137 real,
dimension(16,16,16),
intent(in):: a
138 real,
dimension(16,16,16):: b
140 integer :: ix, iy, iz
147 b(ix,iy,iz) = c*(a(ix,iy,iz) + a(ix,iy-1,iz))
151 END SUBROUTINE ydn_16
154 SUBROUTINE zdn_16 (a, b)
155 real,
dimension(16,16,16),
intent(in):: a
156 real,
dimension(16,16,16):: b
158 integer :: ix, iy, iz
165 b(ix,iy,iz) = c*(a(ix,iy,iz) + a(ix,iy,iz-1))
169 END SUBROUTINE zdn_16
172 SUBROUTINE xup_16 (a, b)
173 real,
dimension(16,16,16),
intent(in):: a
174 real,
dimension(16,16,16):: b
176 integer :: ix, iy, iz
179 b(16,:,:) = a(16,:,:)
183 b(ix,iy,iz) = c*(a(ix+1,iy,iz) + a(ix ,iy,iz))
190 SUBROUTINE yup_16 (a, b)
191 real,
dimension(16,16,16),
intent(in):: a
192 real,
dimension(16,16,16):: b
194 integer :: ix, iy, iz
197 b(:,16,:) = a(:,16,:)
201 b(ix,iy,iz) = c*(a(ix,iy+1,iz) + a(ix,iy ,iz))
208 SUBROUTINE zup_16 (a, b)
209 real,
dimension(16,16,16),
intent(in):: a
210 real,
dimension(16,16,16):: b
212 integer :: ix, iy, iz
215 b(:,:,16) = a(:,:,16)
219 b(ix,iy,iz) = c*(a(ix,iy,iz+1) + a(ix,iy,iz ))