<p><b>mpetersen@lanl.gov</b> 2011-10-17 12:57:45 -0600 (Mon, 17 Oct 2011)</p><p>A little clean-up<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/split_explicit_timestepping/src/core_ocean/mpas_ocn_mpas_core.F
===================================================================
--- branches/ocean_projects/split_explicit_timestepping/src/core_ocean/mpas_ocn_mpas_core.F        2011-10-17 18:29:40 UTC (rev 1096)
+++ branches/ocean_projects/split_explicit_timestepping/src/core_ocean/mpas_ocn_mpas_core.F        2011-10-17 18:57:45 UTC (rev 1097)
@@ -228,7 +228,7 @@
 ! mrp changed to 0
 !              :block % mesh % nVertLevels,iCell) =  -1e34
 
-! mrp 110516 temp, added just to test for conservation of tracers
+! mrp 110516, added just to test for conservation of tracers
          block % state % time_levs(1) % state % tracers % array(3,:,iCell) = 1.0
 
       end do

Modified: branches/ocean_projects/split_explicit_timestepping/src/core_ocean/mpas_ocn_time_integration_split.F
===================================================================
--- branches/ocean_projects/split_explicit_timestepping/src/core_ocean/mpas_ocn_time_integration_split.F        2011-10-17 18:29:40 UTC (rev 1096)
+++ branches/ocean_projects/split_explicit_timestepping/src/core_ocean/mpas_ocn_time_integration_split.F        2011-10-17 18:57:45 UTC (rev 1097)
@@ -159,7 +159,7 @@
 
         block =&gt; domain % blocklist
         do while (associated(block))
-! mrp 110512 not sure if I need the following three.  Leave be, assume I need it.
+
            call dmpar_exch_halo_field2dReal(domain % dminfo, block % state % time_levs(2) % state % pv_edge % array(:,:), &amp;
                                             block % mesh % nVertLevels, block % mesh % nEdges, &amp;
                                             block % parinfo % edgesToSend, block % parinfo % edgesToRecv)
@@ -611,7 +611,6 @@
           else
 
              ! uBtrNew = uBtrOld + dt*(-f*uBtroldPerp - g*grad(SSH) + G)
-! mrp 110926 change to weighted avg of SSH.  revised to:
 
              sshCell1 = &amp;
                (1-config_btr_gam2_SSHWt1)*block % state % time_levs(oldBtrSubcycleTime) % state % sshSubcycle % array(cell1) &amp;
@@ -633,20 +632,6 @@
                       + block % state % time_levs(1) % state % u_diffusionBtr % array(iEdge))
                       ! added del2 diffusion to btr solve
 
-! mrp 110926 change to weighted avg of SSH.  orig:
-!                block % state % time_levs(newBtrSubcycleTime) % state % uBtrSubcycle % array(iEdge) &amp; 
-!              = block % state % time_levs(oldBtrSubcycleTime) % state % uBtrSubcycle % array(iEdge) &amp; 
-!              + dt/config_n_btr_subcycles *( &amp;
-!                        uPerp &amp;
-!                      - gravity &amp;
-!                        *(  block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(cell2) &amp;
-!                          - block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(cell1) ) &amp;
-!                          /block % mesh % dcEdge % array(iEdge) &amp;
-!                      + block % state % time_levs(1) % state % GBtrForcing % array(iEdge) &amp;
-!                      + block % state % time_levs(1) % state % u_diffusionBtr % array(iEdge))
-                      ! added del2 diffusion to btr solve
-! mrp 110926 change to weighted avg of SSH.  end
-
           endif
 
          end do
@@ -685,7 +670,6 @@
 
        end do !do BtrCorIter=1,config_n_btr_cor_iter
 
-! mrp 110922 for extra corrector eta at end
             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
             ! Barotropic subcycle: Compute thickness flux and new SSH: CORRECTOR
             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -753,16 +737,8 @@
                block =&gt; block % next
             end do  ! block
 
-! mrp 110923 print
-!print *, 'sshsubcycleOld 2',minval(block % state % time_levs(oldBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve)), &amp;
-!                            maxval(block % state % time_levs(oldBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve))
-!print *, 'sshSubcycleNew 2',minval(block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve)), &amp;
-!                            maxval(block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve))
-! mrp 110923 print end
-
         endif ! config_btr_solve_SSH2
 
-
             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
             ! Barotropic subcycle: Accumulate running sums, advance timestep pointers
             !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -884,12 +860,6 @@
                 * block % tend % ssh % array(iCell) / block % mesh % areaCell % array (iCell)
          end do
        endif
-         ! Now can compare sshSubcycleNEW (big step using summed fluxes) with
-         ! sshSubcycleOLD (individual steps to get there)
-!print *, 'ssh, by substeps',minval(block % state % time_levs(oldBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve)), &amp;
-!                            maxval(block % state % time_levs(oldBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve))
-!print *, 'ssh, by 1 step  ',minval(block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve)), &amp;
-!                            maxval(block % state % time_levs(newBtrSubcycleTime) % state % sshSubcycle % array(1:block % mesh % nCellsSolve))
 
          ! Correction velocity    uCorr = (Flux - Sum(h u*))/H
          ! or, for the full latex version:

</font>
</pre>