<p><b>duda</b> 2009-09-03 14:52:12 -0600 (Thu, 03 Sep 2009)</p><p>Correct bug in scalar advection: h_edge should come from provisional state, s,<br>
rather than tend.<br>
<br>
M    module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: trunk/swmodel/src/module_time_integration.F
===================================================================
--- trunk/swmodel/src/module_time_integration.F        2009-09-02 21:13:33 UTC (rev 43)
+++ trunk/swmodel/src/module_time_integration.F        2009-09-03 20:52:12 UTC (rev 44)
@@ -481,7 +481,7 @@
                do k=1,grid % nVertLevels
                   do iTracer=1,grid % nTracers
                      tracer_edge = 0.5 * (s % tracers % array(iTracer,k,cell1) + s % tracers % array(iTracer,k,cell2))
-                     flux = s % u % array(k,iEdge) * grid % dvEdge % array(iEdge) * tend % h_edge % array(k,iEdge) * tracer_edge
+                     flux = s % u % array(k,iEdge) * grid % dvEdge % array(iEdge) * s % h_edge % array(k,iEdge) * tracer_edge
                      tend % tracers % array(iTracer,k,cell1) = tend % tracers % array(iTracer,k,cell1) - flux
                      tend % tracers % array(iTracer,k,cell2) = tend % tracers % array(iTracer,k,cell2) + flux
                   end do 

</font>
</pre>