<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) &
= block % state % time_levs(oldBtrSubcycleTime) % state % uBtrSubcycle % array(iEdge) &
+ dt/config_n_btr_subcycles *( &
@@ -806,6 +811,9 @@
- block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(cell1) ) &
/block % mesh % dcEdge % array(iEdge) &
+ 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>