<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("ts_rk4_stage_tend_hu_u_vertTransport")
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("ts_rk4_stage_tend_hu_u_vertTransport")
!
! velocity tendency: pressure gradient
!
+ call timer_start("ts_rk4_stage_tend_hu_u_pgrad")
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("ts_rk4_stage_tend_hu_u_pgrad")
!
! 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("ts_rk4_stage_tend_hu_u_dissipation")
if ( config_h_mom_eddy_visc2 > 0.0 ) then
do iEdge=1,grid % nEdgesSolve
cell1 = cellsOnEdge(1,iEdge)
@@ -604,10 +609,12 @@
deallocate(delsq_vorticity)
end if
+ call timer_stop("ts_rk4_stage_tend_hu_u_dissipation")
!
! velocity tendency: nonlinear Coriolis term and grad of kinetic energy
!
+ call timer_start("ts_rk4_stage_tend_hu_u_coriolis")
do iEdge=1,grid % nEdgesSolve
cell1 = cellsOnEdge(1,iEdge)
cell2 = cellsOnEdge(2,iEdge)
@@ -626,6 +633,7 @@
end do
end do
+ call timer_stop("ts_rk4_stage_tend_hu_u_coriolis")
!
! 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("ts_rk4_stage_tend_hu_u_forcing")
do iEdge=1,grid % nEdgesSolve
! forcing in top layer only
@@ -664,10 +673,12 @@
endif
enddo
+ call timer_stop("ts_rk4_stage_tend_hu_u_forcing")
!
! velocity tendency: vertical mixing d/dz( nu_v du/dz))
!
+ call timer_start("ts_rk4_stage_tend_hu_u_vmix")
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("ts_rk4_stage_tend_hu_u_vmix")
call timer_stop("ts_rk4_stage_tend_hu_u")
end subroutine compute_tend
</font>
</pre>