<p><b>mpetersen@lanl.gov</b> 2012-05-02 11:21:24 -0600 (Wed, 02 May 2012)</p><p>BRANCH COMMIT, imp_vert_mix_error.  Correcting the error: Coefficient A in tridiagonal solve had divisor h(k-1), should be h(k).  This does not produce a bit-for-bit match, as expected.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/imp_vert_mix_error/src/core_ocean/mpas_ocn_vmix.F
===================================================================
--- branches/ocean_projects/imp_vert_mix_error/src/core_ocean/mpas_ocn_vmix.F        2012-05-02 17:09:00 UTC (rev 1859)
+++ branches/ocean_projects/imp_vert_mix_error/src/core_ocean/mpas_ocn_vmix.F        2012-05-02 17:21:24 UTC (rev 1860)
@@ -325,7 +325,7 @@
          do k=2,N
             A(k) = -2.0*dt*vertViscTopOfEdge(k,iEdge) &amp;
                / (h_edge(k-1,iEdge) + h_edge(k,iEdge)) &amp;
-               / h_edge(k-1,iEdge)
+               / h_edge(k,iEdge)
          enddo
 
          ! C is upper diagonal term
@@ -546,7 +546,7 @@
          A(1)=0
          do k=2,N
             A(k) = -2.0*dt*vertDiffTopOfCell(k,iCell) &amp;
-                 / (h(k-1,iCell) + h(k,iCell)) / h(k-1,iCell)
+                 / (h(k-1,iCell) + h(k,iCell)) / h(k,iCell)
          enddo
 
          ! C is upper diagonal term

</font>
</pre>