<p><b>dwj07@fsu.edu</b> 2012-02-27 11:13:59 -0700 (Mon, 27 Feb 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Changing the way epsilon is used to remove NaNs in monotonic advection.<br>
<br>
        This improves compatibility with some versions of intel's compilers.<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-24 03:14:36 UTC (rev 1534)
+++ branches/ocean_projects/monotonic_advection/src/operators/mpas_tracer_advection_mono.F        2012-02-27 18:13:59 UTC (rev 1535)
@@ -224,10 +224,10 @@
             tracer_max_new = (tracer_cur(k,iCell)*h(k,iCell) + dt*(upwind_tendency(k,iCell)+flux_incoming(k,iCell))) * inv_h_new(k,iCell)
             tracer_upwind_new = (tracer_cur(k,iCell)*h(k,iCell) + dt*upwind_tendency(k,iCell)) * inv_h_new(k,iCell)
            
-            scale_factor = (tracer_max(k,iCell)-tracer_upwind_new)/(tracer_max_new-tracer_upwind_new+eps)
+            scale_factor = (tracer_max(k,iCell)-tracer_upwind_new)/max(tracer_max_new-tracer_upwind_new,eps)
             flux_incoming(k,iCell) = min( 1.0, max( 0.0, scale_factor) )
 
-            scale_factor = (tracer_upwind_new-tracer_min(k,iCell))/(tracer_upwind_new-tracer_min_new+eps)
+            scale_factor = (tracer_upwind_new-tracer_min(k,iCell))/max(tracer_upwind_new-tracer_min_new,eps)
             flux_outgoing(k,iCell) = min( 1.0, max( 0.0, scale_factor) )
           end do ! k loop
         end do ! iCell loop

</font>
</pre>