<p><b>dwj07@fsu.edu</b> 2012-10-05 09:28:04 -0600 (Fri, 05 Oct 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Merging trunk to branch.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_monthly_forcing.F
===================================================================
--- branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_monthly_forcing.F        2012-10-05 00:08:19 UTC (rev 2196)
+++ branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_monthly_forcing.F        2012-10-05 15:28:04 UTC (rev 2197)
@@ -176,6 +176,8 @@
 
       if(config_use_monthly_forcing) then
         monthlyForcingOn = .true.
+
+        write (0,'(a)') &quot; Monthly forcing is on.  Make sure monthly forcing variables include iro in Registry, and are in your initial condition or restart file.&quot;
       end if
 
    !--------------------------------------------------------------------

Modified: branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_mpas_core.F
===================================================================
--- branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_mpas_core.F        2012-10-05 00:08:19 UTC (rev 2196)
+++ branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_mpas_core.F        2012-10-05 15:28:04 UTC (rev 2197)
@@ -116,27 +116,27 @@
 
       call ocn_init_split_timestep(domain)
 
-      print *, ' Vertical grid type is: ',config_vert_grid_type
+      write (0,'(a,a10)'), ' Vertical grid type is: ',config_vert_grid_type
 
       if (config_vert_grid_type.ne.'isopycnal'.and. &amp;
           config_vert_grid_type.ne.'zlevel'.and. &amp;
           config_vert_grid_type.ne.'zstar1'.and. &amp;
           config_vert_grid_type.ne.'zstar'.and. &amp;
           config_vert_grid_type.ne.'zstarWeights') then
-         print *, ' Incorrect choice of config_vert_grid_type.'
+         write (0,*) ' Incorrect choice of config_vert_grid_type.'
          call mpas_dmpar_abort(dminfo)
       endif
 
-      print *, ' Pressure type is: ',config_pressure_type
+      write (0,'(a,a10)'), ' Pressure type is: ',config_pressure_type
       if (config_pressure_type.ne.'pressure'.and. &amp;
           config_pressure_type.ne.'MontgomeryPotential') then
-         print *, ' Incorrect choice of config_pressure_type.'
+         write (0,*) ' Incorrect choice of config_pressure_type.'
          call mpas_dmpar_abort(dminfo)
       endif
 
       if (config_filter_btr_mode.and. &amp;
           config_vert_grid_type.ne.'zlevel')then
-         print *, 'filter_btr_mode has only been tested with'// &amp;
+         write (0,*) 'filter_btr_mode has only been tested with'// &amp;
             ' config_vert_grid_type=zlevel.'
          call mpas_dmpar_abort(dminfo)
       endif

Modified: branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_time_integration_rk4.F
===================================================================
--- branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_time_integration_rk4.F        2012-10-05 00:08:19 UTC (rev 2196)
+++ branches/omp_blocks/openmp_test/src/core_ocean_elements/mpas_ocn_time_integration_rk4.F        2012-10-05 15:28:04 UTC (rev 2197)
@@ -236,9 +236,8 @@
            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;
-                                                                        + rk_weights(rk_step) * block % tend % tracers % array(:,k,iCell) )&amp;
-                                                                       / block % state % time_levs(2) % state % h % array(k, iCell)
+                                                                        block % state % time_levs(2) % state % tracers % array(:,k,iCell) &amp;
+                                                                        + rk_weights(rk_step) * block % tend % tracers % array(:,k,iCell)
               end do
            end do
 
@@ -256,6 +255,20 @@
       !  A little clean up at the end: decouple new scalar fields and compute diagnostics for new state
       !
       call mpas_timer_start(&quot;RK4-cleaup phase&quot;)
+
+      ! Rescale tracers
+      block =&gt; domain % blocklist
+      do while(associated(block))
+        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) &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>