<p><b>dwj07@fsu.edu</b> 2012-03-18 12:58:09 -0600 (Sun, 18 Mar 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Cleaning up non-dimensionalization of omega.<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 16:49:17 UTC (rev 1666)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin-template.F        2012-03-18 18:58:09 UTC (rev 1667)
@@ -364,7 +364,7 @@
    allocate(psiVertex(nVerticesNew))
    r0 = 140.0
    delta = 0.05
-   omega0 = 2.0 * r0* pi / (2.0 * 60.0 * 60.0) ! one revolution of a circle in 2 hours 
+   omega0 = 2.0 * pi / (2.0 * 60.0 * 60.0) ! one revolution of a circle in 2 hours 
 
    do iVertex = 1, nVerticesNew
        xin = xVertexNew(iVertex) - xmid
@@ -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/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) = (psiVertex(verticesOnEdgeNew(2, iEdge)) - psiVertex(verticesOnEdgeNew(1, iEdge))) / (dvEdgeNew(iEdge))
+     uNew(1, :, iEdge) = - r0 * (psiVertex(verticesOnEdgeNew(2, iEdge)) - psiVertex(verticesOnEdgeNew(1, iEdge))) / (dvEdgeNew(iEdge))
    end do
 
    do iCell = 1, nCellsNew

</font>
</pre>