<p><b>dwj07@fsu.edu</b> 2012-02-21 21:12:00 -0700 (Tue, 21 Feb 2012)</p><p><br>
-- BRANCH COMMIT --<br>
<br>
Fixing bounds check, for multiple processors.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/monotonic_advection/src/operators/mpas_tracer_advection_mono.F
===================================================================
--- branches/ocean_projects/monotonic_advection/src/operators/mpas_tracer_advection_mono.F 2012-02-21 22:49:18 UTC (rev 1526)
+++ branches/ocean_projects/monotonic_advection/src/operators/mpas_tracer_advection_mono.F 2012-02-22 04:12:00 UTC (rev 1527)
@@ -288,17 +288,17 @@
end do ! iCell loop
!dwj 02/21/12 build min and max of cur and new tracer states to check bounds
- cur_min = minval(tracer_cur)
- cur_max = maxval(tracer_cur)
- new_min = minval(tracer_new)
- new_max = maxval(tracer_new)
+ cur_min = minval(tracer_cur(:,1:nCellsSolve))
+ cur_max = maxval(tracer_cur(:,1:nCellsSolve))
+ new_min = minval(tracer_new(:,1:nCellsSolve))
+ new_max = maxval(tracer_new(:,1:nCellsSolve))
!dwj 02/21/12 do bounds check
-! if(new_min < cur_min ) then
+! if(new_min < cur_min-eps) then
! write(*,*) 'Minimum out of bounds on tracer ', iTracer, cur_min, new_min
! end if
-! if(new_max > cur_max) then
+! if(new_max > cur_max+eps) then
! write(*,*) 'Maximum out of bounds on tracer ', iTracer, cur_max, new_max
! end if
</font>
</pre>