1 c---------------------------------------------------------------------
2 c---------------------------------------------------------------------
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
9 c---------------------------------------------------------------------
10 c addition of update to the vector u
11 c---------------------------------------------------------------------
18 do k = start(3,c), cell_size(3,c)-end(3,c)-1
19 do j = start(2,c), cell_size(2,c)-end(2,c)-1
20 do i = start(1,c), cell_size(1,c)-end(1,c)-1
22 u(m,i,j,k,c) = u(m,i,j,k,c) + rhs(m,i,j,k,c)