1 | do k = 1, nz |
2 | do j = 1, ny |
3 | do i = 1, nx |
4 | do m = 1, nc |
5 | do mp = 1, nc |
6 | a(i,j,k,m) = u(mp,m) * (b(i+1,j,k,mp) + b(i-1,j,k,mp) & |
7 | b(i,j+1,k,mp) + b(i,j-1,k,mp) + b(i,j,k+1,mp) + & |
8 | b(i,j,k-1,mp)) + a(i,j,k,m) |
9 | enddo |
10 | enddo |
11 | enddo |
12 | enddo |
13 | enddo |