<p><b>mpetersen@lanl.gov</b> 2011-10-26 08:18:48 -0600 (Wed, 26 Oct 2011)</p><p>Update new branch with bug fixes from the last week: Correct boundary update for RK4 KE and split explicit tracers.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/split_explicit_mrp/src/core_ocean/mpas_ocn_time_integration_rk4.F
===================================================================
--- branches/ocean_projects/split_explicit_mrp/src/core_ocean/mpas_ocn_time_integration_rk4.F        2011-10-25 21:01:33 UTC (rev 1136)
+++ branches/ocean_projects/split_explicit_mrp/src/core_ocean/mpas_ocn_time_integration_rk4.F        2011-10-26 14:18:48 UTC (rev 1137)
@@ -148,10 +148,10 @@
                                             block % parinfo % edgesToSend, block % parinfo % edgesToRecv)
 
            if (config_h_mom_eddy_visc4 &gt; 0.0) then
-              call mpas_dmpar_exch_halo_field2d_real(domain % dminfo, block % state % time_levs(2) % state % divergence % array(:,:), &amp;
+              call mpas_dmpar_exch_halo_field2d_real(domain % dminfo, provis % divergence % array(:,:), &amp;
                                                block % mesh % nVertLevels, block % mesh % nCells, &amp;
                                                block % parinfo % cellsToSend, block % parinfo % cellsToRecv)
-              call mpas_dmpar_exch_halo_field2d_real(domain % dminfo, block % state % time_levs(2) % state % vorticity % array(:,:), &amp;
+              call mpas_dmpar_exch_halo_field2d_real(domain % dminfo, provis % vorticity % array(:,:), &amp;
                                                block % mesh % nVertLevels, block % mesh % nVertices, &amp;
                                                block % parinfo % verticesToSend, block % parinfo % verticesToRecv)
            end if

Modified: branches/ocean_projects/split_explicit_mrp/src/core_ocean/mpas_ocn_time_integration_split.F
===================================================================
--- branches/ocean_projects/split_explicit_mrp/src/core_ocean/mpas_ocn_time_integration_split.F        2011-10-25 21:01:33 UTC (rev 1136)
+++ branches/ocean_projects/split_explicit_mrp/src/core_ocean/mpas_ocn_time_integration_split.F        2011-10-26 14:18:48 UTC (rev 1137)
@@ -963,20 +963,18 @@
            block =&gt; block % next
          end do
 
-        ! ---  update halos for prognostic variables
+        ! update halo for thicknes for unsplit only
+        if (trim(config_time_integration) == 'unsplit_explicit') then
+           block =&gt; domain % blocklist
+           do while (associated(block))
+              call mpas_dmpar_exch_halo_field2d_real(domain % dminfo, block % tend % h % array(:,:), &amp;
+                                               block % mesh % nVertLevels, block % mesh % nCells, &amp;
+                                               block % parinfo % cellsToSend, block % parinfo % cellsToRecv)
+              block =&gt; block % next
+           end do
+        endif ! unsplit_explicit
 
-        block =&gt; domain % blocklist
-        do while (associated(block))
-           call mpas_dmpar_exch_halo_field2d_real(domain % dminfo, block % tend % h % array(:,:), &amp;
-                                            block % mesh % nVertLevels, block % mesh % nCells, &amp;
-                                            block % parinfo % cellsToSend, block % parinfo % cellsToRecv)
-           call mpas_dmpar_exch_halo_field3d_real(domain % dminfo, block % tend % tracers % array(:,:,:), &amp;
-                                            block % tend % num_tracers, block % mesh % nVertLevels, block % mesh % nCells, &amp;
-                                            block % parinfo % cellsToSend, block % parinfo % cellsToRecv)
-           block =&gt; block % next
-        end do
 
-
         block =&gt; domain % blocklist
         do while (associated(block))
            allocate(hNew(block % mesh % nVertLevels))
@@ -1056,11 +1054,6 @@
              + block % state % time_levs(2) % state % uBcl % array(:,iEdge) 
            end do ! iEdge
 
-         ! mrp 110512  I really only need this to compute h_edge, density, pressure.
-         ! I can par this down later.
-         call ocn_diagnostic_solve(dt, block % state % time_levs(2) % state, block % mesh)
- 
-
          elseif (split_explicit_step == config_n_ts_iter) then
 
            hNew(:) = block % mesh % hZLevel % array(:)
@@ -1086,6 +1079,30 @@
          block =&gt; block % next
        end do
 
+        ! Boundary update on tracers.  This is placed here, rather than 
+        ! on tend % tracers as in RK4, because I needed to update 
+        ! afterwards for the del4 diffusion operator.
+        block =&gt; domain % blocklist
+        do while (associated(block))
+           call mpas_dmpar_exch_halo_field3d_real(domain % dminfo, block % state % time_levs(2) % state % tracers % array(:,:,:), &amp;
+                                            block % tend % num_tracers, block % mesh % nVertLevels, block % mesh % nCells, &amp;
+                                            block % parinfo % cellsToSend, block % parinfo % cellsToRecv)
+           block =&gt; block % next
+        end do
+
+
+         if (split_explicit_step &lt; config_n_ts_iter) then
+            ! mrp 110512  I really only need this to compute h_edge, density, pressure.
+            ! I can par this down later.
+            block =&gt; domain % blocklist
+            do while (associated(block))
+
+               call ocn_diagnostic_solve(dt, block % state % time_levs(2) % state, block % mesh)
+
+               block =&gt; block % next
+            end do
+         endif
+
       end do  ! split_explicit_step = 1, config_n_ts_iter
       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
       ! END large iteration loop 

</font>
</pre>