<p><b>dwj07@fsu.edu</b> 2012-03-18 15:38:10 -0600 (Sun, 18 Mar 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Cleaning up non-divergent flow field again.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin-template.F        2012-03-18 18:58:09 UTC (rev 1667)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin-template.F        2012-03-18 21:38:10 UTC (rev 1668)
@@ -372,7 +372,7 @@
 
        r = sqrt(xin**2 + yin**2)/r0 ! nondimensional radial distance
 
-       psiVertex(iVertex) = omega0*r0/4.0*(r**2 - (r + 1)*sqrt((r-1)**2 + delta**2) + delta**2*log((r - 1) + sqrt((r-1)**2 + delta**2)))
+       psiVertex(iVertex) = omega0*r0**2/4.0*(r**2 - (r + 1)*sqrt((r-1)**2 + delta**2) + delta**2*log((r - 1) + sqrt((r-1)**2 + delta**2)))
    end do
 
    do iEdge = 1, nEdgesNew
@@ -396,7 +396,7 @@
 
      r = sqrt( (xEdgeNew(iEdge) - xmid)**2 + (yEdgeNew(iEdge) - ymid)**2)
 
-     uNew(1, :, iEdge) = - r0 * (psiVertex(verticesOnEdgeNew(2, iEdge)) - psiVertex(verticesOnEdgeNew(1, iEdge))) / (dvEdgeNew(iEdge))
+     uNew(1, :, iEdge) = - 1.0 * (psiVertex(verticesOnEdgeNew(2, iEdge)) - psiVertex(verticesOnEdgeNew(1, iEdge))) / (dvEdgeNew(iEdge))
    end do
 
    do iCell = 1, nCellsNew

</font>
</pre>