bulk_flux itermax

[jcwarner]

Small tiling bug in Nonlinear/bulk_flux due to variable itermax. The variable itermax is set in an "i" loop. Itermax takes a value of 1 or 3, depending on the last setting in the "i" loop. After that i loop completes, we have

DO Iter=1,IterMax

DO i=Istr-1,IendR

.....

Well, the number of iterrations for any given i-tile will depend on the last value of itermax. Thus, the values in an i-tile will depend on how it is tiled. Solution: dimension itermax(private_1d_scratch_array) and compute itermax(i). Then change the order of the do loops as:

DO i=Istr-1,IendR

DO Iter=1,IterMax(i)

.....

--- or--- just set itermax = 3 and be done with it (?).