<p><b>mpetersen@lanl.gov</b> 2011-06-07 17:04:39 -0600 (Tue, 07 Jun 2011)</p><p>Small change to fix bug I introduced into Higdon unsplit h value.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_mpas_core.F
===================================================================
--- branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_mpas_core.F        2011-06-07 20:36:07 UTC (rev 881)
+++ branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_mpas_core.F        2011-06-07 23:04:39 UTC (rev 882)
@@ -83,8 +83,9 @@
       call compute_mesh_scaling(mesh)
  
       call rbfInterp_initialize(mesh)
-      call init_reconstruct(mesh)
-      call reconstruct(block % state % time_levs(1) % state, mesh)
+! mrp temp
+!      call init_reconstruct(mesh)
+!      call reconstruct(block % state % time_levs(1) % state, mesh)
 
       ! initialize velocities and tracers on land to be -1e34
       ! The reconstructed velocity on land will have values not exactly

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-07 20:36:07 UTC (rev 881)
+++ branches/ocean_projects/timesplitting_mrp/src/core_ocean/module_time_integration.F        2011-06-07 23:04:39 UTC (rev 882)
@@ -346,7 +346,8 @@
 
          call compute_solve_diagnostics(dt, block % state % time_levs(2) % state, block % mesh)
 
-         call reconstruct(block % state % time_levs(2) % state, block % mesh)
+! mrp temp
+!         call reconstruct(block % state % time_levs(2) % state, block % mesh)
 
          block =&gt; block % next
       end do
@@ -1079,11 +1080,19 @@
 
       if (trim(config_time_integration) == 'higdon_unsplit') then
 
+         do iCell=1,block % mesh % nCells
            ! this is h_{n+1}
-             block % state % time_levs(2) % state % h % array(:,:) &amp;
-           = block % state % time_levs(1) % state % h % array(:,:) &amp;
-           + dt* block % tend % h % array(:,:) 
+             block % state % time_levs(2) % state % h % array(:,iCell) &amp;
+           = block % state % time_levs(1) % state % h % array(:,iCell) &amp;
+           + dt* block % tend % h % array(:,iCell) 
 
+            ! this is only for the hNew computation below, so there is the correct
+            ! value in the ssh variable for higdon_unsplit case.
+            block % state % time_levs(oldBtrSubcycleTime) % state % sshSubcycle % array(iCell) &amp;
+          = block % state % time_levs(2) % state % h % array(1,iCell) &amp;
+          - block % mesh % hZLevel % array(1)
+           end do ! iCell
+
       endif ! higdon_unsplit
 
 ! printing:

</font>
</pre>