<p><b>dwj07@fsu.edu</b> 2012-02-07 15:21:26 -0700 (Tue, 07 Feb 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Fixing thickness definition.<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-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin-template.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -27,11 +27,11 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 implicit none
 
-integer, parameter :: nx = 300
+integer, parameter :: nx = 225
 ! This needs to be changed for correct periodic boundaries
 ! Lx is the TOTAL domain width, and needs to be exact for correct periodic
 ! boundaries in x.
-real, parameter :: Lx = 600.0
+real, parameter :: Lx = 450.0
 
 ! original grid variables
 integer :: time, nCells, nEdges, nVertices
@@ -319,8 +319,8 @@
  ! setting for three levels - Set h values for isopycnal system
    write(6,*) ' setting three levels for isopycnal system'
    do i = 1, nVertLevelsMOD
-       hNew(1,i,:) = i * h_total_max / nVertLevelsMOD
-       hZLevel(i) =  i * h_total_max / nVertLevelsMOD
+       hNew(1,i,:) = h_total_max / nVertLevelsMOD
+       hZLevel(i) = h_total_max / nVertLevelsMOD
    end do
 
    h_sNew(:) = -h_total_max
@@ -385,8 +385,8 @@
    ymax = maxval(yEdgeNew)
    ymid = (ymax + ymin)/2.0
 
-   u0 = 2.0 * pi * 75.0 / (2.0 * 60.0 * 60.0) ! one revolution of a circle of radius 75m in 2 hours
-   u0 = 1.0 / 1000.0
+   u0 = 2.0 * pi / (2.0 * 60.0 * 60.0) ! one revolution of a circle in 2 hours
+!  u0 = 1.0 / 1000.0
 !  u0 = 1.0
 
    do iEdge = 1, nEdgesNew
@@ -413,8 +413,10 @@
 
      if(r .le. 140.0) then
          uNew(1,:,iEdge) = dotProd * u0 * r
-     else if(r .le. 280.0) then
-         uNew(1,:,iEdge) = dotProd * u0 * (280.0 - r)
+!    else if(r .le. 280.0) then
+!        uNew(1,:,iEdge) = dotProd * u0 * (140.0 - r)
+     else if(r .le. 210.0) then
+         uNew(1,:,iEdge) = 2.0 * dotProd * u0 * (210.0 - r)
      else
          uNew(1,:,iEdge) = 0.0
      end if

Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin.F        2012-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/advective_transport/basin_src/basin.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -27,11 +27,11 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 implicit none
 
-integer, parameter :: nx = 300
+integer, parameter :: nx = 225
 ! This needs to be changed for correct periodic boundaries
 ! Lx is the TOTAL domain width, and needs to be exact for correct periodic
 ! boundaries in x.
-real, parameter :: Lx = 600.0
+real, parameter :: Lx = 450.0
 
 ! original grid variables
 integer :: time, nCells, nEdges, nVertices
@@ -61,7 +61,7 @@
 real, dimension(40) :: dz
 
 ! Step 1: Set the number of Vertical levels
-integer, parameter :: nVertLevelsMOD = 1
+integer, parameter :: nVertLevelsMOD = 4
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! basin-mod
@@ -385,8 +385,8 @@
    ymax = maxval(yEdgeNew)
    ymid = (ymax + ymin)/2.0
 
-   u0 = 2.0 * pi * 75.0 / (2.0 * 60.0 * 60.0) ! one revolution of a circle of radius 75m in 2 hours
-   u0 = 1.0 / 1000.0
+   u0 = 2.0 * pi / (2.0 * 60.0 * 60.0) ! one revolution of a circle in 2 hours
+!  u0 = 1.0 / 1000.0
 !  u0 = 1.0
 
    do iEdge = 1, nEdgesNew
@@ -413,8 +413,10 @@
 
      if(r .le. 140.0) then
          uNew(1,:,iEdge) = dotProd * u0 * r
-     else if(r .le. 280.0) then
-         uNew(1,:,iEdge) = dotProd * u0 * (280.0 - r)
+!    else if(r .le. 280.0) then
+!        uNew(1,:,iEdge) = dotProd * u0 * (140.0 - r)
+     else if(r .le. 210.0) then
+         uNew(1,:,iEdge) = 2.0 * dotProd * u0 * (210.0 - r)
      else
          uNew(1,:,iEdge) = 0.0
      end if

Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/basin_src/basin-template.F        2012-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/basin_src/basin-template.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -322,8 +322,8 @@
  ! setting for three levels - Set h values for isopycnal system
    write(6,*) ' setting three levels for isopycnal system'
    do i = 1, nVertLevelsMOD
-       hNew(1,i,:) = i * h_total_max / nVertLevelsMOD
-       hZLevel(i) =  i * h_total_max / nVertLevelsMOD
+       hNew(1,i,:) = h_total_max / nVertLevelsMOD
+       hZLevel(i) =  h_total_max / nVertLevelsMOD
    end do
 
    h_sNew(:) = -h_total_max

Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/basin_src/basin-template.F        2012-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/basin_src/basin-template.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -321,8 +321,8 @@
  ! setting for three levels - Set h values for isopycnal system
    write(6,*) ' setting three levels for isopycnal system'
    do i = 1, nVertLevelsMOD
-       hNew(1,i,:) = i * h_total_max / nVertLevelsMOD
-       hZLevel(i) =  i * h_total_max / nVertLevelsMOD
+       hNew(1,i,:) = h_total_max / nVertLevelsMOD
+       hZLevel(i) =  h_total_max / nVertLevelsMOD
    end do
 
    h_sNew(:) = -h_total_max

Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/internal_waves/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/internal_waves/basin_src/basin-template.F        2012-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/internal_waves/basin_src/basin-template.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -322,8 +322,8 @@
  ! setting for three levels - Set h values for isopycnal system
    write(6,*) ' setting three levels for isopycnal system'
    do i = 1, nVertLevelsMOD
-       hNew(1,i,:) = i * h_total_max / nVertLevelsMOD
-       hZLevel(i) =  i * h_total_max / nVertLevelsMOD
+       hNew(1,i,:) = h_total_max / nVertLevelsMOD
+       hZLevel(i) =  h_total_max / nVertLevelsMOD
    end do
 
    h_sNew(:) = -h_total_max

Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/lock_exchange/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/lock_exchange/basin_src/basin-template.F        2012-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/lock_exchange/basin_src/basin-template.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -316,8 +316,8 @@
  ! setting for three levels - Set h values for isopycnal system
    write(6,*) ' setting three levels for isopycnal system'
    do i = 1, nVertLevelsMOD
-       hNew(1,i,:) = i * h_total_max / nVertLevelsMOD
-       hZLevel(i) =  i * h_total_max / nVertLevelsMOD
+       hNew(1,i,:) = h_total_max / nVertLevelsMOD
+       hZLevel(i) =  h_total_max / nVertLevelsMOD
    end do
 
    h_sNew(:) = -h_total_max

Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/overflow/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/overflow/basin_src/basin-template.F        2012-02-07 20:43:31 UTC (rev 1474)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/overflow/basin_src/basin-template.F        2012-02-07 22:21:26 UTC (rev 1475)
@@ -316,8 +316,8 @@
  ! setting for three levels - Set h values for isopycnal system
    write(6,*) ' setting three levels for isopycnal system'
    do i = 1, nVertLevelsMOD
-       hNew(1,i,:) = i * h_total_max / nVertLevelsMOD
-       hZLevel(i) =  i * h_total_max / nVertLevelsMOD
+       hNew(1,i,:) = h_total_max / nVertLevelsMOD
+       hZLevel(i) =  h_total_max / nVertLevelsMOD
    end do
 
    h_sNew(:) = -h_total_max

</font>
</pre>