2 c---------------------------------------------------------------------
3 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
8 c---------------------------------------------------------------------
10 c---------------------------------------------------------------------
11 c addition of update to the vector u
12 c---------------------------------------------------------------------
20 do k = start(3,c), cell_size(3,c)-end(3,c)-1
21 do j = start(2,c), cell_size(2,c)-end(2,c)-1
22 do i = start(1,c), cell_size(1,c)-end(1,c)-1
23 u(i,j,k,m,c) = u(i,j,k,m,c) + rhs(i,j,k,m,c)