<p><b>mpetersen@lanl.gov</b> 2012-07-26 07:12:23 -0600 (Thu, 26 Jul 2012)</p><p>BRANCH COMMIT: merge trunk to branch<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/surface_pressure_mrp/src/core_ocean/mpas_ocn_time_integration_rk4.F
===================================================================
--- branches/ocean_projects/surface_pressure_mrp/src/core_ocean/mpas_ocn_time_integration_rk4.F        2012-07-26 00:08:47 UTC (rev 2056)
+++ branches/ocean_projects/surface_pressure_mrp/src/core_ocean/mpas_ocn_time_integration_rk4.F        2012-07-26 13:12:23 UTC (rev 2057)
@@ -258,6 +258,21 @@
!
call mpas_timer_start("RK4-cleaup phase")
+ block => domain % blocklist
+ do while (associated(block))
+
+ ! In main RK loop, tracer array is phi*h. Now divide by h^{n+1} to obtain phi^{n+1}.
+ do iCell=1,block % mesh % nCells
+ do k=1,block % mesh % maxLevelCell % array(iCell)
+ block % state % time_levs(2) % state % tracers % array(:,k,iCell) &
+ = block % state % time_levs(2) % state % tracers % array(:,k,iCell) &
+ /block % state % time_levs(2) % state % h % array(k,iCell)
+ end do
+ end do
+
+ block => block % next
+ end do
+
if (config_implicit_vertical_mix) then
call mpas_timer_start("RK4-implicit vert mix")
block => domain % blocklist
</font>
</pre>