<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(&quot;RK4-cleaup phase&quot;)
 
+      block =&gt; 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) &amp;
+               = block % state % time_levs(2) % state % tracers % array(:,k,iCell) &amp;
+                /block % state % time_levs(2) % state % h % array(k,iCell)
+            end do
+         end do
+
+         block =&gt; block % next
+      end do
+
       if (config_implicit_vertical_mix) then
         call mpas_timer_start(&quot;RK4-implicit vert mix&quot;)
         block =&gt; domain % blocklist

</font>
</pre>