10 SUBROUTINE ddxdn_32 (ds, a, b)
11 real,
dimension(32,32,32),
intent(in):: a
12 real,
dimension(32,32,32):: b
21 b(ix,iy,iz) = c*(a(ix ,iy,iz) - a(ix-1,iy,iz))
25 END SUBROUTINE ddxdn_32
28 SUBROUTINE ddydn_32 (ds, a, b)
29 real,
dimension(32,32,32),
intent(in):: a
30 real,
dimension(32,32,32):: b
39 b(ix,iy,iz) = c*(a(ix,iy ,iz) - a(ix,iy-1,iz))
43 END SUBROUTINE ddydn_32
46 SUBROUTINE ddzdn_32 (ds, a, b)
47 real,
dimension(32,32,32),
intent(in):: a
48 real,
dimension(32,32,32):: b
57 b(ix,iy,iz) = c*(a(ix,iy,iz ) - a(ix,iy,iz-1))
61 END SUBROUTINE ddzdn_32
64 SUBROUTINE ddxup_32 (ds, a, b)
65 real,
dimension(32,32,32),
intent(in):: a
66 real,
dimension(32,32,32):: b
75 b(ix,iy,iz) = c*(a(ix+1,iy,iz) - a(ix ,iy,iz))
79 END SUBROUTINE ddxup_32
82 SUBROUTINE ddyup_32 (ds, a, b)
83 real,
dimension(32,32,32),
intent(in):: a
84 real,
dimension(32,32,32):: b
93 b(ix,iy,iz) = c*(a(ix,iy+1,iz) - a(ix,iy ,iz))
97 END SUBROUTINE ddyup_32
100 SUBROUTINE ddzup_32 (ds, a, b)
101 real,
dimension(32,32,32),
intent(in):: a
102 real,
dimension(32,32,32):: 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_32
118 SUBROUTINE xdn_32 (a, b)
119 real,
dimension(32,32,32),
intent(in):: a
120 real,
dimension(32,32,32):: 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_32
136 SUBROUTINE ydn_32 (a, b)
137 real,
dimension(32,32,32),
intent(in):: a
138 real,
dimension(32,32,32):: 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_32
154 SUBROUTINE zdn_32 (a, b)
155 real,
dimension(32,32,32),
intent(in):: a
156 real,
dimension(32,32,32):: 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_32
172 SUBROUTINE xup_32 (a, b)
173 real,
dimension(32,32,32),
intent(in):: a
174 real,
dimension(32,32,32):: b
176 integer :: ix, iy, iz
179 b(32,:,:) = a(32,:,:)
183 b(ix,iy,iz) = c*(a(ix+1,iy,iz) + a(ix ,iy,iz))
190 SUBROUTINE yup_32 (a, b)
191 real,
dimension(32,32,32),
intent(in):: a
192 real,
dimension(32,32,32):: b
194 integer :: ix, iy, iz
197 b(:,32,:) = a(:,32,:)
201 b(ix,iy,iz) = c*(a(ix,iy+1,iz) + a(ix,iy ,iz))
208 SUBROUTINE zup_32 (a, b)
209 real,
dimension(32,32,32),
intent(in):: a
210 real,
dimension(32,32,32):: b
212 integer :: ix, iy, iz
215 b(:,:,32) = a(:,:,32)
219 b(ix,iy,iz) = c*(a(ix,iy,iz+1) + a(ix,iy,iz ))