<p><b>mpetersen@lanl.gov</b> 2011-06-06 12:53:28 -0600 (Mon, 06 Jun 2011)</p><p>Enforce boundary condition of u=0 within each barotropic subcycle.<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-06-05 17:31:26 UTC (rev 876)
+++ branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_time_integration.F        2011-06-06 18:53:28 UTC (rev 877)
@@ -797,6 +797,11 @@
 !     uBtrSubCycleNew(iEdge) = uBtrSubcycleOld(iEdge) + dt*(fEdge(iEdge)*uBtrPerp(iEdge) 
 !         - grad_ssh + GBtrForcing(iEdge))
 
+          ! mrp 110606 efficiency note: could make this a 1D integer factor instead of an if statement.
+          if (block % mesh % boundaryEdge % array(1,iEdge).eq.1) then
+                block % state % time_levs(newBtrSubcycleTime) % state % uBtrSubcycle % array(iEdge) = 0
+          else
+
                 block % state % time_levs(newBtrSubcycleTime) % state % uBtrSubcycle % array(iEdge) &amp; 
               = block % state % time_levs(oldBtrSubcycleTime) % state % uBtrSubcycle % array(iEdge) &amp; 
               + dt/config_n_btr_subcycles *( &amp;
@@ -806,6 +811,9 @@
                           - block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(cell1) ) &amp;
                           /block % mesh % dcEdge % array(iEdge) &amp;
                       + block % state % time_levs(1) % state % GBtrForcing % array(iEdge))
+
+          endif
+
 !print *, '9.4'
          end do
 !print *, '9.5'
@@ -962,8 +970,8 @@
               enddo
 
 ! mrp temp 110602 testing split: remove correction term.
-             uCorr = (block % state % time_levs(1) % state % FBtr % array(iEdge) - uhSum)/hSum
-! uCorr = 0.0
+!             uCorr = (block % state % time_levs(1) % state % FBtr % array(iEdge) - uhSum)/hSum
+ uCorr = 0.0
 ! mrp temp 110602 testing split end
 
               ! put u^{tr}, the velocity for tracer transport, in uNew

</font>
</pre>