10 SUBROUTINE ddxdn_36 (ds, a, b)
11 real,
dimension(36,36,36),
intent(in):: a
12 real,
dimension(36,36,36):: b
20 b(ix,iy,iz) = c*(a(ix ,iy,iz) - a(ix-1,iy,iz))
25 END SUBROUTINE ddxdn_36
28 SUBROUTINE ddydn_36 (ds, a, b)
29 real,
dimension(36,36,36),
intent(in):: a
30 real,
dimension(36,36,36):: b
39 b(ix,iy,iz) = c*(a(ix,iy ,iz) - a(ix,iy-1,iz))
43 END SUBROUTINE ddydn_36
46 SUBROUTINE ddzdn_36 (ds, a, b)
47 real,
dimension(36,36,36),
intent(in):: a
48 real,
dimension(36,36,36):: b
57 b(ix,iy,iz) = c*(a(ix,iy,iz ) - a(ix,iy,iz-1))
61 END SUBROUTINE ddzdn_36
64 SUBROUTINE ddxup_36 (ds, a, b)
65 real,
dimension(36,36,36),
intent(in):: a
66 real,
dimension(36,36,36):: b
74 b(ix,iy,iz) = c*(a(ix+1,iy,iz) - a(ix ,iy,iz))
79 END SUBROUTINE ddxup_36
82 SUBROUTINE ddyup_36 (ds, a, b)
83 real,
dimension(36,36,36),
intent(in):: a
84 real,
dimension(36,36,36):: b
93 b(ix,iy,iz) = c*(a(ix,iy+1,iz) - a(ix,iy ,iz))
97 END SUBROUTINE ddyup_36
100 SUBROUTINE ddzup_36 (ds, a, b)
101 real,
dimension(36,36,36),
intent(in):: a
102 real,
dimension(36,36,36):: 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_36
118 SUBROUTINE xdn_36 (a, b)
119 real,
dimension(36,36,36),
intent(in):: a
120 real,
dimension(36,36,36):: 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_36
136 SUBROUTINE ydn_36 (a, b)
137 real,
dimension(36,36,36),
intent(in):: a
138 real,
dimension(36,36,36):: 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_36
154 SUBROUTINE zdn_36 (a, b)
155 real,
dimension(36,36,36),
intent(in):: a
156 real,
dimension(36,36,36):: 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_36
172 SUBROUTINE xup_36 (a, b)
173 real,
dimension(36,36,36),
intent(in):: a
174 real,
dimension(36,36,36):: b
176 integer :: ix, iy, iz
179 b(36,:,:) = a(36,:,:)
183 b(ix,iy,iz) = c*(a(ix+1,iy,iz) + a(ix ,iy,iz))
190 SUBROUTINE yup_36 (a, b)
191 real,
dimension(36,36,36),
intent(in):: a
192 real,
dimension(36,36,36):: b
194 integer :: ix, iy, iz
197 b(:,36,:) = a(36,:,:)
201 b(ix,iy,iz) = c*(a(ix,iy+1,iz) + a(ix,iy ,iz))
208 SUBROUTINE zup_36 (a, b)
209 real,
dimension(36,36,36),
intent(in):: a
210 real,
dimension(36,36,36):: b
212 integer :: ix, iy, iz
215 b(:,:,36) = a(:,:,36)
219 b(ix,iy,iz) = c*(a(ix,iy,iz+1) + a(ix,iy,iz ))