<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 < 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) * &
(fnm(k+1) * scalar_new(iScalar,k+1,iCell) + fnp(k+1) * scalar_new(iScalar,k,iCell))
</font>
</pre>