<p><b>duda</b> 2010-10-14 17:43:48 -0600 (Thu, 14 Oct 2010)</p><p>BRANCH COMMIT<br>
<br>
Fix if-test in advance_scalars_mono() so that we correctly revert <br>
to second-order vertical advection for scalars at k=nVertLevels-1.<br>
<br>
<br>
M    src/core_nhyd_atmos/module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F        2010-10-14 21:48:16 UTC (rev 557)
+++ branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F        2010-10-14 23:43:48 UTC (rev 558)
@@ -1562,7 +1562,7 @@
 
             if (k &lt; grid % nVertLevels) then
 
-               if ((config_scalar_vadv_order == 2) .or. (k==1) .or. (k==grid % nVertLevels)) then
+               if ((config_scalar_vadv_order == 2) .or. (k==1) .or. (k==grid % nVertLevels-1)) then
                   do iScalar=1,s_old % num_scalars
                      v_flux(iScalar,iCell,km0) = dt * wwAvg(k+1,iCell) *   &amp;
                           (fnm(k+1) * scalar_new(iScalar,k+1,iCell) + fnp(k+1) * scalar_new(iScalar,k,iCell))

</font>
</pre>