<p><b>ringler@lanl.gov</b> 2010-12-09 16:14:21 -0700 (Thu, 09 Dec 2010)</p><p><br>
final level of granularity.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/profiling/src/core_ocean/module_time_integration.F
===================================================================
--- branches/ocean_projects/profiling/src/core_ocean/module_time_integration.F        2010-12-09 22:24:01 UTC (rev 641)
+++ branches/ocean_projects/profiling/src/core_ocean/module_time_integration.F        2010-12-09 23:14:21 UTC (rev 642)
@@ -432,6 +432,7 @@
       !
       ! velocity tendency: vertical advection term -w du/dz
       !
+      call timer_start(&quot;ts_rk4_stage_tend_hu_u_vertTransport&quot;)
       if (config_vert_grid_type.eq.'zlevel') then
         allocate(w_dudzTopEdge(nVertLevels+1))
         w_dudzTopEdge(1) = 0.0
@@ -456,10 +457,12 @@
         enddo
         deallocate(w_dudzTopEdge)
       endif
+      call timer_stop(&quot;ts_rk4_stage_tend_hu_u_vertTransport&quot;)
 
       !
       ! velocity tendency: pressure gradient
       !
+      call timer_start(&quot;ts_rk4_stage_tend_hu_u_pgrad&quot;)
       rho0Inv = 1.0/config_rho0
       if (config_vert_grid_type.eq.'isopycnal') then
         do iEdge=1,grid % nEdgesSolve
@@ -481,12 +484,14 @@
           end do
         enddo
       endif
+      call timer_stop(&quot;ts_rk4_stage_tend_hu_u_pgrad&quot;)
 
       !
       ! velocity tendency: del2 dissipation, </font>
<font color="black">u_2 </font>
<font color="black">abla^2 u
       !   computed as </font>
<font color="black">u( </font>
<font color="black">abla divergence + k \times </font>
<font color="gray">abla vorticity )
       !   strictly only valid for config_h_mom_eddy_visc2 == constant
       !
+      call timer_start(&quot;ts_rk4_stage_tend_hu_u_dissipation&quot;)
       if ( config_h_mom_eddy_visc2 &gt; 0.0 ) then
          do iEdge=1,grid % nEdgesSolve
             cell1 = cellsOnEdge(1,iEdge)
@@ -604,10 +609,12 @@
          deallocate(delsq_vorticity)
 
       end if
+      call timer_stop(&quot;ts_rk4_stage_tend_hu_u_dissipation&quot;)
 
       !
       ! velocity tendency: nonlinear Coriolis term and grad of kinetic energy
       !
+      call timer_start(&quot;ts_rk4_stage_tend_hu_u_coriolis&quot;)
       do iEdge=1,grid % nEdgesSolve
          cell1 = cellsOnEdge(1,iEdge)
          cell2 = cellsOnEdge(2,iEdge)
@@ -626,6 +633,7 @@
 
          end do
       end do
+      call timer_stop(&quot;ts_rk4_stage_tend_hu_u_coriolis&quot;)
 
       !
       ! velocity tendency: forcing and bottom drag
@@ -633,6 +641,7 @@
       ! mrp 101115 note: in order to include flux boundary conditions, we will need to 
       ! know the bottom edge with nonzero velocity and place the drag there.
 
+      call timer_start(&quot;ts_rk4_stage_tend_hu_u_forcing&quot;)
       do iEdge=1,grid % nEdgesSolve
 
         ! forcing in top layer only
@@ -664,10 +673,12 @@
         endif
 
       enddo
+      call timer_stop(&quot;ts_rk4_stage_tend_hu_u_forcing&quot;)
 
       !
       ! velocity tendency: vertical mixing d/dz( nu_v du/dz))
       !
+      call timer_start(&quot;ts_rk4_stage_tend_hu_u_vmix&quot;)
       allocate(vertViscTop(nVertLevels+1))
       if (config_vert_visc_type.eq.'const') then
         vertViscTop = config_vert_viscosity
@@ -708,6 +719,7 @@
 
       end do
       deallocate(fluxVertTop, vertViscTop)
+      call timer_stop(&quot;ts_rk4_stage_tend_hu_u_vmix&quot;)
       call timer_stop(&quot;ts_rk4_stage_tend_hu_u&quot;)
 
    end subroutine compute_tend

</font>
</pre>