10 SUBROUTINE ddxdn_28 (ds, a, b)
11 real,
dimension(28,28,28),
intent(in):: a
12 real,
dimension(28,28,28):: b
21 b(ix,iy,iz) = c*(a(ix ,iy,iz) - a(ix-1,iy,iz))
25 END SUBROUTINE ddxdn_28
28 SUBROUTINE ddydn_28 (ds, a, b)
29 real,
dimension(28,28,28),
intent(in):: a
30 real,
dimension(28,28,28):: b
39 b(ix,iy,iz) = c*(a(ix,iy ,iz) - a(ix,iy-1,iz))
43 END SUBROUTINE ddydn_28
46 SUBROUTINE ddzdn_28 (ds, a, b)
47 real,
dimension(28,28,28),
intent(in):: a
48 real,
dimension(28,28,28):: b
57 b(ix,iy,iz) = c*(a(ix,iy,iz ) - a(ix,iy,iz-1))
61 END SUBROUTINE ddzdn_28
64 SUBROUTINE ddxup_28 (ds, a, b)
65 real,
dimension(28,28,28),
intent(in):: a
66 real,
dimension(28,28,28):: b
75 b(ix,iy,iz) = c*(a(ix+1,iy,iz) - a(ix ,iy,iz))
79 END SUBROUTINE ddxup_28
82 SUBROUTINE ddyup_28 (ds, a, b)
83 real,
dimension(28,28,28),
intent(in):: a
84 real,
dimension(28,28,28):: b
93 b(ix,iy,iz) = c*(a(ix,iy+1,iz) - a(ix,iy ,iz))
97 END SUBROUTINE ddyup_28
100 SUBROUTINE ddzup_28 (ds, a, b)
101 real,
dimension(28,28,28),
intent(in):: a
102 real,
dimension(28,28,28):: 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_28
118 SUBROUTINE xdn_28 (a, b)
119 real,
dimension(28,28,28),
intent(in):: a
120 real,
dimension(28,28,28):: 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_28
136 SUBROUTINE ydn_28 (a, b)
137 real,
dimension(28,28,28),
intent(in):: a
138 real,
dimension(28,28,28):: 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_28
154 SUBROUTINE zdn_28 (a, b)
155 real,
dimension(28,28,28),
intent(in):: a
156 real,
dimension(28,28,28):: 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_28
172 SUBROUTINE xup_28 (a, b)
173 real,
dimension(28,28,28),
intent(in):: a
174 real,
dimension(28,28,28):: b
176 integer :: ix, iy, iz
179 b(28,:,:) = a(28,:,:)
183 b(ix,iy,iz) = c*(a(ix+1,iy,iz) + a(ix ,iy,iz))
190 SUBROUTINE yup_28 (a, b)
191 real,
dimension(28,28,28),
intent(in):: a
192 real,
dimension(28,28,28):: b
194 integer :: ix, iy, iz
197 b(:,28,:) = a(:,28,:)
201 b(ix,iy,iz) = c*(a(ix,iy+1,iz) + a(ix,iy ,iz))
208 SUBROUTINE zup_28 (a, b)
209 real,
dimension(28,28,28),
intent(in):: a
210 real,
dimension(28,28,28):: b
212 integer :: ix, iy, iz
215 b(:,:,28) = a(:,:,28)
219 b(ix,iy,iz) = c*(a(ix,iy,iz+1) + a(ix,iy,iz ))