<p><b>dwj07@fsu.edu</b> 2012-02-03 11:30:43 -0700 (Fri, 03 Feb 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Working version of monotonic vertical advection.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/advection_routines/src/operators/mpas_tracer_advection_mono.F
===================================================================
--- branches/ocean_projects/advection_routines/src/operators/mpas_tracer_advection_mono.F        2012-02-03 18:23:05 UTC (rev 1462)
+++ branches/ocean_projects/advection_routines/src/operators/mpas_tracer_advection_mono.F        2012-02-03 18:30:43 UTC (rev 1463)
@@ -125,26 +125,26 @@
           ! zero flux at top and bottom
            k = 1
            vert_flux(k,iCell) = 0.
-!          s_max(k,iCell) = max(tracer_cur(k,iCell),tracer_cur(k+1,iCell))
-!          s_min(k,iCell) = min(tracer_cur(k,iCell),tracer_cur(k+1,iCell))
+           s_max(k,iCell) = max(tracer_cur(k,iCell),tracer_cur(k+1,iCell))
+           s_min(k,iCell) = min(tracer_cur(k,iCell),tracer_cur(k+1,iCell))
 
-!          k = 2
-!          vert_flux(k,iCell) = w(k,iCell)*(zWeightK(k)*tracer_cur(k,iCell)+zWeightKm1(k)*tracer_cur(k-1,iCell))
-!          s_max(k,iCell) = max(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
-!          s_min(k,iCell) = min(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
-!            
-!          do k=3,maxLevelCell(iCell)-1
-!             vert_flux(k,iCell) = flux3( tracer_cur(k-2,iCell),tracer_cur(k-1,iCell),  &amp;
-!                                    tracer_cur(k  ,iCell),tracer_cur(k+1,iCell),  &amp;
-!                                    w(k,iCell), coef_3rd_order )
-!             s_max(k,iCell) = max(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
-!             s_min(k,iCell) = min(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
-!          end do 
-!
-!          k = maxLevelCell(iCell)
-!          vert_flux(k,iCell) = w(k,iCell)*(zWeightK(k)*tracer_cur(k,iCell)+zWeightKm1(k)*tracer_cur(k-1,iCell))
-!          s_max(k,iCell) = max(tracer_cur(k,iCell),tracer_cur(k-1,iCell))
-!          s_min(k,iCell) = min(tracer_cur(k,iCell),tracer_cur(k-1,iCell))
+           k = 2
+           vert_flux(k,iCell) = w(k,iCell)*(zWeightK(k)*tracer_cur(k,iCell)+zWeightKm1(k)*tracer_cur(k-1,iCell))
+           s_max(k,iCell) = max(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
+           s_min(k,iCell) = min(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
+             
+           do k=3,maxLevelCell(iCell)-1
+              vert_flux(k,iCell) = flux3( tracer_cur(k-2,iCell),tracer_cur(k-1,iCell),  &amp;
+                                     tracer_cur(k  ,iCell),tracer_cur(k+1,iCell),  &amp;
+                                     w(k,iCell), coef_3rd_order )
+              s_max(k,iCell) = max(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
+              s_min(k,iCell) = min(tracer_cur(k-1,iCell),tracer_cur(k,iCell),tracer_cur(k+1,iCell))
+           end do 
+ 
+           k = maxLevelCell(iCell)
+           vert_flux(k,iCell) = w(k,iCell)*(zWeightK(k)*tracer_cur(k,iCell)+zWeightKm1(k)*tracer_cur(k-1,iCell))
+           s_max(k,iCell) = max(tracer_cur(k,iCell),tracer_cur(k-1,iCell))
+           s_min(k,iCell) = min(tracer_cur(k,iCell),tracer_cur(k-1,iCell))
 
            vert_flux(maxLevelCell(iCell)+1,iCell) = 0.
 

</font>
</pre>