<p><b>mpetersen@lanl.gov</b> 2011-10-11 14:41:41 -0600 (Tue, 11 Oct 2011)</p><p>Branch commit, timesplitting_mrp, fix bug in u initialization on ocean-land boundaries.  This caused incorrect vertical velocities that led to lack of tracer conservation with topography.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_time_integration.F
===================================================================
--- branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_time_integration.F        2011-10-11 19:54:22 UTC (rev 1060)
+++ branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_time_integration.F        2011-10-11 20:41:41 UTC (rev 1061)
@@ -1144,7 +1144,12 @@
           if (block % mesh % boundaryEdge % array(1,iEdge).eq.1) then
               block % state % time_levs(2) % state % u % array(:,iEdge) = 0.0
           else
-              block % state % time_levs(2) % state % u % array(:,iEdge) = uTemp(:) + uCorr
+            do k=1,block % mesh % maxLevelEdgeTop % array(iEdge)
+              block % state % time_levs(2) % state % u % array(k,iEdge) = uTemp(k) + uCorr
+            enddo
+            do k=block % mesh % maxLevelEdgeTop % array(iEdge)+1,block % mesh % nVertLevels
+              block % state % time_levs(2) % state % u % array(k,iEdge) = 0.0
+            enddo
           endif
 
          ! Put new sshEdge values in h_edge array, for the compute_scalar_tend call below.

</font>
</pre>